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ABSTRACT 

The  results  reported  are  part  of  a  continuing  study 
to  improve  numerical  models  for  meso-scale  and  small-scale 
effects  which  influence  global  weather  and  its  modification. 
The  two  major  areas  being  studied  are  the  effects  of  mountain 
ranges  on  energy  and  momentum  transfer  and  the  transient 
interaction  of  solar  radiadon  with  the  Earth's  atmosphere. 

In  the  study  of  orographic  effects  on  global  climate 
two  new  codes  to  investigate  the  effects  of  compressibility 
and  moisture  are  described  and  calculations  are  reported. 
Several  modifications  made  to  the  general  Boussinesq  code 
described  in  the  last  semiannual  are  outlined  and  results  of 
test  problems  reported.  Additional  computations  completed 
with  the  meso-scale  codes  are  shown  including  a  comparison 
of  dalculational  results  and  experimental  data  for  flow  over 
the  Sierra  Nevada  Range. 

A  new  theoretical  and  numerical  scheme  for  solving 
the  radiative  transfer  equation  in  the  Earth's  atmosphere  has 
been  developed.  It  involves  the  approximation  of  the  trans¬ 
mission  function  by  a  sum  of  exponentials,  some  plausible 
assumptions  about  frequency-averaging,  based  on  the  idea  of 
Lobcsguc  quadrature,  and  the  use  of  the  (’.rant  and  Hunt 
algorithm  to  solve  a  formally  equivalent  set  of  monochromatic 
problems  within  each  frequency  group.  In  addition  to  being 
able  to  handle  simultaneous  band  absorption  and  scattering, 
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this  method  allows  the  use  of  zones  of  arbitrary  size  and  is 
non-iterative.  Difficulties  with  numerical  ill-conditioning 
and  non -uniqueness  in  exponential -  sum  approximation  are  de¬ 
scribed  and  a  resolution  of  these  difficulties  proposed. 

Some  improvements  in  the  treatment  of  Mie  scattering  are 
presented,  and  the  overall  status  of  the  radiation  code  is 
reviewed. 
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NOMENCLATURE 


phase  speed 

specific  heat  at  constant  pressure 
sound  speed 

drag  force  on  the  obstacle 


fluid  vorticity  =  —  - 

saturation  vapor  pressure 
advective  flux  across  a  boundary 
Coriolis  parameter 
acceleration  of  gravity 
dry  adiabatic  lapse  rate  =  g/C 
numerical  grid  indices 
horizontal  wave  number 


R/Cp 

latent  heat  of  vaporization  for  water,  or 
a  characteristic  distance 

cloud  water  content 

rain  water  content 

ratio  of  molecular  weight  of  water  and  dry 

water  production  terms 

pressure 

Helmholtz  pressure  function 


3SR-1034 


NOMENCLATURE,  contd. 

q  3  total  water  content  excluding  rainwater 

(j>  3  compressibility  stream  function  defined  in 

liq.  (3.6),  or  angle  between  mountain  range 
and  northern  direction  (Section  4) 

4'  3  stream  function 

ft  3  rotational  velocity  of  Earth 

p  3  defined  by  Eq.  (5.7a) 

R  3  gas  constant  for  air 

Ry  3  gas  constant  for  water  vapor 

r  3  relative  humidity 

r  3  saturated  water  vapor  mixing  ration 

p  3  density 

S  3  static  stability 

s  3  density  stratification 

T  3  temperature 

t  3  time 

0  3  potential  temperature  latitude  (Section  4) 

Vt  3  terminal  velocity  of  water  droplet  in  atmosphere 

/•  A 

U  3  iu  +  kw  =  total  two-dimensional  velocity 
u  3  horizontal  velocity  component  (x- direction) 
v  3  horizontal  velocity  component  (y- direction) 
w  3  vertical  velocity 
x  3  horizontal  Cartesian  coordinate 
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NOMENCLATURE,  contd. 


I 

!...  I 


y  =  horizontal  Cartesian  coordinate 
z  =  vertical  Cartesian  coordinate 
c;  =  compressibility  vorticity  function 

-  IxtP-J  - 


SUBSCRIPTS,  SUPERSCRIPTS 

o  =  value  associated  with  the  static  atmosphere 
'  =  perturbation  quantity 

i,j  =  numerical  grid  indices 
g  =  geostrophic  state 


i 
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NOMENCLATURE  for  Section  6 


P 

T 

z 


T’Ti 


n 


S 

e,<P 

v 


m 


sun 
k  . 


a. 


pressure 

temperature 

vertical  coordinate  (measured  from  top  of 
atmosphere) 

total  optical  depth 

optical  depth  of  zone  boundary  n 

number  of  zones 

scattering  optical  depth 

absorption  optical  depth 

frequency 
unit  vector 

spherical  angular  coordinates  of  & 
cos  6 

one  of  the  set  of  Gaussian  or  Radau  quadrature 
angles 

number  of  angles 

cosine  of  scattering  angle  0g 

cosine  of  solar  zenith  angle 

volume  extinction  coefficient  =  a'  +  B 

v  v 

volume  absorption  coefficient  =a(^ne)  +  (cont.) 

v  v 
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NOMliNCLATURH,  contd. 


a 


(line) 


a 


v 

(cont) 


a' 


v 

h 


v  ,M 


^R 


a 


kCO 

v 


=  line  absorption  coefficient 

=  continuum  absorption  coefficient 

ri  -hv/kT) 

*  (1  -  e  Ja^ 

=  Planck's  constant 

=  volume  scattering  coefficient  =3  ..  +  3 

v,M  Mv,R 

=  Mie  scattering  coefficient 
=  Rayleigh  scattering  coefficient 


v  ^  =  Mie  absorption  coefficient 


smooth  part  of  extinction  coefficient 
c'Ccont)  ♦  6 

V 


=  solar  flux 


v 


Pv  =  scattering  phase  function 


P  =  azimuthal  average  of  P 


v  ,M 
Tv 

j-solar 

v 


=  Mie  phase  function 


azimuthal  average  of  specific  intensity  of 
radiation  I 

v 


=  solar-beam  part  of  I 


v 


diffuse  intensity  (I  with  solar  beam  sub¬ 
tracted  out)  v 

solar-beam  scattering  source  for  diffuse  inten¬ 
sity 
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NOMENCLATURE,  contd. 


P 


v 


a 


1  >a2 


uv 


ai 


k. 


i 

M 


a. 


Uj 


N 


u 

+ 


azimuthal  average  of  the  surface  bidirectional 
reflectivity 

surface  directional  emissivity 
surface  temperature 
transmission  function 


factors  in  scaling  approximation 
(%  “  “1(v)a2(p,T)) 

equivalent  absorber  amount  for  scaling  approxi¬ 
mation 

«Jline>  (P0»T0),  where  pQ  and  Tq  are  STP 

coefficients  in  exponential-sum  fit  to  T^v 

exponents  in  exponential-sum  fit  of 

number  of  terms  in  exponential - sum  fit  to 

absorption  coefficient  corresponding  to 

one  of  a  discrete  set  of  u-’s  at  which  T.  is 
to  be  fitted  by  an  exponen- 
tial  sum 

number  of  u . ’ s 
J 


u"  =  positively-  and  negatively-directed  intensity  vectors 
u*  - 

r,rij»rn»rn  n±1  =  reflection  matrices 

t,tij,tn,tn  n±l  =  transmission  matrices 
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NOMliNCLATURI:,  contd. 


Z±»Zij*En»En,n+l»i:(i,h)  s  source  vectors 
I  =  identity  matrix 

(d  =  albedo  for  single  scaLiering  *  3V/<V 
*  one  of  Gaussian  or  Radau  weights 


Ax 


max 


n 


maximum  allowable  primary- layer  thickness 
-2nAx/u 

9  rgnn 

c  ,  where  Ax  is  the  primary- layer  thick¬ 

ness  in  a  zone 

quantities  involved  in  the  lincar-in-x  inter¬ 
polation  of  » 


Yn  =  intermediate  quantities  involved  in  source 
doubling  for  the  Planck  source 

3  mxni  matriccs  computed  in  the  forward  pass  of  the 
Grant  and  Hunt  algorithm  and  used  on  the  backward 
pass  to  calculate  intensities 

/  *1  \  /  M  \  r  «•  N 

Vn  *vn  ;,nV,i  }  "  ■-vectors  computed  in  the  forward  pass  of 

the  Grant  avid  Hunt  algorithm  and  used  on  the 
backward  pass  to  calculate  intensities 

rc  c  reflection  matrix 


w  =  vector  connected  with  the  surface  boundary  condition 
Fy  *  radiation  flux 

a  =  radius  of  spherical  aerosol  particle 
X  *  wavelength 
a  =  2na/A 

n(a)  =  distribution  of  aerosol  radii 
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XOMI-XCUniBIi,  contd. 


(AO)  «  fundamental  jnigular  interval  in  tabulation  of 
0  Mic  phase  function 

0*  *  angle  below  which  Mic  phase  function  is  trun¬ 

cated 

p*  3  truncated  Mic  phase  function 

*  adjusted  Mic  scattering  coefficient  and  phase 
function  after  truncation 


:jk’cj'cj 


(1) 


c(2) 

ck 


various  quantities  used  in  re-normalizing 
the  phase  function 


xiv 


3SR-1034 


1,  INTRODUCTION 

The  results  reported  herein  are  the  continuation  of 
the  development  of  numerical  codes  which  are  being  utilized 
in  the  study  of  meso-scale  phenomena  related  to  flow  over 
mountain  ranges  and  the  interaction  of  solar  radiation  with 
the  Earth’s  atmosphere. 

1.1  OROGRAPHIC  EFFECTS  ON  GLOBAL  CLIMATE 

The  scope  of  work  during  the  past  six  months  study 
has  emphasized  a  continuing  effort  to  develop  and  expand  nu¬ 
merical.  codes  which  may  be  used  to  understand  the  physical 
processes  which  influence  momentum  transfer.  The  last  semi¬ 
annual  report ^  gave  results  of  flow  over  an  obstacle  ar¬ 
rived  at  by  using  a  numerical  code  based  on  the  Boussincsq 
approximation.  Some  indications  of  codes  which  we  hoped  to 
develop  to  study  compressible  effects  and  moisture  effects 
were  outlined.  This  report  will  further  describe  the  devel¬ 
opment  of  those  codes  and  present  some  test  problems  run  with 
two  new  codes.  In  addition,  a  problem  based  on  experimental 
data  measured  in  the  Owens  Valley  area  of  California  lias  been 
completed  using  our  Boussincsq  code  for  flow  over  the  Sierra 
Nevada  Range.  Comparisons  made  between  the  experimental  and 
calculated  results  show  good  agreement. 

Further  developments  which  have  been  completed  during 
the  past  six  months  include  a  scheme  to  take  into  account 
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variable  grid  sizes  in  the  vertical  direction,  triangular 
zones  which  allow  a  better  description  of  the  obstacles 
representing  mountain  ranges,  a  linear  analysis  of  the  equa¬ 
tions  used  in  the  codes  developed  to  date  with  comparisons 
with  the  complete  compressible  and  incompressible  equations, 
and  an  outline  of  a  code  to  be  developed  which  takes  into  ac¬ 
count  Coriolis  forces. 

1 . 2  HAD I AT I VO  TRANSFER 

The  atmospheric  radiation  effort  has  been  concentrated 
on  a  revamping  of  the  method  of  solution  proposed  in  the  semi¬ 
annual  report  in  order  to  take  advantage  of  a  sophisticated 
new  method  of  Grant  and  lhmt^0^  for  monochromatic  radiative 
transfer  problems.  Grant  and  Hunt's  method  is  the  most  ad¬ 
vanced  available  for  plane-parallel  problems,  and  has  two 
large  advantages  over  other  techniques;  (1)  it  imposes  no 
restrictions  on  zone  size,  except  that  there  will  be  an  in¬ 
evitable  loss  of  accuracy  if  atmospheric  structure  (p,  T, 
etc.)  is  specified  too  coarsely;  (2)  it  is  non-iterative, 
whether  or  not  there  is  scattering.  Since  the  technique  is 
only  applicable  to  monochromatic  problems,  however,  it  is  not 
directly  useful  unless  one  is  willing  to  do  a  line-by-line 
calculation  in  band  absorption  regions.  Barring  this,  some 
approximation  is  necessary.  The  one  chosen  involves  fitting 
the  transmission  function  for  a  frequency  group  &v  (con¬ 
taining  many  lines)  with  a  sum  of  M  exponentials, 

TAv<">  = 

i»l 


M 


a^  c 


-kjU 


(1.1) 
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Then  !»>•  jmoklng  certain  plausible  31s  umpiions  based  on  the 
idea  of  Lebesgtto  quadra  lure,  the  frequency*  averaged  form?  of 
the  radiative  transfer  equation  may  he  broken  down  into  hi 
monochromatic  problems,  each  of  which  may  he  solved  by  ftrant 
and  Hunt's!  method.  The  technique  just  described,  if  It  proves 
accurate  when  compared  with  line*t*y* line  calculation*,,  will  he 
very  powerful  indeed,  and  will  solve  seme  long*# (land Inn  dlffi* 
cutties  which  have  clouded  all)  atmospheric  heating* rate  eaten* 
lotion#  up  to  now.  Tor  example,  It  will  improve  calculation* 
In  all  infrared  absorption  band#  in  which  scattering  cannot 
be  neglected. 


Code  development  has  proceeded  apace  with  the 
isent  of  the  theory.  Further  improvement*  have  beea  **ide  in 
the  treatment  of  *lie  scattering,  and  excessively  large  forwai 
peaks  In  the  phase  function  are  new  truncated  and  replaced 
by  pure  forward  scattering.  (Great  difficulties  were  encoMit 
cred  in  numerical  ly  imp!  emeu  ting  the  exponential  fit  of  iq. 
(1.1),  due  to  numerical  1 1l*condlf lotting  m& 
but  the  resolution  of  that  problem  is  now-  in  oig.1 
the  straightforward  subroutine*  have  been  ceded, 
those  for  calculating  (I)  sene  structure,  angle  structure, 
and  frequency  group  structures  {2}  mgtftMpftlcal  parameters 
such  as  the  solar  flux  S  ,  solar  s«&iiti  angle,  arl  tattle 
Sun  distance;  (%}  ttayleigh  scattering  parameters;  (4)  the 
frequency  average  of  the  flame!  function!  f$>)  transmissie# 
functions  from  tabular  data;  (b|  tbs.  index  of  refract $m  #( 
water  at  any  frequency,  from  extensive  tables,  for  use  in 


Hie  scattering  and  water  surface  reflect 
and  (?)  the  *>20  distribution  of  Hie  scat »t rets. 
Input  or  from  one  of  a  number  of  analytic  forms 
the  literature.  I  hi  (Irani  and  giant 
subroutines  are  partially  complete. 
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fklUng  *b«ul«l  hr  rc4*lf  «wn.  ©*!il  loss 

routine  Mill  eot^ttfic  tbs*  £©4*?,  aitssl  c«mf>arisi«iM$  with  flood* 

ard  tto^rlelfb-scauering)  foluiien*  astd  witf*  8b©  Hims- 
AhUjm  radio!  logs  |»arls»ft<f  nu*'  iLd  tog  in. 
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2*  omsmic  effects?  mi  uitmvm  mm  wm 


irt.i  imst®  ct*H$tlct<n4  «!«*rii$§  Ik?  ini* 
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I'  I®  stud*  i*#  «44it i«a&l  test  mw 

u  line  tfttst  <$#®c rife-*-#  flow  met  m  etoteete 

«ii*»  mi  iiaHtiliiy  eel  tiiiriseetftt  wind 
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feMlttS  $f#  i 


1*1 


l!a#» 


TO.f  tftfflr 


%**  «*llHg  t 

at  IS«  i  !#• 


Helens* 


site  #i 


m  tt* 


§ 


They  are  reiterated  in  Table  2.3.  The  Initial  horizontal 
velocity  profile  is  ghowr.  in  Figure  2.1, 


TABU;  2.1 

PROM.LH 

TI1I.I1 

LAPSE 

run.* 

ca  i  it  character  i  st  i  t:s 

HORIZONTAL 

VELOCITY 

PROFILE 

two  Wave 

T  *  if 

&\  *  1500  neters 

As  *  441  neters 

Obstacle: 

Height  *  441  neters 
Length  *  4S00  neters 

Exponential 
profile  as 
shown  in 
figure  2.1 

Initial  surface  temperature  *  5tHJ*K 
f  *  adiabatic  lapse  rate  *  1©#L/Las 


Ihe  % inear  analyst*  *f  p«in  and  l*eldrUJ2l  indicates 
tl>«t  unde?  these  conditions  Kw  Mares  should  he  present  in 
of  the  mountain  -  one  wave  of  approximately  9.2  Li* 

of  approx  I  stately  2S  In  wave* 

&  cwItfMlai  Sen  was  r m  to  a  tin©  ©f  49St» 

[s  of  the  *t ream  I ines  and  the  vertical 
the  calculation  at  a  tine  ©f  it*9i 
ii  **2.  She  shortet  wave  appears 

te  obstacle  displaying  a  wavelength 
<jt  second  wav©  appears  behind  the 
f  f  to  #  lit  with  a  wavelength  of  $p« 
ftt to  and  tOldtf'il>*s  analysis  indicated 
*f  she  mm  new  Id  appear  at  7.4  to. 
the  mim  features  of  the  I  i dour  theory  ate  ©b* 

I hose  waves  ate 
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ligiirc  2.1.  Initial  velocity  profile  -  two  wave  problem. 
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distinctly  indicated  in  the  vertical  velocity  contour  plots 
also  shown  in  the  figure. 

2.1.1  Wave  Drag  Results 

In  stratified  flow  over  an  obstacle  the  pressure  is 
systematically  higher  on  the  upstream  side,  resulting  in  a 
drag  force  on  the  obstacle,  and  a  corresponding  drag  of  op¬ 
posite  sign  on  the  air  stream. 

This  drag  force  of  the  air  on  the  mountain  is  ob¬ 
tained  by  integrating  the  momentum  transport  along  the 
x-axis.  For  a  symmetrical  mountain  of  height  H, 

/•  00 

Drag  =  -  J  (pu'w*)  |Z=H  dx  , 


where  u',w’  =  velocity  perturbation  quantities. 

An  edit  routine  was  written  to  integrate  the  results 
obtained  in  these  calculations.  The  results  are  presented 
as 


pu'w ' 


pu '  w  * 


dx 


where  L  is  one-half  the  computational  x- interval.  The  drag 
units  puTw*r  ,  as  presented  here,  arc  in  dynes/cm2.  The 
quantity  L  is  48  km.  The  values  presented  may  be  related 
to  other  studies  where  drag  force  is  presented  in  dyncs/cm  by 
multiplying  the  values  shown  here  by  21.. 

The  momentum  edits  pu'w'  ,  located  one  cell  above 
the  mountain  top,  arc  shown  as  a  function  of  time  in  Figure 
2.3.  The  largest  value  of  the  drag  which  was  reached  at 
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1000  seconds,  was  approximately  equal  to  4  dynes/cm2.  This 
value  agrees  qualitatively  with  measured  values  of  the  mo - 
mentum  flux  reported  by  D.K.  Lilly1-  J  for  measurements  at 
Boulder,  Colorado  (7  dynes/cm2).  Palm  and  Foldvik,  and 
Vergeiner^  present  calculated  values  of  7*106  to  2><107 
dynes/cm  for  similar  problems.  The  value  of  4  dynes/cm2 
corresponds  to  3.8*107  dynes/cm  and  thus  we  also  see  quali¬ 
tative  agreement  here. 

The  other  important  feature  of  the  momentum  flux 
edit  is  the  oscillatory  character  of  the  values  with  time. 

It  is  thought  to  be  related  to  the  formation  of  the  individ¬ 
ual  vertical  velocity  cells,  i.e.,  as  new  positive  or  nega¬ 
tive  cells  are  formed,  there  is  an  increase  or  decrease  in 
the  horizontal  average  of  the  vertical  flux  of  horizontal 
momentum. 

Figure  2.4  presents  the  two-wave  momentum  flux  as  a 
function  of  height  at  several  times.  The  momentum  edits  in¬ 
dicate  a  cyclic  character  at  late  times  at  a  height  of  3  to 
5  km.  The  interaction  of  the  long  and  short  waves  takes  place 
in  this  altitude  range  and  is  thought  to  be  responsible  for 
this  phenomenon. 

2.2  SIERRA  NEVADA  WAVE  PROBLEM 

A  test  problem  using  the  HAIFA  code  was  formulated 
based  on  experimental  data  (see  Figure  2.5)  gathered  from 
flights  over  the  Owens  Valley  on  February  16  ,  1952.^  These 
data  included  measurements  taken  on  two  sailplane  flights 
which  extended  over  a  time  period  of  3-1/2  hours.  An  alti¬ 
tude  of  33,000  feet  was  attained.  The  data  are  especially 
interesting  because  of  its  completeness  and  the  presence  of 
a  strong  lee  wave. 
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Figure  2.4.  Momentum  flux  as  a  function  of 
height  for  two-wave  problem. 
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The  numerical  calculation  was  made  using  constant 
zone  sizes  with  Ax  =  2.0  km  and  Az  =  0.5  km.  The  total 
height  of  the  grid  was  17.5  km  while  the  horizontal  length 
was  128  km.  The  horizontal  velocity  was  treated  as  a  con¬ 
stant  above  9.5  km  and  the  temperature  was  assumed  constant 
above  14  km.  The  Sierra  Nevada  range  was  treated  as  a  tri¬ 
angular  barrier  with  a  front  face  assumed  perpendicular  to 
the  Earth's  surface,  with  a  maximum  height  of  2  km  and  a 
horizontal  length  of  8  km. 

Figure  2.6  gives  the  measured  streamline  pattern 
directly  in  the  lee  of  the  mountain.  The  wavelength  of  the 
major  lee  wave  above  the  mountain  is  21  to  22  km.  A  rotor 
was  also  observed  in  the  Owens  Valley  in  the  lee  of  the 
mountain  range. 

The  results  of  our  calculation,  to  a  time  of  approxi¬ 
mately  2000  seconds,  are  shown  in  Figures  2.7  and  2.8.  The 
results  indicate  a  small  rotor  was  formed  in  the  lee  of  the 
obstacle  and  the  calculated  lee  wavelength  is  22  km.  These 
comparisons  arc  especially  significant  as  the  problem  pre¬ 
sents  a  very  severe  test  of  the  HAIFA  code.  The  temperature 
profile  is  not  smooth  and,  in  fact,  a  small  inversion  layer 
exists  at  an  altitude  of  20  km.  A  large  portion  of  the  hori¬ 
zontal  velocity  profile  is  treated  as  uniform  flow  and,  as 
noted  from  our  previous  results, this  type  of  flow  excites 
a  continuous  spectrum  of  waves.  Also,  numerical  stability  as¬ 
sociated  with  uniform  flows  adds  to  the  difficulty  of  the 
computer  calculation.  A  direct  comparison  with  the  experi¬ 
mental  streamline  results  shown  in  Figure  2.6  is  presented 
in  Figure  2.9  at  a  time  of  2000  seconds.  The  presence  of  the 
two  small  rotors  at  that  time,  in  contrast  to  the  larger 
single  rotor  observed  in  the  experiment,  may  be  a  consequence 
of  the  fact  that  the  calculation  has  not  yet  reached  a  steady 
state. 
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5.  CODE  TiODIFICATIOriS 


TVo  najor  code  mod i fi cat ions  have  been  completed  in 
the  past  six  months.  A  code  similar  to  HAIFA  but  which  in¬ 
cludes  compressibility  of  the  atmosphere  was  completed  and 
the  two  wave  test  problem  was  run  with  it,  and  compared  with 
the  Koussi itesq  flow  calculations.  A  second  code  which  in¬ 
cludes  moisture  effects  was  also  completed,  and  a  test  cal¬ 
culation  nade.  both  of  these  codes  are  discussed  in  detail 
in  this  section. 

In  addition  to  the  two  major  developments  noted, 
several  ©the*1  general  improvements  in  the  codes  were  com¬ 
pleted.  These  included  modifications  to  allow  the  use  of 
variable  grid  sires  in  the  vertical  direction,  the  use  of 
triangular  rones  to  improve  the  description  of  mountain 
ranges,  a  test  of  the  Crowley ^ 1 aJveetion  scheme  and  a 
isodi ficalion  to  improve  its  accuracy,  and  filially  an  edit 
routine  to  calculate  pressure  throughout  the  flow  field, 
particularly  for  the  compressible  code. 

I.ach  of  these  modifications  is  also  described  in  this 
section.  In  addition,  a  listing  of  the  compressible  code  is 
Included  as  Appendix  f*  of  this  report. 
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3.1  COMPRESSIBILITY  EFFECTS 

Modifications  to  the  basic  HAIFA  equations  to  take 
account  of  the  compressibility  of  the  atmosphere  have  been 
previously  reported . ^ ^ ^  The  resulting  system  of  compressi¬ 
bility  equations  (3.1  through  3.4)  have  been  put  in  finite 
difference  form,  and  a  computer  code  has  been  written  which 
incorporates  these  in  a  two-dimensional ,  time-dependent  cal¬ 
culation. 

The  compressibility  equations  are: 


3pu  3pu2 
3t  3x 


3pw  .  3puw  ^ 
3T"  +  3x  + 


3puw  3p 
3z  3x 


ipwl  +  lE  = 

3z  3z 


gP 


(3.1) 

(3.2) 


3pu  3pw 
3x  3z 


(3.3) 


3pT'  3pT 1 u  3  pT 1 w 

3t  3x  3z 


(3.4) 


Here  we  recall  the  assumption  that  the  density  at 
every  position  can  be  determined  from  the  perfect  gas  equa¬ 
tion  of  state,  in  which  the  pressure  takes  the  value  associ¬ 
ated  with  the  static  atmosphere,  pQ  ,  through  the  relation 


P  = 


Po 

RT 


(3.5) 


The  system  of  equations  (3.1  through  3.4)  have  a  form 
similar  to  the  Boussinesq  equations  of  standard  HAIFA,  and 
can  be  solved  in  a  similar  manner. 
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A  stream- function- like  quantity  <p  can  be  introduced 
such  that  Bq.  (3.3)  is  satisfied  and 

Pu  =  §f  ,  pw  =  -  ||  .  (3.6) 

If  a  vorticity-like  function  c  is  defined, 

-  !p_u  3pw  (x  7-, 

^  “  9z  "  9x  ^  J 

then  the  same  Poisson  equation  as  for  the  Boussinesq  approxi¬ 
mation  results, 


C 


3  2 

Tx~ 


•f 


9  2  (J) 
9Z* 


(3.8) 


The  prognostic  equation  for  C  is  obtained  by  cross  differen 
tiating  Bqs .  (5.1)  and  (3.2)  and  subtracting, 


H 

9t 


9x 


(uO 


9 

9z 


9_ 

~9z 


'  9  (|>  9w 
1 9x"  9x" 


(V.C) 

+ 

9  /9(f)  9u  +  9 (p  9u\ 
9x\9x  9x  Tz  Tzf 

9w' 

1  =  c3p  -  _  SP 

9T ' 

9  z 

9  Z ) 

f  *»9x  T  +  T' 

n 

9x 

(3.9) 


Bq.  (3.9)  replaces  the  vorticity  equation  of  the  Boussinesq 
system  of  equations,  differing  principally  in  having  the 
additional  terms  containing  the  derivatives  of  <J> ,  u,  and  w. 
These  additional  terms,  which  for  convenience  we  call  the 
chi  terms,  will  be  discussed  in  detail  in  a  following  section 
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3.1.1  Method  of  Numerical  Solution 

3. 1.1.1  Caleulati.onal  Logic  —  The  sequence  of  cal¬ 
culations  for  one  cycle  of  the  compressible  code  is: 


SUBROUTINE 

FUNCTION 

UPDATE 

CD 

Eq.  (3.4)  is  solved  for 
new  values  of  (pT') 

(2) 

Eq.  (3.9)  is  solved  for 

C 

(3) 

As  adjuncts,  p  is  computed 
from 

o=o-  ^T') 

M  Mo  T 

o 

and  T'  follows  from 

T»  _  CpT') 
p 

LAPLAC 

(4) 

Eq .  (3.8)  is  solved  for 
new  values  of  <j> 

VELOC 

(5) 

From  Eq.  (3.6) 

u  is  derived  from 

u  =  H /p 

(6) 

v  is  derived  from 

V  *  -  %  /P 

3.1. 1.2  Finite  Difference  Scheme  -  The  centering  of 
the  dependent  variables  in  the  compressible  code  is  identical 
to  the  centering  of  the  analogous  variables  in  standard 
HAIFA,  viz., 
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The  difference  scheme  used  in  the  standard  HAIFA 
equations  is  applied  dire  ;tly  to  the  compressible  system  of 
equations,  and  need  not  be  reviewed  here.  The  differencing 
of  the  chi  term,  which  is  unique  to  the  compressible  code, 
is  discussed  below. 

To  form 


3H  + 

_9£  9u*l 

.  3  p<!> 

3w 

9<J) 

9wl 

[dx 

3x 

'dZ  3Z  J 

9  z  [_9x 

9x 

9  z 

9  z  J 

we  define  X...  X.  .  Z...  and 


vm  iom 


that 


4  (tel 


3.1*2  till  I  Sal  sm* 


<3*  will*  IIAII'A*  llir  talur*  ef  si#|f|  *t$4  I  |$|  are 
*JKPCif3^«l  toy  iisp'ii I  t*»  ttor  c@4r.  in  a»«ta3 1 1  i «phm 9 

the  initial  pnptSM'tt  #w**t  to#  ^prclfisfj. 

fwnwi iiing  initial  $r***m%$  at#  fmmiS  m*i«g  tfe#  ^tasint 


i##s 


tl.tij 


m>  that 


^  4§  A/ 

*  Pe(i^)  W|*|*  |  j  t 

Tilt*  initial  ilcii*it)'  firoTili*  (lien  l**  1 1  trv  *-  Train 

'V'-l  *  iuou)  • 

The  s  t  re ;m*  f unc  t  i  on  -  I  i  he  quan  l  i  i  v  4  1 $  fo  rue J  by  i  n  t  ej:  r:i  t  * 

ing 


H 

dz 


R 


0o(-.) 


w  ( z) 

U  * 
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and  the  verticil} ‘MU*  quantity  t,  is  initialised  fro®  the 
appli cat  ion  of  Cq.  (3.7)  to  yield 

4  *  hK<£)  «„(=))  • 

Matidard  HAIFA,  initially  there  is  no  x* dependence  of 
any  qua®  lily* 


1% 


SMifUal  low  of  Advert  Ion  Scholar 

quantities  i@  he  adverted  In.  t! 

are  t  ,  li|.  (3.9), 
nitMiSon  of  continuity 


system  of 
Cot*)  » 
is  th#  for 


# 


It  'v-'i*  ftrtfif' 


Svt  m&'iff  tiwK toy ** 

to  n*o 


pt  -§tm 

ffoolo* 


I  he  no  Jar  cnnytitationul  difference  hetfeeen  the 
hoMss i r.esq  and  compressible  codes  Is  the  so  Uniat  ion  of  the 
chi  term*  llisinjg  studies  of  the  two  code  versions  mere  per* 
lotted,  based  on  »  grid  of  61  x*di  reel  ion  cells  and  3$ 
y-direction  cells,  and  incorporating  an  interior  boundary 
condition  (uhich  requires  two  passes  per  cyajle  in  the  Poisson 
iqu.it  ion  Solver),  The  results  are  sttsana  r  i  zed  iti  the  table 
be  1 ow : 
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COMI'KHSSHMi.l: 

STANDARD 

UI’DATI;  TOTAL 

2,,t*  Order  Scheme 

Chi  Terra 

1.S8  sec 

0.33  sec 

1 .02  sec 

0.51  sec 

0.33  sec 

LA PL AC  TOTAL 

1.14  sec 

1.14  sec 

VI1L0C  TOTAL 

0.14  sec 

i.ll  sec 

TINSTI*  101 AL 

0.03  sec 

o.t>3  &#c 

»asi t  umijoii.h 

“  -  ■  —  -  ■■  1  —  fc— ■  ^ 

2. $$  sec 

■k.  *  — n  ii  —if  fftriiir'--  r-  -  r  -yi 

m€ 

tmzs^.>,*ami;is&!a-  ssw-sr.-.., 

iitr  a*. !  i * .  i  %*>$*§«*»  #f  iiit  ill 8 til 


(•Iri. 

twijiarai lv*  piPts  tHguw*  $J  anJ  3.2)  *sf  the  rcsiilts 
from  I  lit  coupr  !<  ible  i*«*4«*  and  ill! PA  i«l  m  sltwr  of  jipjifVAiifean^* 
1>*  IS9Q  h*vo»i!.«  iiitlli'nu^  «  single  tfaw  (in  each  east)  wis.li  n 
10*12  In  wavelength.  Hie  vcriical  velocity  pattern  indicates 
that  only  the  single  wave  Is  forming  in  tin?  compressible  cause 
while  a  second  wave  is  appearing  hi  the  boussinesq  case. 
However,  since  t lie  compressible  calculation  was  run  to  a  tine 
of  only  2000  seconds,  these  results  are  inconclusive.  The 
first  wavelength  in  each  case  agrees  with  the  shorter  of  the 
two  waves  predicted  by  Palm  and  1‘oldvik. 
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3. 1.5. 2  Xunerical  Stability  —  The  stability  of  the  comp res  - 
sihle  finite  difference  equations  is  not  well  understood  at 
this  tine.  Originally,  the  line-step  criterion  of  the 
Rintss  i  nes»|  f ©rani  1  a t  i  on , 


M 


<  Odin 


*****  I ©»’©d  f#i  the  cciiprrsfiitle  e«|ua  siting,  In  tsetse  In* 
ft  incrf *  calrii1»ti«tiiiil  in*sal  i ti sir*  developed*  IjicI*. S«»^  n 
stubmil*  I*  illicit  in  lutce  *tev  tim 

•fltp  tfltfi  lot •  §  *%;i l|-  r^ftpf tarte  |^rf*%r40:dl  In  aw 

<#l  l  #*»$*!  I#  iuMithl  tiie  #iat4©  rmge  far  tin?  cf^ihrf&SMe  r#4*:, 

A  iMMiiif*  *  «  |j» 


M 


■td4ii 


»»*•«*»  #*« 


iMtvtMhMd §4  m  m  jttitwiqpt  t©  gtapllfy  ik  l*wt#i8ji»ti©#* 
lt§i$  a^nfpiels  did  *1**5  yield  tiM-ful  refill  t$9  a*  m  prered  t© 
tu?  feiufely  pretties*  depeodc?it •  Indeed*  value*  a*  M§tll  j» 
n  *  0.2  have  fit  I ©4  l«  ptmlde  a  fffltite  Ci*npfv*s*ll*t©  feint  lest 
lit  $«IC  €•**#*• 

■*  I  stirprifingly ,  the  i«aj«»r  *4*14 1  i  If  pretoletif  have 
been  eitferved  In  the  regi®hf  ©f  highest  fpeed  flov,  vhieh 
typically  are  the  upperswst  port  I*****  of  the  ealculat §©nal  grid 
f©r  tuwy  a  si)*®*  pit  ef  le  tv?!  problem,  ihi*  if  true,  for  ex* 
an  pie ,  ©f  the  tv©  cave  test  j»  rob  least  ttf  thif  feet  Son*  lurther, 
the  effects  ©f  nest  interest  in  cede  e«upar»f©nf  for  the  tv© 
wave  problem  tale  place  near  the  obstacle,  far  removed  fret* 
the  upper  atwsphere.  litis*  fugges*  ted  another  approach;  the 
in t redoes ion  of  s 1 1 ©ng  artificial  diffusion  in  the  upper  cells 
of  the  grid,  forcing  the  faint  i©i»  there  to  becorn*  spatial  ly 
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smoothed.  While  the  use  of  diffusion  proved  successful  in 
this  application,  not  enough  is  known  yet  about  the  stability 
characteristics  of  the  compressible  formulation,  further 
theoretical  and  numerical  research  is  warranted  in  this  area, 
both  to  improve  the  accuracy  and  economics  of  the  compressibl 
formulation,  and  to  advance  understanding  of  atmospheric 
tics© -scale  codes. 

3.2  hoistuki.  i.h  u ts 

Atmospheric  water  in  the  fern  of  water  vapor,  cloud 
water,  and  j.revipH.i?  lea  any  have  Important  effects  on  the 
district  eristic*  of  gravity  waves  canned  by  n;oon  mitts,  bee 
«*mses  are  frepitttl ly  aww^aiiled  by  clouds  which  can  be  ex¬ 
pected  to  modify  the  stnMlIt)  of  the  air  through  the  pres* 
e*icc  of  the  latent  heat  «tf  condensation  which  the  cloud  water 
add*-  so  she  air.  t'owse«|swt9S ly ,  she  terms  re?*ttl sing  in 
chatty;?*  of  stability  ©f  the  air  in  which  clouds  are  forcing 
arc  «f  primary  interest  • 

Hie  development  «f  $  cade  to  determine  She  effects  of 
moisture  on  she  ions  for  a  hotissinesq  fluid  is  essential 

ly  complete*  the  differential  equal  Sons  describing  iho  physi 
cal  problem  have  been  reported  previously .  ^  Hie  HAN  A 

equations  bate  been  modified  to  Incorporate  the  following 
changes : 

(I)  she  lament  tm  equation  incorporates  she 
effects  of  mils  sure  In  the  buoyancy 
ter**, 

(21  the  equation  of  stale  for  air  tabes 
account  of  a  w4aicr  vapor  component , 
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(3)  the  energy  equation  includes  energy 
changes  equivalent  to  the  latent  heat 
of  water  being  given  to  or  taken  from 
the  air, 

(4)  an  equation  is  added  to  account  for 
the  conservation  of  all  moisture  ex¬ 
cept  min  water,  and 

(5)  a  conservation  equation  is  added  which 
governs  the  rain  water  content  in  the 
atmosphere  including  sources  and  sinks 
at  the  boundaries. 


The  differential  equations  are  outlined  beloK.  All  terms  arc 
defined  in  the  nomenclature  list. 


(5.12) 


9T" 

at 


,DT"  .  3T" 

i -  +  - 


as 


-wr 


(3.1 3) 


ill  ♦  u2l!  ♦  wiJl  ,  .  JSL-n  ♦  f  +  £  lil! 
at  uax  \n  tqu  c  Vax 


(3.14) 


♦  £-  1 
O  LP 


3r 


.  !M 

Max  ax  /  » 


tr  c(J  *  ;'c ♦  +  « 


a;.,.  avt.  £  v.  .  ov 

_ L  .  p  _  J  _» 1  op  . 

at  r  3x  ■  P  ax  r  dz 


t 


£  ,  V.  . 

r  t  Op 

P  r)Z 


(3.15) 


Mr  3£  r3£. 

+  Ut — -  ♦  wr--  -  V 
ax  oz 


p: .  !!ri . 
tlax  3z  J  1  r  » 


39 


3SK-1034 


t 


+ 


3x 


+  w 


M  = 

3z 


-A 


r 


Pr 
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where  pr  =  source  terms  for  the  conversion  of  water  sub¬ 
stance  to  rain  water.  These  source  terms  are  described  in 
detail  in  a  following  section  of  this  report. 


3.2.1  Integration  of  the  Finite  Difference  Equations 

The  Eqs.  (3.12)  through  (3.16)  are  written  in  finite 
difference  form  and  integrated  numerically  in  a  similar  man¬ 
ner  to  the  HAIFA  equations. 

The  basic  scheme  used  in  the  integration  is  shown  in 
Figure  3.7.  The  locations  of  the  major  variables  with  re¬ 
spect  to  the  grid  cells  are  shown  in  Figure  3.8. 


3. 2. 1.1  The  AdA'ection  Scheme  —  The  advection  terms  for  vortic 
ity,  water  vapor  content  and  rain  water  arc  calculated  using 
the  second  order  scheme  of  Crowley.  The  scheme  is  written  in 
conservation  form  and  is  based  on  a  forward  time  difference 
and  centered  space  differences.  Test  calculations  performed 
by  Crowley  indicated  that  for  the  same  order  of  accuracy, 
the  conservation  form  produced  more  accurate  solutions  than 
the  advection  form. 

In  the  conservation  form,  the  time  derivative  and  ad¬ 
vection  terms  of  the  vorticity,  temperature,  or  moisture 
equations  may  be  written  as 


3a  x  D(ua)  ,  3(va)  c 

3t  3x  3  z  b  » 

where  a  represents  T,  n,  q,  or  A  and  S  is  the  source 
term. 
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STEP  1  Initial  values  of  temperature,  vorticity, 

velocities  and  water  content  are  speci¬ 
fied  . 

STEP  2  Set  water  vapor  content  calculated  and 

initial  water  content  values  are  redis¬ 
tributed  as  rain  water,  clouds  and  vapor. 

STEP  3  New  values  of  temperature,  rain  water, 

water  content  other  than  rain  water,  and 
vorticity  arc  obtained  accounting  for  the 
advection  terms  only. 

STEP  4  Vorticity,  water  content  and  rain  water 

n  1 

are  updated  to  time  t  by  evaluating 
source  terms  with  the  quantities  obtained 
in  Step  3. 

STEP  5  The  Poisson  equation  is  solved  for  ip 

,  r  .  .n+1 

using  values  of  n  at  t 

STEP  6  Velocities  are  updated  using  \p  values 
from  Step  4.  Initial  data  for  a  new 
cycle  are  now  available  for  edit,  or  con¬ 
tinuation  of  the  calculation  beginning 
at  Step  2. 


Figure  3.7.  HAIFA '"MOISTURE  calculation  sequence 
used  in  numerical  integration  of 
Eqs .  (5.12)  through  (5.16). 
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In  two  dimensions  a  splitting  technique  is  used;  the 
calculational  scheme  calls  for  solving  two  one-dimensional 
equations  sequentially.  First,  the  net  flux  of  vorticity  or 
temperature  is  solved  for  in  the  horizontal,  the  quantity 
solved  for  in  the  zone  being  updated  due  to  this  flux.  The 
procedure  is  then  repeated  in  the  vertical  direction  using 
the  partially  updated  values  from  the  horizontal  calculation. 

The  temperature  equation  is  transformed  from  the  form 
shown  in  liq .  (3.13)  in  order  to  write  it  in  conservation  form. 
A  new  variable  3  is  defined,  at  T"  +  Tz  .  Since  r  is  at 
most  a  function  of  z  only,  Eq.  (3.13)  then  becomes 


3x 


+  w 


33 
3  z 


0 


In  this  manner,  only  the  advection  scheme  is  required 
to  update  the  temperature  since  all  source  terms  have  been 
removed.  Solving  the  temperature  equation  in  conservation 
form  should  increase  the  accuracy  of  the  solution. 

3 . 2 . 1 . 2  Changes  in  Vorticity  and  Water  Content  Due  to  Other 
Terms  in  the  Conservation  Equations  -  All  source  terms  in 
Eqs.  (3.14)  through  (3.10)  are  evaluated  as  zone  centered 
quantities  containing  centered  space  differences.  Update 
takes  place  after  the  effects  of  advection  have  been  calcu¬ 
lated  . 


3.2.2  Moisture  Equation  Source  Terms 

The  source  terms  in  the  conservation  equations  (3.12) 
through  (3.16)  depend  on  complicated  physical  processes  which 
have  been  parameterized  from  experimental  measurement  data. 
These  include  the  terminal  velocity  of  rain  drops,  the 
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equilibrium  vapor  pressure  of  water,  and  several  eon  stall  is 
dealing  with  the  water  production  term.  Parametric  represen¬ 
tations  of  these  terms  are  given  below;  we  have  attempted  to 
obtain  the  most  comprehensive  results  available  but  we  have 
not  yet  performed  critical  evaluations  of  these  terms. 

3. 2. 2.1  Terminal  Velocity  -  The  terminal  velocity  of  rain 
drops  is  given  by 


Vt  =  5.. 32 (tr) 1/5 


where  is  in  m/scc  and  in  gm/Kgin  . 

r  g  i 

This  expression,  derived  by  Liu  and  Orvill  c,OJ  has 
been  compared  to  data  from  the  Smithsonian  Meteorological 
Tables ^  and  an  expression  developed  by  Srivastava. 


3, 2. 2. 2  Saturation  Vapor  Pressure  -  The  saturation  vapor  pres 
sure  is  defined  as 


e 


s 


6.11  x  10' 


7.5(T0-2.73)/(To-36) 


ini  1  lib  ars 


where  TQ  =  Tq  -  gz/C  and  rg  ,  the  saturated  water  vapor 
mixing  ratio  in  gm/Kg,  is  given  by 


R  es(T0)  L(T«  -  L/C  ) 
r_  =  » — rr-r-rr  exp - jt — L~ 


Rv  *0™ 


o 


Both  expressions  arc  obtained  from  the  work  of  Orville  ,  8  *  ^ 
Ogura,^12^  and  Ogura  and  Phillips. 
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S.2.2.?  teaser  I'tiwlwct  ion  I  mu  «  Kin?  W'asier  j*ff»»I*is,s  i#-w  c*siisia&l(S 
«f  llirfe  physlul  i^awaa  which  can  add  10,  ssihkFssrt  fruM, 
or  change  I  hr  state  of  the  water  la  I  hr  at  •^sphere.  Ihis^  in* 
eludes  (I)  I  hr  evaperiat  ion  of  rail*  water  out wide  the 
(2)  the  conversion  of  eiussd  water  tu  rain  water,  s»§$i2  (3)  tin* 
growth  ©f  rain  through  coalescence.  those  f  enm>  were  orlgi* 
nnlly  expressed  in  Orville's  work  as 

l»r  *  -Mr-r^)  ♦  &UC'*C  1  *  l.o  *  I0‘*  lc<4-f)P‘"n  • 

During  otir  initial  t«.<t  calculations,  it  was  discovered  that 
treating  0  ,  the  evaporation  parameter,  as  a  constant  led  to 
computational  difficulties,  this  term  was  modified  so  agree 
with  that  used  in  Kessler's  work*  ‘  1  which  expressed  II  as 

8  *  3. *07  *  10*S(^?r)0*rS  . 

3. 2. 2. 4  Density  of  Air  Containing  tfater  Vapor  -  The  equation 
of  state  utilized  is  that  of  a  perfect  gas  consisting  of  air 
and  water  vapor.  The  expression  is 

*  *i>tU  ♦  r/n  I 
p  *  pUTj-j-— --j 

where  m  is  the  ratio  of  the  no locular  weight  of  water  and 
dry  air.  The  pressure  p  is  treated  as  p  (z),  the  initial 
atmospheric  distrihut .ion ,  while  T  and  r  are  allowed  to 
vary  in  both  the  horizontal  and  vertical  directions. 
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Figure  3.14.  Computer  plot  of  vertical  velocity  contours  -  Tuo 
wave  p rob  let:. 
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its  temperature  and  a  stronger  downdri f t  is  created.  These 
phenomena  effect  the  velocity  profiles  in  different  ways  at 
various  times. 

The  momentum  flux  results  shown  in  figure  5.J5  indi¬ 
cates  only  minor  differences  in  drag  between  the  moist  and 
dry  cases.  However,  since  the  problem  has  only  been  run  a 
short  time,  the  results  do  not  necessarily  reflect  those  ex¬ 
pected  near  steady  state. 

5c  5  HAIFA  trim  VAR i ABU:  VliRTl  CAL  ZONING 

The  modification  discussed  in  the  previous  semiannual 
(to  the  basic  HAIFA  code)  enabling  it  to  operate  using  a  grid 
mesh  of  variable  spacing  in  the  vertical  direction  was  com¬ 
pleted.  Three  test  problems  were  performed  using  a  625  meter 
high  mountain  similar  to  the  single  wave  problem  described  in 
Ref.  1.  The  initial  input  to  this  problem  was  (1)  a  tempera¬ 
ture  lapse  rate  equal  to  one -half  the  dry  adiabatic,  (2)  an 
obstacle  measuring  625  meters  high  by  1500  meters  Jong  and 
(5)  an  exponential  horizontal  velocity  profile  described  by 
the  equation 


u  »  9.04  exp  (1.795  *  10'4  z) 


Two  of  the  test  problems  used  constant  vertical 


grid  he  i  gilts 


ami  the  standard  HAIFA  code.  The.  grid  heights  were 
meters  and  208.5  meters.  The  results  were  compared 
third  problem  run  using  a  variable  vertical  grid  si 
grid  varied  from  125  meters  near  the  lower  boundary 
meters  at  the  upper  boundary.  The  total  number  of 
zones  in  all  three  calculations  was  the  same. 


512.5 
with  n 
ze .  Th  i  s 
to  758 
vertica  1 


for 


the 


Figures  5.16  and  5.17  contain  the  detailed  comparison 
Az  =  512.5  meter  problem  and  the  variable  Az  problem 
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of  the  temperature  and  vorticity  at  times  of  approximately 
1500  seconds.  Quantities  are  shown  as  a  function  of  z  for 
three  locations;  just  upstream  of  the  obstacle  (1=16),  just 
downstream  of  the  obstacle  (1  =  24)  ,  and  far  downstream  of  the 
obstacle ( 1=36) . 

It  can  be  noted  that  the  trends  in  the  results  are 
the  same  for  both  problems.  However,  some  differences  in  the 
computed  values,  particularly  at  small  heights,  are  large. 
Part  of  this  discrepancy  was  thought  to  be  due  to  the  smaller 
grid  near  the  lower  boundary.  As  a  result,  the  problem  was 
re-run  using  constant  Az  equal  to  208.3  meters.  The  results 
were  closer  to  the  variable  Az  case  indicating  that  the  grid 
size  was  playing  an  important  role  in  the  detailed  results. 
Figure  3.18  shows  the  momentum  flux  at  a  height  of  781  meters 
for  the  three  problems  as  a  function  of  time.  The  good 
agreement  between  the  small  constant  Az  case  and  the  variable 
Az  problem  should  be  noted.  In  particular,  the  cyclic  be¬ 
havior  of  the  drag  values  is  not  as  pronounced  as  in  the 
larger  zoned  case. 

These  phenomena  will  continue  to  be  investigated  as 
time  for  more  computations  becomes  available. 

3.4  TREATMENT  OF  TRIANGULAR  ZONES  IN  HAIFA 

In  order  to  characterize  mountain  slopes  by  zones 
other  than  of  rectangular  shape,  two  modifications  of  the 
HAIFA  code  were  required;  First,  cell  centered  variables 
were  assumed  to  be  located  at  the  centroid  of  the  zones  and 
Crowley's  advection  scheme  was  derived  fur  flow  between  tri¬ 
angular  and  rectangular  cells.  Second,  the  solution  of  the 
Poisson  equation  requires  the  value  of  the  vorticity  at  a 
grid  point  which  is  normally  obtained  from  an  arithmetic 
average  of  the  vorticity  of  the  four  surrounding  cells.  This 
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H.  CORIOLIS  TERiiS  Hi  2-D  iiESO-SCALE  EQUATIONS 


One  of  the  major  tanks  to  be  undertaken  in  the  next 
contract  period  will  he  to  consider  the  effect  of  Coriolis 
forces  on  the  flow  over  mountain  ranjj.es . 

In  the  2-D  calculations  previously  reported,  lee  waves 
are  formed  over  the  mountain  in  a  time  interval  of  about  one 
hour,  during  which  dynamic  effects  from  the  Coriolis  force  are 
small,  (The  geostrophic  wind  used  for  the  unperturbed  flow, 
of  course,  is  strongly  influenced  by  the  Coriolis  terms.)  If 
the  mountain  range  were  more  extensive  or  if  the  fate  of  the 
waves  radiated  by  the  mountain  were  followed  to  larger  dis¬ 
tances  it  is  to  be  expected  that  larger  effects  from  the 
Coriolis  terms  would  he  realized. 


Ihe  characteristic  distance  scale  for  the  Coriolis 
force  is  Lj.  =  u/f  ,  where  u  is  a  typical  wind  speed  and 
f  is  the  Coriolis  parameter  (if  u  ^  10  m/s  ,  we  obtain 
!*£>  ^  100  km).  When  the  mountain  range  is  comparable  to  L^., 
an  appreciable  modification  of  the  gravity  waves  will  result 
to  form  a  complex  system  of  gravity- inertia  waves.  An  inves¬ 
tigation  of  these  waves  when  they  are  small  in  amplitude  and 
the  wind  speed  is  uniform  has  been  carried  out  by  Qucney.^1^ 
Since  the  distance  l.(.  is  comparable  with  the  extent  of 
mountain  ranges  which  we  must  consider  in  climate  dynamics, 
it  appears  to  he  important  to  take  account  of  the  Coriolis 
force  in  our  subsequent  investigations. 
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Since  the  ('oriolis  force  induces  a  turning  of  the  wind 
it  is  necessary  to  take  into  account  several  new  factors  in 
the  calculations: 

(.1)  the  component  of  the  wind  parallel  to 
the  mountain  must  be  included; 

(2)  pressure  gradients  in  directions  par¬ 
allel  and  perpendicular  to  the  mountain 
must  be  included  to  establish  geo- 
strophic  balance  in  the  unperturbed 
flow;  and 

(3)  the  vertical  atmosphere  structure  is 
slightly  modified  to  account  for  the 
Coriolis  contribution  to  the  hydro¬ 
static  balance  condition. 

In  the  following  formulation  we  attempt  to  parallel  the  nu¬ 
merical  treatment  of  the  HAIFA  code  as  closely  as  feasible 
in  order  to  be  able  to  compare  the  effects  of  the  Coriolis 
terms  with  those  pertaining  to  a  non- rotating  Earth. 

4.1  FORMULATION 

The  differential  equations  of  the  dry  atmosphere  are 
formulated  in  a  system  of  reference  fixed  to  a  rotating 
Earth.  As  discussed  by  Thompson  wc  incorporate  the  cen¬ 

trifugal  terms  into  the  definition  of  the  local  gravity  to 
obtain 


dU 

dt 


+ 


2  SI  x  U  = 


-  VP 
P 


(4.1) 
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where  U  is  the  velocity  relative  to  the  Earth's  surface, 
ft  is  the  rotational  velocity  of  the  Earth,  P  is  the  pres¬ 
sure,  p  is  the  density  of  the  atmosphere,  g  is  the  local 
acceleration  of  gravity  which  is  assumed  to  act  in  the  verti¬ 
cal  direction  Ik  .  The  time  derivative  is  that  evaluated 
following  the  fluid  motion.  The  equations  of  mass  and  energy 
conservation  are  not  affected  by  the  Coriolis  force. 

Neglecting  the  curvature  of  the  Earth's  surface  (and 
the  resulting  centrifugal  terms  associated  with  the  relative- 
velocity)  ,  the  Eq .  (4.1)  can  be  resolved  into  components.  We 

choose  a  Cartesian  coordinate  system  in  which  the  x-axis  lies 
in  the  surface  and  forms  an  angle  <jJ  with  the  eastward  direc¬ 
tion,  the  y-axis  lies  in  the  surface  at  the  same  angle  c)> 
with  the  northward  direction,  and  the  z-axis  is  perpendicular 
positive  upward.  We  shall  subsequently  assume  that  the  x-axis 
is  perpendicular  to  the  2-D  mountain  range  which  is  oriented 
at  the  angle  <j>  with  the  northward  direction.  Denoting  x, 
y,  and  z  components  of  the  velocity  by  u,  v,  w,  the  compo¬ 
nent  equations  arc 

-  2 SI  sin0  v  +  2ft  cosG  coscj)  w  =  -  i  —  , 

+  2ft  sinO  u  -  2ft  cos 0  sine))  w  =  ~  |y  >  (4.2) 

-jY  ‘  m  cos0  (,l  C0-S(i’  "  v  s  1  n (;) )  =  -  ^  -  g  , 


where  ft  is  the  magnitude  of  the  rotational  velocity,  0 
is  the  latitude  of  the  position,  and  the  time  derivatives 
arc  those  formed  following  the  fluid  motion. 
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In  the  absence  of  tlic  mountain  barrier  we  consider 
the  atmosphere  to  be  in  gcostrophic  balance;  the  motion  :i  s 
unaccc  .Lerated ,  the  vertical  velocity  component  w  is  zero, 
and  the  pressure  gradients  arc  just  balanced  by  the  Coriolis 
and  gravity  terms.  Denoting  the  gcostrophic  state  by  sub¬ 
scripts  g  we  obtain 


-21  sin0  v  =  -  — 
g  Pg 

2  si  sin0  ii  =  -  — 

fi  Pn 


-21  cos6(u  cosef) 


vg  sllKj)) 


DP 

Dx  > 

DP 

£ 

Dy  ’ 

2 _ 8. 

-f-  p  = 

P„ 

(4.3) 


Since  tlie  velocity  of  the  unperturbed  state  is  taken  to  be 
independent  of  x  and  y  ,  Oq ,  (4.5)  shows  that  the  pressure 
at  most  need  depend  linearly  on  x  and  y  . 

When  the  mountain  is  present  the  pressure,  density 
and  velocity  arc  perturbed  from  their  gcostrophic  values.  It 
will  be  convenient  to  introduce  the  deviation  P '  of  the 
pressure  from  the  gcostrophic  value: 

P'  =  P  -  P  .  (4.4) 

6 


We  also  introduce  the  Boussinesq  approximation, 
departure  of  density  from  the  gcostrophic  value 
account  only  in  the  gravity  term  of  hq .  (4.5). 
terms,  we  use  the  gcostrophic  density: 


in  which  the 
is  taken  into 
In  all  other 
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-pjr  -  2 fi  sin0(v-v  )  +  2ft  cos0  coscj)  w  =  -  —  , 

Pg  X 

+  s inG (u-u  )  -  2ft  cos0  sin<J>  w  =  -  —  ,  (4.5) 

Pg  * 

*  2ft  cos0[(u-ug)cos<j)  -  (v-vg) sincj)]  =  '  ~  +  (l  -  ~~)g  . 

g  g 

These  equations  constitute  the  Boussincsq  approximation  for 
the  equations  of  motion  when  the  geos trophic  flow,  assumed 
steady  and  independent  of  x  and  y  ,  is  perturbed.  Clearly, 
they  are  only  approximately  satisfied  in  a  local  region,  since 
u  and  v  arc  not  constant  on  the  synoptic  scale. 

o  o 

We  now  consider  the  special  case  in  which  the  initial 
and  boundary  equations  arc  independent  of  the  y- coordinate , 
corresponding  to  a  uniform  but  obliquely  incident  wind  en¬ 
countering  a  two-dimensional  ridge,  the  topography  of  which  is 
independent,  of  y  .  In  this  case,  the  initial  and  boundary  con¬ 
ditions  and  the  equations  depend  only  on  the  coordinates  x 
and  z  .  The  resulting  2-D  equations  derived  from  Eq.  (4.5) 
are : 


9u 

at 


+ 


u 


,9u 

9x 


+  w 


9u 
9  z 


2ft 


sinO  (v-v 


+ 


2ft  cos0  cost))  w 


1  9P ' 

Pg  8x  • 


gY  +  u-—  +  w|y  +  s.inO(u-u  )  -  2ft  cosO  sin<J)  w  =  0  , 

A  o 

(4.6) 

9w  9w  9w  ~r.  „  ,  ,  ~  ,  . 

dt  U9x  W9'z  "  c°s0  [  (u-u  )cos<J>  -  (vv)sin  +  j  - 
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Eqs.  (4.6)  are  to  he  supplemented  with  the  equations 
of  incompressibility  and  the  temperature  equation.  The  equa¬ 
tions  are  unchanged  from  those  in  Ref.  1;  they  arc  to  be 
found  in  that  document  as  Eqs.  (2.3)  and  (2.10). 

There  arc  two  major  modifications  of  the  equations 
which  have  resulted  from  the  treatment  of  the  rotation  of  the 
Earth : 


(1)  The  component  of  the  wind  parallel  to  the  range 
influences  both  the  x  and  z  momentum  equations  through 
the  Coriolis  terms.  The  y  momentum  component  equation  is 
only  weakly  coupled  to  the  others  and  may  be  solved  in  a 
similar  way  to  the  temperature  equation. 

(2)  Additional  terms  from  the  Coriolis  force  enter 
the  momentum  equations  giving  rise  to  new  terms  in  the  vortic- 
ity  equation. 

4.1.1  Difference  Equations 

The  difference  equation  formulation  corresponding  to 
Eq .  (4.6)  and  supplementary  equations  can  be  chosen  to  parallel 
that  of  the  1IA.IEA  code.  The  equation  for  the  y-component  of 
vorticity  is  obtained  by  cross -differentiation  of  the  x  and 
z  components  of  the  momentum,  thereby  eliminating  the  pres¬ 
sure  from  the  equations  entirely.  The  resulting  equations 
are : 


2  n 


sinO 


2.Q  cosO  sin(J)T,—  - 

1  t)X 


fi_  3£ 
p  3x 
8 


> 


Dv 
3 1 


+ 


siu0(u-.u  ) 


20  cosO  SiiKf)  w  =  0 
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where  the  y - component  of  vorticity  n  is  defined  by 

_  3_u  9w 
n  9z  “  3x  ‘ 

Ihe  difference  approximation  for  the  above  equations  also  fol¬ 
lows  that  of  MAlbA.  The  time-dependent  equations  are  solved 
in  an  explicit  two-level  formulation.  The  advection  terms  are 
obtained  by  using  a  high  order  conservative  scheme  in  which 
the  two  directions  are  integrated  by  the  splitting  technique. 
Additional  terms  are  centered  in  space  from  quantities  avail¬ 
able  at  the  current  time  cycle.  All  of  these  terms  result 
fiom  the  Coriolis  force  acting  on  the  v-component  of  velocity. 
Ihe  additional  equation  for  the  v-velocity  is  similar  to  the 
temperature  equation  in  structure.  Changes  in  v  are  prin¬ 
cipally  from  advection;  the  forcing  terms  are  the  result  of 
the  Coriolis  force  acting  on  the  u-  and  v-components  of 
velocity . 
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5,  A  LINEAR  ANALYSIS  OF  NUMERICAL  APPROXIMATIONS 


The  application  of  the  equations  of  fluid  dynamics  to 
geophysical  phenomena  has  prompted  the  development  of  simpli¬ 
fied  equations  which  emphasize  the  salient  features  of  strongly 
sub-sonic  flow.  ihis  investigation  of  these  equations  treats 
only  the  limiting  case  of  small  amplitude  disturbances  about  an 
initial  steady  state. 

Since  sound  waves  play  a  very  minor  role  in  motions  of 
the  atmosphere  and  oceans,  it  is  convenient  to  neglect  terms 
which  contribute  principally  to  their  propagation.  The  effect 
of  buoyancy,  on  the  other  hand,  has  a  profound  influence  on 
the  flow  field;  highly  stable  stratification  gives  rise  to 
wave  motions  which  are  frequently  ubiquitous  features  of  at¬ 
mospheric  phenomena.  Less  stable  atmospheric  configurations 
may  result  in  the  dissipation  of  waves  through  instabilities 
which  aiise  from  the  amplification  of  wave  motions,  and  an  un¬ 
stable  atmosphere  gives  rise  to  spontaneous  motions  of  a  con¬ 
vective  nature  which  provide  a  mechanism  for  restoring  margi¬ 
nal  stability.  In  much  of  the  motion  of  the  atmosphere  and 
oceans  the  fluid  behaves  as  though  it  were  substantially  in¬ 
compressible.  Consequently,  methods  of  approximation  have 
been  designed  around  the  assumption  of  incompressibility  and 
make  greater  or  lesser  use  of  it. 

Ihe  fully  compressible  equations  are  treated  in  this 
report.  Comparisons  are  made  with  results  from  Dr.  M.  G. 
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Wurtele  at  the  University  ol.'  California  at  Los  Angeles  and  the 
differences  arc  discussed.  Following  thi  s ,  the  results  arc 
modified  to  the  incompressible  equations.  These  two  cases  are 
the  limiting  solutions  of  our  problem. 

Three  approximations  to  the  fluid  dynamic  equations 
are  then  solved.  These  include  (1)  the  strict  Boussinesq 
approximation,  (2)  the  "compressible"  approximation,  and 
(5)  the  anelastic.  equations.  The  Boussinesq  approximation 
is  based  on  the  assumption  that  the  fluid,  as  far  as  the  dy¬ 
namical  terms  are  concerned,  is  incompressible  and  has  a  uni¬ 
form  density.  The  "compressible"  approximation^1^  was  devel¬ 
oped  to  take  account  of  the  effects  of  the  stratification  of 
the  atmosphere  on  waves  which  propagate  through  vertical  dis¬ 
tances  comparable  to  the  atmospheric  scale  height.  The  third 
approximation  studied  is  the  anelastic  equations  introduced 
by  Ogura  and  Phillips.  These  equations  are  based  on  small 
departures  from  a  neutrally  stratified  atmosphere  in  which 
the  potential  temperature  has  a  constant  value. 

5.1  EULER'S  EQUATIONS 

Initially,  the  equations  for  the  motion  of  a  fluid 
obeying  the  perfect  gas  law  is  considered,  as  is  appropriate 
for  the  dry  atmosphere.  We  assume  that  the  gravitational 
force  is  constant,  that  viscous  or  turbulent  transfer  can  be 
neglected,  that  the  motion  :i.s  adiabatic,  and  that  the  Coriolis 
terms  have  a  negligible  effect.  In  addition,  the  motion  is 
assumed  to  be  confined  to  a  vertical  (x,z)  plane  in  which  x 
and  z  are  horizontal  and  vertical  Cartesian  coordinates. 

The  resulting  simplified  equations  arc: 
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Mass  Conservation  -- 


dp 

dt 


+ 


/Du 

N 


o 


(5.1a) 


llor  izontal  Momentum  -- 


du 

at 


+ 


1  DP 
p  Dx 


(5.1b) 


Vertical  Momentum  -  - 

dw  1  3JP  = 

dt  p  Dz  b  > 

Energy  Equati  on  for  Adiabatic  Motion  -  - 

dT  =  J_  dp 
dt  pC^  dt  ’ 


(5.1c) 


(5.  Id) 


Perfect  Gas  Equati on  of  btate  -  - 


p  =  pRT  ,  (5 .  le) 

where  the  substantive  derivative  is  defined  as 

d  _  D  9  D 

dtf  =  DT  U9x  Vdz  ' 

In  the  above  equations  there  are  five  dependent  vari¬ 
ables;  the  horizontal  component  of  velocity  u,  the  vertical 
component  of  velocity  w,  the  density  p,  the  pressure  p,  and 
the  temperature  T.  The  specific  heat  at  constant  pressure  C 
is  assumed  to  be  constant. 


70 


.vSR"  111 3  I 


The  unperturbed  state  oi‘  the  atmosphere  is  assumed 
to  he  horizontally  strati  lied  and  in  hydrostatic  equilibrium. 
IVc  consider  perturbations  to  this  system  in  which  the  dis- 
tui banco  amplitude  is  small  enough  that  only  terms  linear  in 
the  perturbed  quantities  need  be  retained.  The  resulting 
equations  admit  solutions  which  arc  periodic  in  the  horizontal 
variable  and  the  time.  Due  to  the  linearity  of  the  equations, 
these  wave  solutions  can  be  superimposed  to  describe  the  be¬ 
havior  of  an  arbitrary  disturbance. 

The  assumed  solution  takes  the  form 


p  =  p  +  pf  , 

(5.2a) 

u  <=  U  +  flf  , 

(5.2b) 

w  =  wf  , 

(5.2c) 

p  =  p  +  pf  , 

(5. 2d) 

T  =  T  +  Tf  , 

(5 . 2e) 

where  f  -  c  ^  ^  in  which  k  is  the  horizontal  wave 

number  of  the  disturbance  and  c  is  its  phase  speed.  The  re¬ 
maining  quantities  depend  only  on  the  vertical  coordinate  z, 
with  the  first  quantity  in  each  sum  representing  the  unper¬ 
turbed  solution. 

The  unperturbed  quantities  must  satisfy  the  equation 
of  state  and  the  hydrostatic  condition.  This  leads  to 

P  =  P]^  >  P'  =  -gp"  ,  (5.3) 

where  differentiation  with  respect  to  z  is  denoted  with  a 
prime . 
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in  the  equations  governing  the  perturbed  quantities, 
wc  retain  onJy  tcri!’.->  linear  in  the  perturbed  quantities;  t lie* 
unperturbed  terms  are  removed  by  subtracting  Hq.  (5.3). 
Factoring  out  the  common  factor  f  we  obtain 

ik(U  -  c)p  +  p'w  +  p(iku  +  w’)  =  0  ,  (5.4a) 


ik(U  -  c)u  +  U'w  +  ~  p  =  0  , 

P 


(5.4b) 


ik (U  -  C)w  +  =  0  , 

P  P 


(5.4c) 


ik(U  -  c) (pRT  -  <p)  +  w(pRT'  -  <p ' )  =  0  , 


JL 

P 


4"  =  0  • 


(5 . 4d) 
(5.4c) 


whe  re 


1  f 


From  Eq .  (5.3),  the  derivative  of  the  unperturbed  den¬ 
sity  and  temperature  can  be  related  as  follows: 

_H  +  JLl  b  _  ££  ; 

T  p  0 
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In  terms  o!  the  slr.it  i  !'i  cations  of  density  and  potent  ia 
pern  Mire  0 


s  - 


P ' 


and  S  = 


the  above  equation  becomes 


s  -  s  , 

s 


in  which  the  sound  speed  C  is  defined  as 


c2  =  -i.  _il_ 

s  .1  - «  - 

P 


and 


6  =  T(P/P0) 


-  K 


Eliminating  f  and  u  from  liq .  (5.4) 


ik(U  -  c)p  -  s  +  n~ — 7  pw  +  pw ' 


U  -  c 


U  -  c 


P  =  0  , 


gp  +  ik  (II  -  c)  pw  +  p'  =  0  , 


-ik(U  -  c)p  +  pSw  +  AiLCJJ.  -  c) 


c2 

s 


p  =  0 
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!;qs  . 

A  A  .  A 

p  ,  w  and  p 
gives 


(5.0)  are  three  equations  for  the  three  unknowns 
hi imina t i ng  the  density  from  these  equations 


w ' 


S  -  s 


U1 


u  - 


w  - 


ik 


U  -  c 


U3  =  0 


Sg 

(U  -  cp 


1-2 


W 


ikSB  .  ikB 1 
U  -  c  U  -  c 


(5.7) 


where 


B  =  p/p  and 


(5.7a) 


In  all  but  the  most  violent  storms  \i  1  .  Finally,  elimi¬ 
nating  B  from  bq ,  (5.7)  and  after  some  reduction,  the  equa¬ 
tion  for  w  is  obtained: 


'  (s  +  +  ['|Jk2  *  iir?arp'(s  +  r) 


+  u' 

i  s  +  ii 11 

.  U"  I 

U  -  c 

y  J 

U  -  cj 

bq.  (5.8)  is  the  same  as  that  of  IVurtcle  1S^  who  derives  the 
equation  from  the  formulation  in  terms  of  the  potential  tem¬ 
perature  and  the  Helmholtz  pressure  function 


It  is  of  interest  to  consider  the  limiting  form  of 
bq .  (5.8)  when  the  velocity  is  small  compared  with  the  sound 
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speed,  1  Ills  app roxi motion  is  basically  obtained  by  replacing 

M  by  1»  but  we  shall  retain  certain  terms  arising  from  the 
differentiation  of  p  .  Neglecting  terms  which  contain 


U  -  c 


we  obtain 


w 


"-[■ 


2.(s  -  S)(l)  -  cjU'lA, 


8 


r...2  ,  rs  .  y 

L  w~^y  it^t(2S  -  s) 


(5.9) 


UM  .  2 (s  -  S)(U’)21 . 


U  -  c 


w  =  0 


This  result  differs  from  that  of  Wurtele 


sw '  + 


[•“ 


+  fiS  +  sU'  _ 
(u  -  c)2  U  -  c  U 


Lv] 


w  -  0  ,  (5.9a) 


principally  in  the  coefficient  of  the  U’  term  which  contains 
an  additional  contribution  from  gp'/(U  -  c) 2  in  Eq .  (5.8). 
Ihe  remaining  additional  terms  in  Eq .  (5.9)  are  small  compared 
with  Eq .  (5.9a)  unless  extremely  large  wind  shear  is  present. 

1  he  Euler  equations  have  been  considered  above  in  con¬ 
siderable  detail  in  order  to  provide  a  basis  for  comparison 
with  the  analogous  equations  from  the  several  approximations. 

J he  subsequent  derivations  resulting  from  the  approximate 
equations  also  give  rise  to  equations  for  w  which  have  the 
same  form  as  Eq.  (5.9).  By  comparison  of  the  coefficients  of 
w  and  w'  it  is  possible  to  determine  how  the  vertical  wave¬ 
length  and  wave  amplitude  compare.  It  was  found  convenient 
to  transform  Eq.  (5.9)  in  these  evaluations  to  remove  the 
first  derivative  term.  The  transformation  to  accomplish  this 
is 
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linns  St  win  S*e  useful  in  formulate  the  hitler  equations  In 
len»»  of  the  dependent  variables  ©  nn«l  m  .  l!ol  lowing 
hurt  win's  tier  I  vat  Inn  the  equations  hoof  «j 


In 

* 


* 


« 


0 


* 


dll 

Tf 


4 


.  a« 
vSf* 


o  , 


CS. 11) 


dw 

3t 


*  *  *8 


» 


do 

Iff 


* 


0 


The  equation  of  stale  has  been  incorporated  into  the 
above  equations  so  that  the  density  does  not  appear.  If  the 
density  is  desired  as  an  auxiliary  quantity  it  can  be  obtained 
frora  the  equation  of  state.  Introducing  unperturbed  and  per¬ 
turbed  equations  for  0  and  r,  analogously  to  liq.  (5.2),  we 
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obtain  ll»i4  ettrfv*|«oii»ltii|*  to  the  hydroisuit  Sc  balance 

condition 


o 


The  I  ht  curbed  pcrioriMtSciii  equal  ion$  are 

F|H(U  *  cjl  *  f»wj  ♦  C;Cikn  ♦  w»)  *  o 

o' 

ik(ll  •  c)ti  #  II* w  ♦  Ikon  -  U  , 


ik(U  *  c)w  ♦  of*  -  *  q  f 

0 


(5.12) 


ik(IJ  -  c)0  ♦  o*w  »  0  , 

ill  which  Cs  depends  on  T  ,  while  in  Eq.  (5.1J)  it  depends 
on  T  .  Eliminating  u  and  0  from  Hq.  (5.12)  we  obtain 

[u-  +  cT  (ll  ’  C)J  w  -  (U  -  C)w»  =  -ikOnp  , 

(5.13) 

[rr-:\  '  k2(u  ■  c)]w  a  ■  • 

Ihe  elimination  of  tt  from  the  above  equations  leads  again 
to  Eq .  (5.9).  Ill  is  requires  differentiation  of  the  first  equa¬ 
tion  giving  rise  to  a  term  in  p'  . 

5.2  INCOMPRESSIBLE  EQUATIONS 

We  consider  the  system  composed  of  an  incompressible 
fluid  (d  /dt  =  0)  which  is  initially  stratified  in  density. 
The  equations  for  t lie  system  are 
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tltl  . 

1  Si*  . 

ai  * 

oaf-  <jzm  a  Sc 

p 

*lw 

at * 

1  §p 

ar-  is  #  e  «r 

P 

ii 

!t*“ 

iH  ♦ 

t)K 

T*  * 

‘  32 

The  corresponding  pcrturhnt ion  equations  arc 

ik(U  -  c)u  ♦  U‘w  +  p  «  0  , 

P 

ik(U  -  c)iv  +  £—  +  ^  =  0  , 

P  P 

ik(U  -  c)p  +  p'w  =  0  , 

iku  +  w'  =  0  . 


(5.14) 


(5.15) 


The  equations  corresponding  to  Lq.  (5.7)  resulting 
from  eliminating  u  and  p  are 


_U_'w_ 
U  c 


+ 


0 


» 


iks_(3_  ikg ' 

U  c  U  -  c 


(5.16) 
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lliesi*  w|ii«ti«ns*  differ  fniH  li«| •  (5.7)  in  two  respects  and  can 
lie  derived  froti  them  by  the  stibst  I  tut  ions  u  I  ami  H  •*  s  . 
The  equation  for  w  becomes 


h 


sw1 


[■' 


nr 


*K 
C  p 


all' 


r 


ii 


c 


m 


i"  v 

•’cJw 


0  .  (5.17) 


As  is  to  lie  expected,  all  of  the  sound  wave  terms  are  missing 
and  the  term  which  is  frequently  the  largest  is  changed.  The 
coefficient  of  g/(ll  -  c)s  ,  which  in  the  exact  equation  is 
the  static  stability  S  ,  is  replaced  by  s  . 


5 . 3  BOUSS 1  NliSQ  APPROX  IMAT J  ON 

In  the  Roussincsq  approximation,  the  variability  of 
the  density  in  the  atmosphere  is  neglected  in  all  the  terms 
except  the  buoyancy  term.  This  formulation  allows  tbc  equa¬ 
tions  to  be  transformed  to  introduce  the  vorticity  and  stream 
functions,  a  property  used  in  the  S3  numerical  solution.  The 
equations  are 


(5.18) 
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where  I**  ,  T*  sire  the  pcrlnrluit ions  of  the  pressure  and 
temperature  f run  the  corresponding  static  si tmo sphere  quanti¬ 
ty**  denoted  hy  the  suhscript  “aero."  The  density  sippenring 
in  the  momentum  equations  is  takes  to  he  constant.  In  the 
above  equations  there  has  been  an  approximation  made  in  the 
continuity  equation  through  the  neglect  of  the  change  of  den¬ 
sity  due  to  changes  in  both  the  static  and  dynamic  pressures. 
In  the  momentum  equations,  in  siddition  to  the  approximation 
of  a  constant  density,  the  density  perturbation  in  the  buoy¬ 
ancy  term  has  been  replaced  by  the  temperature  pevturbatiou; 
the  pressure  perturbation  has  been  neglected,  in  the  tempera¬ 
ture  equation,  the  effect  of  the  pressure  perturbation  on  the 
temperature  change  lias  also  been  neglected, 

lo  obtain  the  equation  for  small  disturbances  in  a 
stratified  atmosphere,  we  linearise  the  Boussincsq  equations 
liq.  (5.18)  about  the  unperturbed  state.  The  substitutions 
are 


u  =  U  +  uf  , 
w  =  wf  , 

r  =  ff  , 

p'  =  pf  . 

The  perturbation  equations  arc 

ik(U  -  c)u  +  U'w  +  —  p  =  0  ,  (5.19a) 

^  r\ 
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A  fit 

n.(u  -  c)w  •  ^  -  y  » «  , 

po  'o 

ik(U  -  c)T  ♦  (TQ  ♦  l*)w  «  0  , 

iku  +  v?*  »  0 

A  * 

lili  in  inn  ting  u  mul  V  wc  obtain 


-(I) 

-  c) w * 

+  U'w 

♦  iM  =  o  , 

po 

-[k!(U 

-  c)  - 

fiS  ' 

1 ♦  Ml .  o  , 

>  po 

U  -  c. 

where 


S  = 


T,  +  r 

T_ 


The  resulting  equation  for  w  is 


u"  i  T  k 2  +  fiS  U"  lu-  -  0 

w  +  L k  (u  -  cp  ir^rr  u 


(S.lDb) 


(5. 19c) 


(5. 1!UI) 


(5.20) 


(5.21) 


Eq.  (5.21)  can  also  be  derived  from  the  vorticity  and  stream 
function  form  of  the  equations. 

The  first  two  coefficients  of  w  in  Eq .  (5.21),  which 
in  most  cases  are  the  largest,  also  arc  present  in  the  corre¬ 
sponding  Euler  equation,  Eq.  (5.19).  The  third  coefficient  oi 
Eq.  (5.21)  also  occurs  in  both  equations.  Consequently,  the 
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verticil!  wnve  length  of  the  ltuiu*sincsq  equn I  i on  should  he  ap¬ 
prox  i  mate  ly  correct.  The  w*  term  of  Hq.  (5.9)  is  miss  in}' , 
however.  As  indicated  by  the  traits  format  ion  of  Hq.  (5.10), 
the  absence  of  this  term  corresponds  to  a  seating  of  the 
velocity  amplitude  by  a  factor  of  (p)  .  Thus,  the  ampli¬ 

tude  of  the  Houssincsq  approximation  velocity  is  expected  to 
fall  off  more  rapidly  with  altitude  by  the  above  factor,  which 
becomes  significant  when  waves  arc  propagated  through  an  at¬ 
mospheric  scale  height. 


By  comparison  with  Eqs.  (5.4)  and  (5.6)  it  is  clear 
that  the  w  coefficients  are  obtained  correctly  from  Eq. 
(5.19)  because  the  temperature  perturbation  terms  in  the 
energy  equation  and  the  equation  of  state  are  correctly  re¬ 
tained.  The  emitted  pressure  perturbation  terms  contribute 
principally  to  smaller  coefficients  of  w  .  The  absence  of  a 
term  in  w'  results  from  the  omission  of  the  density  strati 
fication  term  in  the  continuity  equation  and  from  the  approxi¬ 
mation  of  constant  density  in  the  momentum  equations. 


5.4  THE  "COMPRESSIBLE"  APPROXIMATION 

r  i  7i 

The  "compressible"  approximation1  *  1  was  developed 
to  take  account  of  the  effects  of  the  stratification  of  the 
atmosphere  on  waves  which  propagate  through  vertical  distances 
comparable  to  the  atmospheric  scale  height.  The  equations  are 
formulated  in  such  a  way  that  it  is  possible  to  reduce  them  to 
a  Poisson  equation  for  a  quantity  similar  to  the  stream  func¬ 
tion  and  a  prognostic  equation  for  a  vorticity- like  quantity. 

The  equations  are: 


3pu  3pw 
3  x  3  z 


(5.22a) 
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3pT* 

at 


Dpu  4 
at. 

;»pu2  apuw  .  ap  .. 

3x  uz  +  ax  0  * 

(5 .  22lt  j 

Dpw  Opuw  apw2  ap 

Tt“  ax  +  *  7z  R  '«p  * 

(S.22c) 

+  idul 

Dx 

*  sf  *  -(P0r  *  pT')k  . 

(5 . 2  2d  j 

P  -  r, 

p  RT  . 

(5 . 22e) 

The  quantity  T*  is  the  temperature  perturbation  from 
f  ,  the  initial  temperature  of  the  stratified  atmosphere.  As 
compared  with  the  Euler  equations  there  are  three  approxima¬ 
tions.  I  lie  continuity  equation  omits  the  term  containing 
Jt  *  In  ener8y  equation,  the  pressure  perturbation  term 
is  omitted  from  the  work  term.  And  in  the  equation  of  state, 
the  pressure  takes  its  static  value,  neglecting  the  pressure 
perturbation.  The  linearized  equations  corresponding  to  Eq, 

(5 . 22)  are 


ikpu  + 

ikUp  +  pw'  + 

p  '  w  =  0  , 

(5 . 23a) 

ik(U  -  c )  u 

-  ikUc-?  *  U 

■i  *  iM  ,  o  _ 

(5.23b) 

p 

P 

A 

A 

SR  + 

i  k  ( U  -  c )  w  + 

L-»  o  , 

(5.23c) 

p 

P 

ik(U  - 

c)T  +  T'w  + 

rw  =  o  , 

(5.23d) 

A 

A  M1 

£+1=0 

P  T 

• 

(5 . 23e) 
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l.l  Vitiating  T  am!  u  ,  we  obtain  an  equation  corresponding 
to  l!t| .  (5.6): 

iklVp  »  |s(U  -  c)  «  U'J iw  ♦  F(b  -  c)w'  *  i i.ft  ■  o  , 

P P  *  ik(U  -  e)pw  ♦  p*  «  0  ,  (5.24) 

•ik(U  -  c)p  +  pSw  ■  0 


By  comparison  with  liq.  (5.6)  we  find  that  liq.  (5.24)  lacks 
several  terms:  in  the  first  equation  the  coefficient  oT  p 
contains  U2  rather  than  (U  -  c)?  ,  and  in  the  third  equa¬ 
tion  the  coefficient  of  p  is  zero.  Corresponding  to  liq. 
(5.7)  we  eliminate  p  arriving  ac 


where 


The  most  serious  difference  between  liq.  (5.25)  is  in  the 
second  equation  where  t lie  coefficient  of  3  contains  s  in¬ 
stead  of  S  .  Eliminating  3  we  obtain  the  equation  for  the 
vertical  velocity  amplitude 
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17  H  1 

r«- 

2sj 

£•  *  Jh"  ♦  |.j 

% 

p.  * 

IS  -  cjp 

II" 

x\at  .«aa:'^as.-  : 

ir  -  c 

ip  i 

\r=-t\vr 

U  |ei 

zxl.  ar:,aa  I  IJi 

-  cj‘* 

(2  -  ff * rcj  * 

■  (irH) 

sS  (.**.26) 

(r^r) 

[ 

S'  - 

s'] 

•*> 

w  *  l)  . 

The  leading  terns  in  the  coefficient  of  w  it)  l£«| • 
(5.26)  are  in  agreement  with  liq.  (5.0),  thereby  assuring  that 
the  vertical  wavelength  will  he  suhslnnt  ial  ly  correct.  A 
number  of  additional  terms  also  appear  in  the  coefficient  of 
w  but  in  most  cases  these  will  he  small.  The  major  discrep¬ 
ancy  of  Kq .  (5.26)  as  compared  with  Kq,  (5.11)  is  the  disagree* 

meat  in  the  leading  coefficients  of  w"  ,  in  which  (when  eE0 

and  the  smaller  IP  terra  is  neglected)  the  correct  coeffi¬ 
cient  -s  is  replaced  by  -2s  ♦  S  .  In  order  to  compare  the 
effect  of  this  difference  in  the  w '  coefficient  It  is  useful 
to  perform  the  transformation  of  w  in  Kq.  (5.26)  to  a 

quantity  x  obeying  an  equation  of  the  form  x"  ♦  v2x  *  0  . 

The  required  transformation  (when  c»0)  is 


w 


TW* '  =nWx 


(5.27) 


Since  the  x- equation  is  substantially  the  same  as  that  for  the 
transformed  liulcr  equation ,  as  given  in  the  equation  following 
liq.  (5.10),  the  two  formulations  may  he  compared  entirely  on 
the  basis  of  the  traits  forma  lions  of  liq.  (5.27)  and  liq.  (5.10). 
The  "compressible"  approximation  velocity  amplitude  contains 
an  additional  factor  of  relative  to  the  litilcr 

equation.  Consequently,  in  the  atmosphere  the  "compressible" 
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lift mK*  a* 8  I  S  Increase  with  8 mc* r*?«s§ 8 «g  silt  limb*  by  a  factor' 
o  f  ( p/  |»t| )  *  **  *  ^ '  re  I  a  1 1  vi*  1 1*  l  lie  lin  I  e  r  resu  1 1 ,  which  in  isarei 
Increases  wills  altitude  re  I  a  1 1 ve  to  (ho  ll©usslne$q  approx I 
nation  h>°  a  factor  of 


o 


) 


i  n 


s.s  tiiii  AXKUKrit:  Amttsimiio# 

Tin?  "aneJastic"  approx  I  nail  on  to  the  fluid  dynastic 
equal  lens  was  introduced  In  Ogura  and  Phillips.  This  ap¬ 
proximation  I s?  formulated  in  terms  ©f  the  potential  tempera® 
lure  and  Helmholtz  function  which  appear  in  the  form  of  the 
holer  equal  ions*  given  in  Iq.  (S.li).  The  ancSostic  approx  I  na¬ 
tion  is*  baaed  on  small  departures  from  a  neutrally  stratified 
atmosphere  in  which  the  potential  temperature  has  the  constant 
value  ©  .  The  pressure  gradient  terms  of  the  momentum  equa¬ 
tions  are  approximated  by  neglecting  higher  order  terms*  in 
the  perturbed  quantities  ©*  and  s*  defined  as 


©  *  ©  ♦  o' 


SI  a  HI  *  US* 


The  moment  sea  equations  then  become 


CS.iti) 


where  is  defined  from 


dti 

JIT  * 


,  m  . 

[£  i.  aa.  it 


(I 


dii 


3X 


%  *bs 

-  ar 


?5x~ 


dw 

m 


9  V**6 


dw 

iff 


-Sb1 


dw 

in 


Bll 


'1  « 


♦ 


g  a  0 
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tilt*  rgy  equation  Is*  solved  exactly,  whi  U*  tin*  tsass  const r* 
vat  ion  equation  ha  approx i ip  Sod : 


d§ 

iff 


S  (I 


» 


0 


(5.29) 


According  to  Ogttra  and  l*lti  1 1  Spa ,  the  density  ha  that  assoei* 
tiled  it  it'h  the  neutrally  sti  till  fled  atmosphere.  This*  ssys  tea 
of  equations  etisi  also  he  applied  to  tin  atmo*  pliere  which  is 
not  initially  neutrally  stratified.  lit  this  case  there  will 
he  a  static  so  hi  lion  in  which  C  and  iff*  have  non* zero 
values,  w  and  I  ,  which  defend  on  z  alone  satisfying  the 
relation 


(5.50) 


fte  now  consider  the  sisal  I  amplitude  behavior  of  liqs. 
(3.2$)  and  (5.29)  by  linearizing  the  solution  about  a  steady 
state  in  which  the  atmosphere  is  arbitrarily  stratified  in 
Itarallei  flow.  Ihe  dependent  variables!  are  written 

it  *  IS  *  tif  , 


w  «  wf  , 

c  ■  *  t  Jr  , 

E  =  i  *  «f  , 
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llss i #»|*  l»t| .  (5. SO)  for  tin*  glntig  $30liitS&n  lit4  @hUiin 


ik(ll  »  c)ii  •  [V$  •»  ikCJffi  *  a  , 


Jk(U  -  c)fi  ■*  Ot»  •  *  0  , 


i k (II  *  a) 6  ♦  Pw  5  0  , 


(S. SI) 


Iku  *  «•  4  S_  „■  ,  0  , 


'  0 
o 


ill  which  tli  fferen  tint  ion  with  respect  to  z  is  denoted  by  :i 
prime*.  1:1  imiiKitinj;  u  ;unl  6  wo  obtain 


•«’  *  [-0  *  U~c]''  *  iFc 


i  l;0«t 


p’So  t  ?l  -  .  il;Oa  * 


*  0  , 


Ko  .A* 
;cg”^“cf  ’  1  r 


♦  r^c  - 0  • 


(5.32) 


where  the  parameters  appearing;  are  defined  by 


P 1 
Ko 


o  G 


o 


(5 . 33) 


The  equation  for  w  is 


w"  -  s^w  ’  »  £-h: 


RS 


o 


S«U’ 


ir 


(II  -  c)*'  U  -  c  II  -  c 


•  si]«  ■ 


(5.34) 
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Ccmparhig  Hij.  (5.54)  with  the  correct  )U| .  (5.9)  we  find  that 
Cn«t«l«r  f run  small  terns)  the  coom  cirnts  or  w  mu!  w*  arc 
suhstaitt  I  a)  ly  in  the  hsihc  form.  The  correspondi  ng  terms 
differ,  liowovei;  SQ  is  defined  in  terms  of  0  instead  of 
'  ♦  y  and  s0  is  formed  from  pQ  ,  the  neutrally  stable  den¬ 
sity,  which  differs  from  p  the  density  of  the  actual  atmos¬ 
phere. 


Fortunately,  the  anelastic  static  stability  S  does 
not  differ  by  a  1  a rge  amount  from  the  actual  value  of  S  . 
Moreover,  the  atmospheric  density  gradient  term  s  can  be 
chosen  to  have  the  correct  value  s  rather  than  that  associ¬ 
ated  with  the  neutrally  stratified  atmosphere  as  specified 
hy  Ogui a  and  Phillips.  Jhc  density  enters  only  in  the 

mass  conservation  equation  where  it  occurs  as  a  parameter; 
consequently  there  is  no  difficulty  in  assigning  the  actual 
static  value  to  it. 


In  contrast  to  the  "compressible”  approximation,  the 
anelastic  approximation  thus  can  be  adjusted  to  give  the  cor¬ 
rect  value  of  the  w*  coefficient.  Inspection  of  the  deri¬ 
vation  through  liqs .  (5.31)  and  (5.32)  shows  that  this  takes 
place  through  the  compensation  of  two  approximations  of  the 
anelastic  equations.  First,  a  term  proportional  to  S  is 
omitted  when,  on  differentiation  of  the  horizontal  momentum 
equation,  the  coefficient  of  ir  is  approximated  by  the  con¬ 
stant  G  .  The  missing  term  just  compensates  for  the  differ¬ 
ence  between  sq  ,  resulting  from  the  appearance  of  pQ  in 
the  continuity  equation,  and  the  correct  term,  g/C*  ,  arisin; 
from  the  density  perturbation,  which  according  to  Eq.  (5.5), 
incorporates  an  S  term. 
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5.6  CONCLUSIONS 

Several  approximations  to  the  inviscid  fluid  dynamic 
equations  for  treating  buoyant  low- speed  flow  have  been  ex¬ 
amined.  In  each  case  the  equation  governing  the  small  ampli¬ 
tude  vertical  velocity  perturbation  lias  been  derived  and 
compared  with  that  for  the  Euler  equation.  When  waves  are 
propagated  through  vertical  distances  comparable  to  the  den¬ 
sity  scale  height  of  the  atmosphere  the  familiar  Boussinesq 
approximation  gives  rise  to  amplitude  errors.  The  "compres¬ 
sible"  approximation  introduces  a  correcting  term,  but  the 
coefficient  is  too  large,  permitting  waves  to  propagate  up¬ 
ward  too  strongly.  A  readily  accomplished  modification  of 
the  "anelastic"  approximation  gives  rise  to  the  correct  coef¬ 
ficient,  while  retaining  a  substantially  correct  vertical 
wavelength  as  well. 

The  modified  approximation  appears  to  be  an  attrac¬ 
tive  candidate  for  treatment  of  problems  of  deep  convection 
or . propagation  of  gravity  waves  through  appreciable  fractions 
of  the  atmospheric  scale  height.  The  behavior  of  the  anelas¬ 
tic  approximation  under  conditions  in  which  large  amplitude 
disturbances  are  produced  is  not  currently  known. 


96 


3SK-1034 


6.  RADIATION  III  THE  EARTH'S  ATMOSPHERE 


The  radiative  transfer  problem  in  the  Earth's  atmos¬ 
phere  reduces  to  the  solution  of  the  seemingly  simple  equa¬ 
tion 


XT  =  Jv  '  Kv  !v 


which  states  that  the  radiant  intensity,  in  traversing  the  ele¬ 


ment  of  length  ds  ,  will  be  augmented  by  sources  in  the  amount 


Jv  ds 


and  diminished  by  extinction  in  the  amount 


Kv'v  ds 


In  general,  1^  ,  the  radiant  intensity,  and  J  ,  the  source 
function,  depend  on  both  a  spatial  coordinate  r  and  an  angu¬ 
lar  coordinate  (direction)  $  at  the  point  r  ,  as  well  as 
upon  the  frequency  v  .  The  time  dependence  of  these  quanti¬ 


ties  is  ignored  because  the  radiative  state  of  the  atmosphere 
is  established,  for  all  practical  purposes,  instantaneously. 

Kv  is  the  extinction  coefficient,  which  describes  the  rela¬ 
tive  depletion  in  the  intensity  of  the  beam,  dl v/ I v  ,  upon 
traversing  the  element  of  distance  ds  .  is  in  general 

the  sum  of  an  absorption  part  and  a  scattering  part.  Jv  de¬ 
scribes  the  additions  made  to  the  beam  intensity  along  ds  by 
thermal  (I'lanckinn)  emission  and  by  scattering. 
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6.1  INTRODUCTION 

A  substantial  revision  in  our  fundamental  numerical 
technique  for  solving  the  radiation  transport  equation  has 
taken  place  since  the  previous  semiannual  report  of  this  con¬ 
tract.  ^  This  revision  has  led  to  some  new  numerical  prob¬ 
lems  and  has  delayed  the  intended  completion  date  of  the 
radiation  code,  but  we  consider  the  new  method  sufficiently 
superior  in  both  accuracy  and  computational  economy  to  more 
than  justify  the  delay. 

The  first  part  of  the  new  method  involves  fitting 
the  transmission  function  by  a  sum  of  exponentials  and  is 
discussed  in  Section  6.2.  The  second  part  involves  the  solu¬ 
tion  of  the  monochromatic  equation  of  vadintivc  transfer, 
which  is  the  subject  of  Section  6.3.  In  Section  6.4  the  treat¬ 
ment  of  Mie  scattering,  in  particular  the  treatment  of  the 
forward  scattering  peak,  is  described.  Finally,  a  block 
diagram  of  the  code  is  given  in  Section  6.5,  along  with  brief 
descriptions  of  the  various  blocks. 

Minor  modifications  have  been  made  in  the  equations 
presented  in  the  previous  report.  First,  the  origin  of 
coordinates  is  now  taken  at  the  top  of  the  atmosphere;  this 
facilitates  the  method  of  Section  6.3.  Secondly,  a  new  dif¬ 
fuse  intensity  i^  has  been  defined  as 


i 


v 


y  solar 


-1  <  U  <  0 
0  <  p  <  1 


(6.1) 


where 
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T  solar  .  *v  c‘T(i!j/MsuM  4(  , 

v  2n  Hsuir 


(6.2) 


T(2) 


/V* 


)  dz'  . 


With  this  transformation,  the  monochromatic  equation  of  trans 
for  derived  in  the  previous  report  assumes  the  form 

di  0 

v7z~  *  Kviv =  aiBv  *  2~  J i  V.,y,y)ivc*,y)di.' 


(6.3) 


♦  Q^U.m) 


where 


8vSv  r-  *t(z)/m 

Qv,C  =  ^)  *  -ITT  Fv(z»M*y«»n)  e 


sun" 


The  integral  form  of  iiq.  (6.3)  is 


iv(ZiU)  ■ 

jf) .  i  r  (ju) .  j(2) . 

V  \l  J,  1  V  V  V  1 

*zo 

(6.4) 

where 

j'0)  - 

MV°  exp[-  ijx(2)  -  t(20){] 

(6.5a) 

j(i)  . 

=  Bv|TU")]  OX|>[-  A|t(2)  -  1(2")}  j 

(6.5b) 
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(2)  . 


v 


M2">  r1  _ 


jT  rv(2",  p,p')  i  v  (  z** »  M  *  )  dv»  * 
•  exp£-  £Jt(z)  *  T(z")jJ 


(6.5c) 


(3) 


sv3v(zu)  ^ 

•^4?—  Fv<2"’^sun> 


cxp['^T(z,,)  ‘^T(2)  ‘T(2,,)!] 


(6 . 5d) 


The  zq  of  Eq.  (6.4)  need  not  lie  on  a  boundary  of  the  sys¬ 
tem;  it  need  only  be  less  than  z  if  p  >  0  and  greater 
than  z  if  p  <  0  . 


6.2  EXPONENTIAL -SUM  FIT  OP  THE  TRANSMISSION  FUNCTION  AND 

REDUCTION  TO  MONOCHROMATIC  PROBLEMS 

Let  us  consider  the  frequency- averaging  of  Eq.  (6.4) 
over  a  frequency  interval  Av  which  may  contain  many  (typi¬ 
cally  20-100)  absorption  lines  of  1120,  C02,  etc.,  but  which 
is  narrow  enough  that  B  ,  F  ,  3X. ,  a'/cont',  and  S  do  not 
change  appreciably  across  Av  .  (a^  J  is  the  smoothly 

varying  part  of  the  absorption  coefficient,  due  to  Mic  absorp¬ 
tion,  the  H20  continuum,  the  N2  continuum,  etc.)  Let  the 
subscript  Av  denote  the  frequency  average  of  a  quantity. 

For  example, 


Av 


l  *v 


dv 


(6.6) 


In  frequency-averaging  Eq.  (6.4),  the  frequency  integration  is 
commuted  with  any  factors  containing  only  the  smoothly  varying 
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quantities  l’v  ,  etc.  by  evaluating  these  quantities  at  the 

mid-point  v  of  Av  .  The  two  exceptions  to  this  rule  are 

the  Planck  function  and  solar  flux  S  ,  for  which  B 

v  Av 

and  SAv  arc  used  instead  of  B-  and  S-  .  The  reasons  for 
these  exceptions  are  (1)  if  the  Earth  radiated  as  a  blackbody 
the  total  radiation  output  to  space  would  be  given  correctly 
by  the  scheme  using  BAv  and  not  by  the  scheme  using  B-  , 
and  (2)  the  total  incident  solar  flux  is  given  correctly  by 
the  scheme  using  SAv  and  not  by  the  scheme  using  S—  .  B^v 
is  computed  correctly  to  about  5  decimal  places  by  a  combina¬ 
tion  table  look-up  and  asymptotic  formula  subroutine  called 
PLKAVG  which  has  been  developed  especially  for  this  code. 

SAv  is  computed  from  the  data  of  Thckackara. ^ 1 


According  to  Kirchhoff's  Law, I19)  which  says  that  media 
which  absorb  strongly  at  a  certain  frequency  emit  strongly  at 
that  frequency,  the  intensity  iy  will  tend  to  follow  the 
rapid  variations  in  the  absorption  coefficient  a»0*ne)  , 

If  sufficiently  many  (>  1-2)  mean  free  paths  of  line  absorp¬ 
tion  intervene  between  the  point  z  and  the  surface,  the 
variation  of  i^(z,p)  with  v  will  follow  the  variation  of 
a^(linc)  v  fairly  faithfully;  that  is,  lines  in 


a^(linc)  Wiu  ,)r(Hiucc  bright  lines  in  the  spectrum  of  inten¬ 
sities.  Nearer  the  surface,  lines  in  m.,y  produce 

dark  lines  in  the  upward  directed  intensities  rather  than 
bright  lines,  but  the  statement  that  iy  varies  as  rapidly  as 
still  applies.  Because  of  this  rapid  variation,  i 


.(line) 


cannot  be  commuted  with  the  frequency  integration  in  the  terms 
and  .  In  the  previous  report,  it  was  proposed  to 

make  the  approximation  that  this  commutation  is  valid,  in 
order  to  facilitate  the  proposed  method  of  solution.  The 
argument  was  made  that  the  errors  so  produced  would  have 


little  impact  on  total  fluxes  because: 
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(1)  in  the  case  of  ,  only  errors  near 

the  boundary  zQ  would  be  produced  and 
these  would  damp  out  rapidly  as  we  pro¬ 
ceeded  away  from  the  boundary ; 

(2)  in  the  ease  of  the  scattering  source 

(2) 

Jv  ,  only  a  negligible  fraction  of  the 
frequency  groups  which  experience  line 
absorption  will  have  comparable  amounts 
of  scattering  (J^2^  *  J^1))  ,  so  that 
even  a  large  error  in  J^2^  will  not 
be  felt. 

These  hypotheses  arc  qualitatively  correct,  but  a  quantita¬ 
tive  verification  would  demand  a  prohibitively  expensive  line- 
by-line  calculation  through  the  entire  IR  spectrum,  which  we 
were  not  prepared  to  do. 

Therefore,  it  was  fortunate  that  a  literature  search 
following  the  publication  of  the  previous  report  uncovered 
a  method  for  performing  frequency- averaging  which  eliminates 
the  need  for  the  above  hypotheses  and  instead  requires  cer¬ 
tain  others  which  arc  both  more  intuitively  appealing  and 
more  easily  verifiable.  The  technique  is„  simply  stated, 
the  fitting  of  the  transmission  function  with  an  M-term  expo¬ 
nential  sum  followed  by  the  solution  of  M  monochromatic  prob¬ 
lems,  the  sum  of  whose  solutions  gives  the  frequency-averaged 
intensity  iAv  .  We  give  the  details  below. 

After  performing  the  commutation  of  the  frequency 
integration  with  the  slowly-varying  functions  of  v  in 
through  ,  the  following  kinds  of  integrals  remain; 


*Ta<zl»z2 


)/0 


dv 


(0.7) 
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T„(2i  »  2  2  )  /  lJ 

312  dv 


(6.8) 


*a(Z| »2^) 


dv 


(6.9) 


where  the  line  absorption  optical  depth  t  is  defined  as 

a 


Ta(: 


1  ’  2 


. (line) 

v 


(z)dz 


TAv  (***•  *-s  callc^  a  transmission  function.  A  compu¬ 

tationally  fast  and  reliable  computation  of  TAv  has  been 
made  available  to  us  in  the  form  of  the  subroutine  I.OTKAN 
developed  at  AFCRL.^2^  LOT  KAN  was  described  in  the  semi¬ 
annual  report,  and  prior  to  that  date  had  been  put  on  our 
computer  but  not  adapted  to  the  particular  needs  of  the  radia¬ 
tion  code.  Since  then  a  substantial  amount  of  coding  having 
to  do  with  Rayleigh  scattering  and  Mic  scattering  and  with 
atmospheric  curvature  (which  is  being  neglected)  has  been  re¬ 
moved  from  LOTRAN  and  the  numerous  associated  data  cards  have 
been  converted  to  FORTRAN  DATA  statements  to  avoid  having  to 
read  data  cards  every  time  the  code  is  executed.  Further 
modifications  of  LOTRAN  will  be  necessary  in  connection  with 
the  exponential  fitting  scheme  discussed  below.  Consequently, 
it  can  be  assumed  that  accurate  values  of  for  arbitrary 

zl»z2  arc  available. 
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6.2.1  The  Scaling  Approximation 

LOTRAN  assumes  that  optical  depths  arc  calculated  in 
the  scaling  approximation,  according  to  which  a^(l*ne)  js 
separable  as  follows; 


a. (line) 


(■ 


-hv/k'A 

c  1  ax(v)  a 2 (p  »T) 


(from  here  on  the  stimulated  emission  correction  factor 
1  -  exp  (-hv/kT)  will  be  evaluated  at  v  =  v  since  it  is  a 
slowly  varying  function  of  v  ).  This  approximation  gives  less 
accuracv  in  optical  depth  computations  than  the  familiar 
Curtis -Godson  approximation,1  *  1  but  Zdunkowski1  J  has  re¬ 

cently  shown,  using  exact  line-by-line  transmission  function 
computations,  that  the  scaling  approximation  is  accurate  to 
at  least  two  decimal  places  for  a  wide  range  of  vertical 
atmospheric  paths.  Also,  it  proves  much  easier  in  practice 
to  carry  through  the  analysis  of  Sections  6.2  and  6.3  with  the 
scaling  approximation  than  with  the  Curtis-Godson  approxima¬ 
tion.  In  the  latter,  the  coefficients  and  exponents  k^ 

in  an  exponential  sum  fit  to  T^  become  functions  of  the 
Curtis-Godson  mean  pressure  p  ,  which  in  turn  is  a  function 
of  the  levels  Zj  and  z^  between  which  T^v  is  desired; 
in  the  former,  the  a^  and  k^  depend  only  on  the  frequency 
group  Av  ,  simplifying  the  analysis  and  subsequent  programming 
considerably. 


The  optical  depth  in  the  scaling  approximation  is, 
following  Goody,  ^  ^ 
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JfZ7  -hv/kT 

(1  -  e  )  a2(p,’J')dz 

Z1 


=  a(li«e) 


(P0»V 


(1 


-hv/kT 

e  ) 


a2 


(P,T) 

TiTTi^! 


■dz 


(line) 

= 


(Po  -V 


uv  (Z1'Z2: 


where 


(6.10) 


u?(zi,z2) 


-hv/kT 

e  ) 


a? (P  >T) 


(6.11) 


u-  is  called  the  "equivalent  absorber  amount."  p  and  T 
V  o  o 

arc  a  reference  pressure  and  temperature  (e.g.,  STP) .  Most 
authors,  including  McClatchey ,  make  the  further  assump¬ 
tion  that  a2  is  a  power-law  function  of  pressure  and  temper¬ 
ature  ; 


ct2(p,T)  =  ApY  T6 

McClatchey,  et.al.  pick  values  of  y  and  6  which  lead  to 
an  optimal  fit  of  the  exact  transmission  function.  Different 
values  of  y  and  6  are  used  for  water  vapor,  for  ozone,  for 
the  uniformly  mixed  absorbers  (C07,  N20,  CO,  CH^,  etc.),  for 
the  nitrogen  collision-induced  continuum  around  4p,  and  for 
the  water  vapor  continuum  between  8y  and  13p.  (The  last  two 
properly  fall  into  the  province  of  Ci^(cont)  f  but  are  men¬ 
tioned  here  for  completeness.) 
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McClatchey,  ct.nl.  do  not  I uc halo  the  stimulated  emis¬ 
sion  correction  factor  in  their  definition  of  u-  ,  so  that 
their  u-  docs  not  in  fact  depend  on  frequency  at  all.  The 
reason  for  this  may  he  that  1.0TKAX  was  not  conceived  as  a 
tool  for  radiative  transfer  calculations  per  se ,  hut  rather 
as  a  means  of  obtaining  transmission  functions  for  laser  ap¬ 
plications,  IK  photography,  etc. 

6.2.2  lixponcn t ini -Sum  Pitting  as  a  hcbesr.ne  Quadrature  Rule 

The  transmission  function  (6.7)  can  he  written  in 
the  scaling  approximation  us 


a 


-kvti-/u 


dv 


where 


kv  5 


(1 inc) 


(IVV 


Presume  now  that,  at  a  sequence  of  values  u^(iEl,  ...,  N'u) 
of  u-  ,  we  fit  l'Av  by  a  sum  of  exponentials: 


•k.u/p 

e 


(6.12) 


The  fit  could  he,  say,  in  the  least  squares  sense.  M  is  en¬ 
visioned  as  being  fairly  small,  on  the  order  of  5-10.  Such 

ric  *>(-,  ?7  ■>»!  . 

fits  have  been  discussed  in  the  literature1  “  *“  * **  1  for  at 

least  two  decades,  although  the  application  is  usually  to  a 
whole  absorption  band  rather  than  a  fairly  narrow  group  such  as 
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Ls  being  considered  heie.  Methods  for  computing  the  a j  and 
k|  ,  knowing  TAv  ,  will  be  discussed  in  Section  6.2.4. 

Kewr i  t i ng  Hq .  (6.12)  us 


1 

Sv 


-k  u/u 
e 


-k.ii/ii 

e 


(6.13) 


it  becomes  clear  that  the  sun  on  the  right-hand  side  of  I5q . 
(6.13)  may  be  interpreted  as  a  quadrature  rule.  Since  M  is 
to  be  roughly  5-10,  the  quadrature  cannot  lie  of  the  familiar 
Gaussian  or  New ton -Cotes  type,  because  such  quadrature  rules 
would  require  literally  thousands  of  terns  to  accurately 
integrate  the  line  structure  iu  e'*V,/u  (this  is  why  line- 
by-line  integration  of  transmission  functions  is  so  expen¬ 
sive).  The  rule  is  rather  of  the  l.ehesgne  type,  in  which  the 
ordinate  in  the  graph  of  the  integrand  is  partitioned  rather 
than  tile  abscissa  (for  a  fuller  discussion  of  l.ebesgue  quadra¬ 
ture,  refer  to  Appendix  A).  The  Uj  are  interpreted  as  the 
fraction  of  the  interval  Av  over  which  kv  is  near  k.  . 
Prom  this  interpretation,  it  is  clear  that  we  shall  wish  to 
restrict  the  tij  to  he  positive  iu  any  numerical  fitting 
scheme . 

As  pointed  out  in  Appendix  A,  the  weights  (here  the 
aj)  in  a  Lebcsguc  quadrature  rule  depend  on  the  function  be¬ 
ing  integrated.  However,  large  classes  of  functions  may  be 
quite  accurately  integrated  by  the  same  set  of  a-  .  based 
on  our  earlier  arguments  that  iv  and  kv  tend  to  vary 
together,  we  make  the  conjecture  that  iv  e‘*iv,,^u  belongs  to 
the  same  class  as  c  kvu/u  |„  {)K,  ni|0V£,  sense  and,  there¬ 
fore,  may  he  integrated  by  the  same  quadrature  rule: 
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1 

Av 


-k  u/m 
c 


dv  = 


M 


ai 


•k. u/m 
c 


(6.14) 


We  further  conjecture  that  the  intensity  iv  alone,  and  the 
integrand  in  Eq.  (6.‘)),  may  be  similarly  quadratured: 


i 


Av 


i . 


1 


(6.15) 


fc/. 


exp 


V°'ZP  'ralzl-z2) 


sun 


dv 


•fcl 


exp 


-  k 


(^(O.zp  u-(Sj,2p 


V  M 


sun 


dv  (6.16) 


V'  [,  (Uv(0’ZP  uv(zl'z2} 

ti  l  <  Msun  " 


The  intensities  i^  are  as  yet  unspecified.  However,  since 
each  i^  is  associated  with  an  "absorption  coefficient"  k^  , 
it  is  a  reasonable  conjecture  that  i^  is  the  solution  of  a 
monochromatic  transfer  problem  with  absorption  coefficient 
.  This  conjecture  is  verified  in  Section  6.2.3. 

6.2.3  Splitting  into  Monochromatic  Problems 

Using  the  quadrature  formulas  of  Eqs.  (6.13)  through 
(6.16),  the  frequency  integrals  of  .1^  through  (see 

Eqs.  (6.4)  and  (6.5))  may  be  approximated  as 
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jO»  •k 
Av 


=  exp  -  ~  J  k^C)  (z')dz 


ElM  "Mrr(z  0»z)/M 

ai  h  e 

i  =  1 


‘  vBiv(T(z"))  exp[-  \£  8-(z')dz'j 


d 

3z" 


exp 


•  l  £, 


M 


£  ai 0 
i  =  l 


-k.u-(z’\z)/p 


! 


(2)  *  *v 


3-(zM)  r+i 


'Av 


/; 


dp*  P-(z",p,p') 


•  exp 


i  r  ki' 

P  Jz  „  V 


^-c)  (z')dz 


M 


i  =  1 


i.(z’',p*) 


-k.u-(z",z)/p 


(6.17) 


(6.18) 


(6.19) 
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j(3)  ‘k  Sav^"(2  j  v  } 

J 1  5-  *vu  ,li’MsunJ 


'Av 


Tit 


•  c*p|-  IT1-  fZ  kvc),1::,  -  5  f  „  kvc),lz 

|  ,  [V’-2")  VZ">Z>|) 

•  L,  \  cx"j-  ki|“~  ~  tj 


(6.20) 


where  the  continuum  part  of  the  opacity  is 

k(c)  ,(cont)  , 

V  V  V 


Now  suppose  that  the  approximations  of  liqs.  (6.17)  through 
(6.20)  arc  put  into  the  frequency-average  of  the  transfer 
equation  (6.4): 


(2»M) 


(z»m)  *  J 


( 0 ) 

Av 


♦ 


♦ 


♦ 


The  summations,  and  the  coefficients  a ^  (which  are  altitude- 
independent  in  the  scaling  approximation),  may  he  commuted 
with  the  altitude  integration  on  the  right-hand  side.  Then 
if  we  regard  the  i^  (i»l,  ...»  M)  as  unknowns,  the  simplest 
(hut  by  no  means  the  only)  way  to  determine  them  from  the 
last  equation  is  to  equate  the  sums  on  each  side  term-hy- term: 
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ij (2  •  m)  =  ijU^iO  exp  |-  TiC2o,2)/u|  *  j  d2"  BAv(T(z" 
•  cxp|-  T(S)(2",2)/pj  —TT  CX|>  |  -  T[a)(z",z)/w| 


i  rz  mz")  rl  ~ 

'  \  L3  ^  -  v-  L 

•  ii(2",JI,)dM'  0X|>|-  Ti(2'\2)/pj 

♦  i  jf2  d2-  sj»Wz2  Vz,,^^su„) 


(6.21) 


(0  »***) 
^sun 


tiCs-.a) 


M 


(i»l,  . ..,  M) 


in  which  we*  have  defined  the  scattering  optical  depth 
absorption  optical  depth  Tp'*f  and  total  optical  depth  Ti 
as 


T(s) 


i)(z1.z2)  *  / 


0-(z)dz 


(cont) 


(z)dr  ♦  ki  u-(z1,22)  (6.22) 


Ti(21,z2)  a  '^(si»22^  +  Ti*,^:il,z2^ 


% 
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If  wc  further  define 


a. (2")  5 


5  •  I?1 


0>_(cont)(l„j  .  k.  u_(2..,2) 


(6.23) 


cil(coiit)(2„)  .  ki|l  .  c 
02(p(S"),l(2’')) 


-hv/kT(z”) 


then  the  second  (thermal  emission)  term  on  the  right-hand  side 
of  liq.  (6.21)  can  be  rewritten  as 


^(z") 


bAv(T(*”))  cxp["  Tit2"*2)/*1] 


which  completes  the  formal  identification  of  liq.  (6.21)  with 
the  monochromatic  transfer  equation  (6.4).  The  changes  from 
liq,  (6.4)  arc  simply  that  Fy  ,  a^cont^,  and  0V  arc  evalu¬ 
ated  at  v  ,  Bv  and  Sy  are  averaged  over  Av  ,  be¬ 

comes 


i( 


k,  11  - 


-hv/kT\  «,(|>,T) 

°  F&F7 


and  iv  becomes  i^  .  The  solutions  of  the  M  monochromatic 
liqs.  (6.21)  arc  combined  according  to  liq.  (6.15)  to  yield  the 
frequency -averaged  intensity  iAv  . 
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The  basic  idea  of  this  approach  is  due  to  Yamamoto, 
et.al.^l  and  to  Grant  and  llnnt.^2®^  However,  these  authors 
do  not  explicitly  make  the  Lcbcsgue  quadrature  identification 
embodied  in  Eqs.  (6.14)  and  (6.16),  nor  do  they  present  any 
evidence  for  the  correctness  of  those  quadratures.  An  im¬ 
portant  task  remaining  to  be  performed  is  to  investigate  the 
accuracy  of  the  frequency  quadrature  scheme  by  comparison 
with  simple  line-by-line  solutions  of  the  transfer  equation. 


6.2.4  Exponential -Sum  fitting 

Obtaining  the  exponential -sum  fit  to  the  transmission 
function  (see  Eq.  (6.12))  in  an  automatic  fashion,  with 
guaranteed  convergence  from  initial  guesses  for  the  a^  and 
,  and  for  an  arbitrary  T  ,  has  proved  to  be  a  major 
stumbling  block  in  the  code  development.  What  follows  is  a 
brief  chronicle  of  our  experience  with  several  fitting  schemes. 

The  approach  taken  was  to  minimize 


,  (ai . aM»  ki»  •  •  ’  kMJ  “ 


N. 


E  £ 


j-i 


M 


U-i 


-k.u.  12 

a.  c  J  -  T.  (u.f 
l  Avv 


(6.24) 


subject  to  the  constraints 


a.  >  0  and  k^  >  0  (i=*l,  . . . ,  M) 


(6.25) 


A  further  constraint. 


1 
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was  added  in  the  beginning  to  ensure  T^(0)  =  1  exactly, 
but  it  proved  difficult  to  handle  and  was  eventually  rejected 
in  favor  of  using  u  =  0  as  one  of  the  iij  .  The  weights 
Wj  were  originally  taken  to  be  unity,  but  the  various 
algorithms  tried  seemed  to  converge  more  rapidly  when 


wj  ' 


IT 


A^V1 


which  amounts  to  summing  relative  rather  than  absolute 
residuals. 

The  first  types  of  methods  surveyed  were  those  of  un¬ 
constrained  minimization  of  a  general  scalar  function  F:(x)  . 
Of  various  standard  methods * 30 *31  »32 *331  the  one  used  by 
most  authors  is  the  Flctcher-Powcll -Davidon  (FPD)  method. 

It  is  guaranteed  to  be  "stable, "  in  the  sense  that  the  func¬ 
tion  value  F  decreases  with  each  iteration  stop,  and  it  is 
computationally  economical . 

As  pointed  out  by  Bard, ^ 35 ^  the  method  will  fail  if 
the  ll-matrix  involved  becomes  singular,  causing  the  scavch 
to  be  restricted  to  a  sub-space  of  the  full  parameter  space. 
To  avoid  this  contingency,  our  version  of  HIM)  re-initialized 
II  to  the  identity  matrix  I  whenever  dc t (II)  <  K  ,  where  K 
was  a  pre-selected  constant  (11=1  gives  the  method  of  steepest 
descents).  Nevertheless,  the  suspicion  remains  that  the 
failure  of  1-1*1)  for  some  of  our  test  problems  is  due  to  the 
near-singularity  of  II  . 

Since  the  FPI»  method  must  be  unconstrained,  the  con¬ 
straints  of  lit|.  (6.25)  were  incorporated  into  a  transforma¬ 
tion  (suggested  by  Box^3^) 


a.  =  sin' 


k.  =  s? 

i  l 
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ill  which  y^  and  arc  new  unconstrained  parameters.  (Ail 

alternative  would  be  to  use  the  constrained  version  of  I-'Pl) 
developed  by  Goldfarb,^^  hut  the  complexity  of  Goldfarb's 
method  hardly  seemed  warranted  by  the  simple  constraints 
ai»ki  >  0  .) 

Two  versions  of  FPI)  were  tested;  one  used  cubic 
interpolation  for  the  one-dimensional  minimization  (as  sug¬ 
gested  by  Fletcher  and  Fowcll  ^4^)  t  and  the  other  employed 
parabolic  interpolation.  The  two  versions  performed  similarly 
on  several  test  problems  in  which  a  4- term  exponential  fit 
was  required  for  the  functions 

(I)  TAv(u)  -  0.4c'6u  ♦  0.3c"4u  ♦  0 . 2e  u  ♦  0.1c‘0,25u 
(II)  TAv(u)  - 

1..3/2 

(111)  TAv(u)  *  0.7c"J/u  ♦  0.3c'  7U 


Problem  II  represents  the  strong-line  limit  of  the  Goo<l)y  band 
model  transmission  function. ^ 24 ^ 

The  FPD  method  (and  the  non-linear  least  squares 
methods  discussed  later)  requires  initial  guesses  for  the  a. 
and  k^  .  Our  experience  shows  that,  unless  these  initial 
guesses  arc  fairly  good  (that  is,  F  is  small  initially) 

FPI)  docs  not  converge,  or  converges  to  an  unsatisfactory  local 
minimum.  Therefore,  some  effort  was  made  to  devise  a  satis¬ 
factory  initialization  procedure.  The  best  one,  as  judged 
by  its  ability  to  approximate  the  coefficients  and  exponents 
of  the  l'Av(u)  of  problem  1,  involved  "peeling  off"  a's  and 
k's  starting  from  the  largest  Uj's  and  working  downward. 

If  0  *  Uj  <  u,  <  . . .  u..  ,  then  presume  u„  ,  is  large 
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enough  that  one  term  e"^lu 
u  >  uN  .  Knowing  l'Av ( ux  ) 

us  to  solve  for  a^  and  . 
TAv(u)  ,  we  may  now  peel  off  a 
Checks  must  be  made  that 


dominates  the  sum  for 

and  TAv(u^  )  then  enable 

^  *  k  it 

Subtracting  a^  c  1  from 
2  and  k2  ;  and  so  forth. 


(u) 


-k^u 


>  0 


and 


E 


ai 


< 


1 


at  every  step  I  ;  when  either  criterion  first  fails  to  he 
satisfied',  the  procedure  must  be  terminated.  The  remaining 
unevaluated  a's  may  be  distributed  uniformly  in  order  to 
satisfy 

M 

a.  s  1  and  a^  >  0 

1 

The  remaining  k's  may  be  distributed  uniformly  between  the 
last  evaluated  k  and  km.|X  *  kVJ  .  To  estimate  km{|X, 
assume  u0  is  small  enough  that  we  may  approximate  every  ex 
ponential  by  its  two-term  series  expansion: 
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r(u2) 


‘kiu2 


ai 


(1  -  kju2) 


3  1 


u. 


aMkM 


whereupon 


1  -  T(u2) 
U2aM 


Using  the  above  initialization,  the  FP1)  method  con¬ 
verges  for  test  problems  1  and  11,  and  diverges  for  test 
problem  111.  The  divergence  indicated  that  FPU  would  be 
unsatisfactory  for  fitting  general  transmission  functions. 


To  find  an  alternative  to  FPU ,  a  further  literature 
search  was  undertaken  into  methods  which  capitalized  on  the 
particular  form  of  I*  in  bq.  (6.24),  the  so-called  non- 1  incar 
least  squares  (J.S)  methods.  The  simplest  of  these,  usually 
called  the  Gauss -Newton  method^ was  programmed.  A  calcu¬ 
lation  of  problem  1,  resulted  iu  rapid  divergence.  The  fail- 
uie  was  due  to  a  gross  overshoot  in  estimating  the  position 
of  the  one-dimensional  minimum  on  one  of  the  iterations. 
Further  literature  search  indicated  overshooting  is  a  rather 
common  problem  in  LS  fitting,  and  that  l.evcnbcrg^ 3!>J  and 
later  Marquardt^0!  have  devised  remedies.  The  subroutine 
CNL1.S,  which  is  a  highly  sophisticated  I.S  fitter  using  a 
modified  Levenberg  method,  was  therefore  adapted  from  the 
Oak  Ridge  numerical  analysis  library^11*  and  applied  to  the 
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test  problems.  The  results  were  the  same  as  with  l-TI),  ex¬ 
cept  that  convergence  was  much  faster  on  problems  1  ami  II. 
Divergence  was  again  experienced  on  problem  III. 

There  remain  three  possibilities  for  a  successful  LS 
exponential- sum  fitting  routine: 

(1)  a  subroutine  of  Lang  and  MKllcr^42^  which 
has  been  requested  from  a  program  library 
in  Ireland  and  which  has  still  not  arrived; 
it  purports  to  be  able  to  do  exponential- 
sum  fits; 

(2)  the  Marquardt  subroutine  from  Ref.  (41), 
which  is  currently  being  implemented; 
and 

(3)  CXLLS  run  in  its  more  sophisticated 
mode,  in  which  pre-spccif icd  sub-sets  of 
the  parameters  are  held  invariant  for  a 
pro-specified  number  of  iterations  or 
until  convergence  for  the  active  sub- set 
is  attained. 


Particular  interest  attaches  to  possibility  (1)  because  the 
authors  claim  that  their  subroutine  suppresses  divergences 
of  the  type  we  encountered  in  problem  III. 


However,  we  felt  that  we  could  not  rely  exclusively 

on  an  LS  method.  Therefore,  a  literature  search  for  non-I.S 

methods  of  exponential -sum  fitting  was  undertaken,  and  a 

substantial  number  of  algorithms  were  found^5,<*  *  >45,40,47,48] 
1 47) 

Lanezos1  1  discusses  the  fundamental  problem  in  exponential 
approximation,  which  is  that  the  exponentials  are  highly 
non-orthogonal  functions.  This  implies  that  algorithms  such 
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as  Hudson's ,  Rice's,  ^nd  Prony's^^  are  extremely 

ill-conditioned  numerical ly .  The  values  of  the  a.  and  k. 
turn  out  to  be  extremely  sensitive  to  the  data  T/\v (uj )  ; 
small  changes  in  the  T^  ' s  »  even  in  the  fifth  or  sixth 
decimal  place,  can  produce  large  changes  in  the  a^  and  k.  . 
This  is  a  consequence  of  the  fact  that  there  are  many  fits, 
with  widely  varying  a^'s,  k^'s,  and  M's  which  can  approxi¬ 
mate  a  given  to  the  accuracy  to  which  it  is  known. 

Lanczos  and  Hildebrand give  examples  of  this  phenomenon. 

This  non-uniqueness  raises  the  question  of  whether 
any  physical  interpretation  can  be  attached  to  the  k^  or 
whether  the  whole  procedure  of  Sections  6.2.2  and  6.2.3 
can  only  he  judged  by  the  correctness  of  i ^  ,  the  intervening 
steps  being  regarded  purely  as  a  mathematical  device.  The 
comments  of  Lanczos  together  with  our  own  experience  cause 
us  to  lean  to  the  latter  interpretation.  It  then  also  be¬ 
comes  incumbent  that  any  two  equally  good  fits  of  T^ 
should  lead  to  the  same  i^  ,  a  requirement  that  remains 
to  be  checked. 

The  available  algorithms  divide  rather  naturally  into 
those  with  prc-specified  ^143,45,46,47]  aruj  t]losc  in  wlri cli 
the  algorithm  determines  jq,  [44,48]  Since  there  is  no  a  priori 
basis  on  which  to  choose  M,  having  it  determined  by  t lie 
algorithm  seems  to  be  a  desirable  property.  Gardner's 
method,  which  requires  a  numerical  Fourier  transform 
g^v(k)  of  T^v ,  would  lie  extremely  difficult  to  program,  since 
it  involves  a  good  deal  of  adjusting  of  the  limits  of  integra¬ 
tion  of  two  infinite  integrals  in  order  to  distinguish 
spurious  peaks  in  (k)  from  real  peaks  (representing  the 
k|).  As  for  the  fixed-M  methods,  they  suffer  a  further  dis¬ 
advantage  in  that  they  all  solve  for  the  k^  as  roots  of  a 
polynomial  whose  coefficients  are  determined  by  a  generally 
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ill-conditioned  linear  least  squares  procedure.  The  instabil¬ 
ity  of  polynomial  roots  to  coefficient  variations  is  well- 
known,  and  hence  there  is  a  substantial  possibility  of  grossly 
mis -estimating  the  ki  in  all  the  fixcd-M  methods.  Further¬ 
more,  there  is  nothing  to  prevent  negative  or  complex  k.  . 
Thus,  these  methods  are  poor  choices. 

The  remaining  candidate  is  the  Cantor  and  Evans 
method ^  ^  which  has  several  advantageous  properties  in  addi¬ 

tion  to  determining  its  own  M.  Of  these  the  most  important 
is  unconditional  convergence  (to  arbitrary  accuracy).  In 
addition,  the  size  of  M  is  controlled  by  the  desired  accu¬ 
racy  and  >  0  and  k^  >  0  .  A  disadvantage  is  that  the 
Uj  must  be  equally  spaced.  But,  all  things  considered, 
the  method  of  Cantor  and  Evans  seems  eminently  suited  to  onr 
problem . 

Professor  Evans,  who  is  at  the  University  of  Califor¬ 
nia  at  San  Diego,  has  agreed  to  consult  with  us  on  his 
method.  Wc  anticipate  a  working  version  shortly. 


6.3  THE  MONOCHROMATIC  EQUATION  OE  TUANS PER 

The  method  of  solution  of  the  radiative  transfer 
equation  due  to  Grant  and  Hunt  ^  ^ ^ ’  ■'’*'*  ^ ^  has  been  chosen 
to  solve  the  monochromatic  problems  arising  in  Section  6.2. 

Grant  and  Hunt's  procedure  requires  no  iteration  in  the  pres¬ 
ence  of  scattering,  allows  zones  of  arbitrary  size,  and  is 
computationally  economical.  It  brings  together  several  lines 
of  research  in  radiative  transfer.  On  the  one  hand,  it  repre¬ 
sents,  for  1-1)  plane  parallel  problems,  the  culmination  of 
the  discrete  space  formulation  of  transfer  theory  associated 
with  Prciscndorfcr^ ’t’5l  and  van  de  Mulst.^*’^  On  the  other 
hand,  it  reduces  in  the  limit  of  small  zone  size  to  the  dif¬ 
ferential  formulation  of  Rybieki  and  Usher  J*r,7l  which 
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represented  nit  elegant  formulation  of  the*  method  of  i nvnri ai.t: 
imbedding.  l.astly,  it  improves  upon  the  Sn  melhoil 1 ^  (by 
eliminating  the  need  for  zone-centered  intensities)  to  the 
point  where  previous  Sn  methods  for  plane  parallel  geom¬ 
etry  are  now  virtually  obsolete. 

The  elimination  of  scattering  iteration  is  a  very 
significant  feature  of  the  Grant  and  Hunt  method.  To 
iterate  in  the  previously  proposed  fashioiJ  1 1  when  optically 
thick  clouds  are  present  would  have  been  prohibitively  expen¬ 
sive.  Other  attractive  features  of  the  method  are: 

(1)  a  single  inequality  imposed  on  the 
"primary  layer  thickness"  At](  (cf. 

•Section  6.5.1)  assures  both  positiv¬ 
ity  of  intensities  and  computational 
stahil i ty ; 

(2)  flux  is  conserved  provided  only 
that  the  phase  function  Fv  is  cor¬ 
rectly  normalized  (cf.  Section  6.4); 
mid 

(3)  error  estimates  are  available. 

More  generally  the  Grant  ami  Hunt  method  possesses  a  certain 
elegance,  simplicity,  and  freedom  from  ad  hoc  assumptions 
which  alternative  methods  Jack. 

The  basis  of  the  Grant  and  Hunt  method  is  the  inter¬ 
action  principle.  1  1 
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u+(t2)  -  t(r2,T1)  u+(i1)  +  r^,!,,)  u“(t2)  +  E+(Ti*T2^ 

(6.26) 

u'Ctj)  =  r ( t  2 »  t  ^ )  U+(T1)  +  t(r1,T2)  u‘(x2)  +  e'(t1,t2) 

Figure  6.1  illustrates  the  various  quantities  involved 
in  describing  the  interactions  in  the  slab  between  and 

x2  .  The  total  optical  depth  coordinate  is  t  and  p  =  cosG 
is  the  cosine  of  the  angle  between  the  ray  and  the  vertical. 
The  u  's  are  positively  directed  (p  >  0)  intensities,  the 
u  's  negatively  directed  (p  <  0)  intensities  (the  p-depcndcnce 
is  implicit).  The  r's  and  t's  are  integral  operators 
representing  reflection  from  and  transmission  through  the  slab 
[xltT2],  and  the  Z's  represent  effective  internal  sources  in 
[Ti»t2].  Wc  now  approximate  the  angular  integrals  in  the 
interaction  principle  by  naif-range  Gaussian  or  Radau  quadra¬ 
tures  with  weights  and  angles  0  <  p2  <  ...  <  p  <  l 

(Pm  a  3  for  Radau  quadrature).  (By  using  half-range  quadra¬ 
ture,  an  idea  originally  due  to  Krook,^59^  in  which  separate 
but  symmetric  quadrature  formulas  arc  used  for  0  <  p  <  1 
and  - 1  £  p  <_  0  ,  discontinuities  in  the  intensity  across 
p  =  0  (0=90°)  which  may  occur  physically  arc  allowed  for.) 

The  r's  and  t's  become  matrix  operators,  and  the  u's 
and  E's  become  vectors: 

’^(Tj  ,t2  ,±pj) 

^I('fl>T2)  = 

(6.27) 

where  iv  has  been  written  as  a  function  of  t  rather  than 
liq .  (6.26)  will  now  be  taken  to  be  in  this  matrix-vector  form. 
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l*oi  the  sake  of  computational  economy,  we  intend  to 
use  only  a  relatively  few  (m  <  7)  angles.  Therefore,  the  de¬ 
sirability  of  solving  for  the  diffuse  intensity  i  rather 
the  full  intensity  I v  ,  which  includes  the  solar  flux,  is 
apparent;  if  we  used  ,  many  more  angles  \i ^  would  he 
needed  to  resolve  the  solar  peak,  and  these  could  no  longer 
be  distributed  in  the  optimal  (Gaussian)  way.  We  are,  of 
course,  ignoring  a  possible  specular  reflection  peak  in  i  , 
which  has  n£t  been  subtracted  out.  (This  would  be  of  concern 
only  over  fairly  smooth  water  surfaces.)  Tests  will  he  per¬ 
formed  for  a  variety  of  situations  to  ensure  that  the  angular 
discretization  is  adequate  to  determine  fluxes  accurately. 

6.3.1  Layer  Composition 

We  now  consider  how  r,  t,  E  for  each  zone  arc  to  he 
obtained.  Lach  zone  is  divided  into  sub-layers  called  pri- 
nlar>'  layers,  which  arc  optically  thin  enough  that  r,  t,  X  for 
each  primary  layer  can  be  obtained  from  the  physical  quanti- 
tics  Bv,  l’v,  etc.  of  the  differential  formulation,  Eq.  (6.3). 
Then  the  primary- layer  r’s,  t’s,  Z»s  are  "added  up,"  by  a 
process  wo  shall  call  layer  composition,  to  obtain  r’s,  t’s, 
and  E’s  for  the  entire  zone. 

The  appropriate  primary- layer  quantities,  derived  by 
comparing  the  limits  as  At  =  t2  -  Tj  ->  0  of  Eqs.  (6.26) 
with  Eq.  (6.3),  are[50^ 

t(T2’T.l)  =  t(Tl’T2}  =  1  "  M _1[J  ’  P++(C)c]at  +  0 ( Ax 2 ) 

r(T2’Tl)  =  r ^T1  ’ t2-)  =  M"1  ,j+'CC)cAt  +  0 (At 2 )  (6.28) 

^±(Ti>t2)  =  M'1  B±(C)At  +  0 (At2 ) 
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To  O(At),  t(T2»Tj)  and  t(Tj»T2)  arc  identical;  however,  they 
nay  differ  in  0(Atz)  terns.  A  similar  remark  applies  to  r. 
The  various  scalars,  matrices,  and  vectors  occurring  in  l;qs. 
(6.28)  are  defined  ;;s 


M  “  T2  ‘  T;  »  Ti  -  ^  -  t2  * 

M  “  lpi  6ijJ  »  C  *  lc.  6i.]  , 


-  f**v(T ,-iii#-Hj) J 


1,+  ‘(T) 


W(T) 


Bv(t) 

s(?) 


(6.29) 


B(t,w)  a 


U  *  w(T)J  Bv(T(t))  ♦ 


Sv  «(T) 
‘fi  " 


VT-"’"SU„ 


)e'T/,,-n 


B4(t) 


IKt.iPj) 


It  is  worthwhile  to  note  that  the  error  terms  in  Hcjs.  (6.2S) 
will,  in  fact,  be  0(Ai3)  if  £  is  chosen  as  the  inid-point 
of  l T i * T 2 J  follows  from  Taylor's  theorem. 

In  order  to  guarantee  positive  intensities,  all  ele¬ 
ments  of  the  primary- layer  r-  and  t-matriccs  in  l-q.  (6.28) 
must  be  positive.  Ignoring  0(Atz)  terms,  this  reduces  to 
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the  requirement  that  the  diagonal  elements  of  t  be  positive, 
which  gives  us  the  following  restriction  on  At  : 


/lT  <  AVi*  "  miM 
ajx  l<i<m 


1  -  Wci  V'-'V'P 


(6.30) 


It  is  also  true  that  this  property  of  positive  elements  in 
r  and  t  is  preserved  by  layer  composition,^^  so  that  no 
restrictions  beyond  liq.  (6.30)  are  necessary. 

The  question  of  truncation  error,  thal  is,  the  error 
incurred  by  neglecting  0(At2)  terms  in  Hqs.  (6.28),  is  also 
of  concern.  Grant  and  Hunt  have  used  as  a  practical  critc* 
rio„l5‘J 


At 


1 

2 


At 


max 


Grant*  has  also  made  some  estimates  of  the  growth  of  trunca¬ 
tion  error  during  the  process  of  layer  composition,  lie  finds 

that  if  T  and  R  arc  the  correct  transmission  and  ro¬ 
ll  n 

flection  matrices  for  a  layer  of  thickness  iiAt  ,  and  if 
this  is  built  up  by  composing  primary  layers  of  thickness 
At  ,  then  the  effect  of  small  truncation  errors  n  in  Tj 
and  p  in  R^  is  to  give  approximate  matrices  tn  ,  rn 
such  that 


I  I  »n  -  tj|  <  c:nAn(||p||  ♦  |  |n||) 

MR,  *  rn 1 1  <  2,1(1  |P|  I  ♦  |  |n|  |) 


(6.31) 
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where  C  and  A  are  defined  such  that 

lltjl  <  a®  (u  <  *  <  i) 

and  ||  ||  is  the  norm  of  Ref.  (51).  Thus,  for  fixed  At  , 

the  absolute  error  in  It  and  the  relative  error  in  T 

n  n 

grow  about  linearly  with  n  .  If  we  use  an  approximation  in 
which 


I  Ml  *  I  ln|  |  -  0(Atatl) 


and  write  n  =  I  I/At  for  fixed  II,  then 

I  | T(ll)  -  t (II)  |  |  <  KCII  An  At“ 
|  1 11(11)  -  r(II)  |  |  <  2KII  Axa 


(6 . 52) 


so  that  the  errors  tend  to  zero  as  At  -►  0  .  Currently,  we 
intend  to  ignore  0(At2)  terms  in  l;qs.  (6.28),  so  that 
a  ■=  1  (or  a  =  2  if  f,  is  taken  as  the  mid-point).  If  time 
permits,  however,  we  may  investigate  the  advantages  of  carry¬ 
ing  more  terms  in  the  expansion. 

The  layer  composition  formulas  are  derived  as  follows. 
Suppose  t.  <  t.  <  .  Write  down  the  interaction  principle 

(6.26)  Cor  the  layer  I T i , T 2 J  >  ;,,u*  again  for  the  layer 
|t2,t5J  .  eliminate  ii+(t,)  and  u“(t,)  from  two  of  these 
equations  using  the  remaining  two.  Put  the  first  two  equa¬ 
tions  in  the  form  of  an  interaction  principle  for  the  layer 
ITj,t5)  .  Then  identify  i*j .  r(Tj,T-)  ,  etc.  from  this 
principle.  The  results  are 
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who  rc 


r13  =  r23  *  l52r  r12  *23 
r31  '  r2i  *  ‘ijf  r32  *21 


*31  ’ 

*32'' 

*21 

*13  = 

*12*' 

*25 

*32r 

(r12 

SM  * 

♦ 

4  *23 

13 

*12f 

<r52 

S12  ‘ 

X23) 

4  £n 

r  ■  ('  ■  r,2  r^)'1 


(6.33) 


(6.34) 


In  the  special  case  when  the  layer  [t^Tj]  is 
homogeneous .  (i.e.,  when  w(x)f  I,+  +(t),  and  1,+  "(t)  do  not  vary 
significantly  across  (t^t-J)  and  if  the  two  layers 
and  It2,t5J  arc  of  equal  thickness  At  ,  then  formulas  (6.33) 
slrapl i fy  somewhat : 


rl  -  1 15  ’  *51  ‘  ro  *  *o  ro  ro  *< 


*13 


t..  -  t  r  t 

.>1  o  o  o 


Y. *  £  Z* 


1  ^13  8  *o  ro<ro  >:o  +  V  +  *o 


y]  2  >:'  =  t  r  (r  r.*  ♦  r.~ )  ♦  >;" 

1  1.*  o  o v  o  o  o'  o 


(6. 35) 
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where 


and 


ro  -  *  ro  ro> 


-1 


ro  :  r12  =  r21  “  r23  a  r32 


lo  5  l12  “  l21  =  *23  =  *32 


Eo  "  £12  =  ^23 


(6.36) 


Relations  (6.56)  follow  from  Cqs.  (6.28)  by  the  assumptions 
of  homogeneity  and  equal  Ax’s  .  The  single-subscript  nota¬ 
tion  indicates  that  r,  t,  l  depend  on  only  one  argument 
instead  of  two,  that  argument  being  the  corresponding  layer 
thickness.  Defining 

=  r(2n  At)  ,  etc. 

it  is  possible  to  continue  the  above  process,  called  doubling , 
through  any  homogeneous  region, 


t 


r 


11+ 1 

tn 

‘n 

11+ 1 

rn 

*  rn  rn 

*n 

+ 

11+1 

tn 

Vrn  sn  * 

*  < 

11+1 

tn 

*  £,'. 

(6.37) 
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Doubling  is  much  faster 
in  liqs ,  (6.53).  To  add 
gcncotus  layer  requires 
doublings . 


than 


simple  layer  addition  as  embodied 
primary  layers  comprising  a  homo- 
1  layer  additions  versus  only  N 


In  the  future  the  code  will  allow  for  arbitrarily 
larger  zones,  which  may  or  may  not  be  homogeneous.  If  a  zone 
is  homogeneous,  doubling  will  be  used;  if  it  is  not,  the  code 
will  partition  it  into  sub-zones  which  are  homogeneous,  double 
within  cacli  of  the  sub -zones,  and  add  the  sub- zones  using 
Eqs.  (6.33).  For  now,  however,  the  code  will  be  constructed 
on  the  assumption  that  each  zone  is  homogeneous.  This  is 
fairly  restrictive,  since  it  implies  that  neither  the  scatter¬ 
ing  nor  the  absorption  coefficient  vary  substantially  across 
a  zone.  No  more  than  40  zones  should  be  necessary,  however, 
and  this  docs  not  seem  overly  burdensome  (Grant  and  Hunt  run 
witli  about  40  zones*). 


6.3.2  Source  Doubling 

The  problem  of  an  inhomogeneous  source  2  4  in  a  zone 
which  is  otherwise  homogeneous  requires  special  attention. 

An  inhomogeneous  source  region  will  occur  frequently  on  ac¬ 
count  of  the  factor  c  T^sun  in  B(t,jj)  and  on  account  of 
the  rapid  variation  of  BV(T)  with  T  in  many  spectral  re¬ 
gions.  It  turns  out,  however,  that  the  most  frequently  occnr- 
ing  types  of  source  inhomogcncity  still  admit  a  modified  form 
of  doubling.  Inhomogeneous  sources  can  be  doubled  right  along 
with  r  and  t  with  the  resulting  saving  of  computing  time. 

Let  us  begin  by  developing  a  general  formula  for 
source  composition  when  the  source  variation  across  an  other¬ 
wise  homogeneous  zone  is  described  by  f(-r)  .  Referring  to 

* 

Private  communication. 
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Figure  6.2  and  Hqs.  (6.28), 
of  the  primary  layers  he 


let  the  source  vector  for  any  one 


£•  a  4  C(C.) 


i  =  1, 


,N 


(6.38) 


where  2"  arc  independent  of  t  ,  and  is  the  mid-point 

of  primary  layer  i  , 

fi  ■  >0  '  -  r-  (0.39J 

Adding  the  sources  in  layers  1  and  2  according  to  Uqs.  (6.33), 
it  is  clear  that  the  source  EJ  +  2  ^or  ^lc  ined  layer 
will  be  a  linear  combination  of  f(£  )  and  f(£  ): 


Zone 


Figure  6.2.  Homogeneous  zone  composed  of  2^ 
primary  layers. 
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,;1.2  *  ''I0’  +  V2a)  C«2> 

(The  expressions  for  the  V’s  arc  not  of  interest  here.) 
Similarly,  the  source  for  the  combined  layer  3+4  involves  only 
a  change  in  the  S’s: 


v*(i) 

1 


fU5)  ♦  v 


±(D 

2 


fu4) 


Then,  upon  combining  layers  1+2  and  3+4,  a  linear  combination 
of  f(Cj)  through  f(C4)  results: 

^ J  +  2+3+4  -  T  VJ'2)  f(V 
1=1 

Suppose  that  any  2n  adjacent  primary  layers  have  been  super¬ 
posed  in  the  above  manner.  Using  the  notation 

+  + 

Z(i,n)  =  i  +  (i  +  1) +  ...  +(i  +  2n-l) 

where  i  denotes  the  first  primary  layer  in  the  grouping, 
it  is  possible  to  prove  by  induction  that 

2n 

*  E.'l00  f'Wl)  (6'40) 

1=1 

h’c  now  specialize  our  considerations  to  the  solar 
source  (f(t)  -  e  T^,Jsun)  and  the  Planck  source 
(f(t)  =  B  (T(t))J  .  By  the  linearity  of  the  source  composi¬ 
tion  formulas  (6.53)  it  is  possible  to  compose  the  solar 


152 


3SR-1034 


source  and  Planck  source  separately  and  then  simply  add  the 
results  at  the  end  to  obtain  the  full  source. 


For  the  case  of  the  solar  source,  the  quantities 
£p  of  Eq.  (6.38)  are,  from  Eqs.  (6.28)  and  (6.29), 


ZP,sol 


Sv  w(t) 
4tt 


At 


V(T  ,±tJm,tJsun^-i 


(6.41) 


Since  w  and  P  are  assumed  not  to  vary  across  the  zone, 
the  argument  t  refers  to  the  mid-point  of  the  entire  zone 
in  Figure  6.2.  Now  assume  that,  beginning  with  the  E*  gol 
and  Eq.  (6.38),  the  first  2n  primary  layers  have  been  com- 
posed,  so  that  the  are  known.  The  adjacent  2U  pri¬ 

mary  layers  will  then,  by  Eq.  (6.40),  have  source  vectors 


>n 


E 


yiCll)  £/■£-  n-j 

(l  +  2n  ,n)  ft  *  (^+2) 


v±(n)  -(t  +  (£  +  2n)Ar  -  ~)/d 

v  o  v  J  2  J  sun 


i-1 


'Z"4T/"s„n  * 

bl.n) 


(6.42) 


h 


n 


Z(l,n) 
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where 


h 


n 


At/m 


sun 


(6.45) 


Remembering  the  definitions  rn  =  r(2n  At)  ,  etc.,  the  sources 

E(l,n+1)  ^or  t*lc  comkincd  layer,  containing  2n+^  primary 
layers,  arc 


£(l,n*l) 


£(1 »n+l) 


ln  rn(rn  Z(l*2%)  +  SU.n))  +  ^l  +  2\n) 

ln  rn(Jn  hn  Z(l,n)  +  Z(l,n))  +  hn  X(l,n) 

tn  rn(rn  E(l,n)  +  r'(l  +  2n,n))  +  E(l,n) 

tn  rn(rn  E(l,n)  +  hn  E(l,n))  +  2(l,n) 


(6.44) 


Ihese  are  adaptations  of  the  source  composition  formulas  in 
hcjs .  (6,33).  From  the  definition  of  h  ,  it  is  obvious  that 

hn.l  =  hn  (G.4S) 


logether ,  Eqs .  (6.44)  and  (6,45)  form  a  doubling  scheme  for 
the  solar  (or  any  exponential)  source.  The  scheme  is  initia¬ 
lized  by 


Z(l,0) 


,sol 


+  ^-)  /  M  . 

•  sun 


2 


h  = 


-At/m 


sun 


(6.46) 
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and  iterated  from  n  «  0  to  n  =  -  j 

In  order  to  derive  doubling  formulas  for  the  Planck 
source,  it  is  necessary  to  have  a  more  tractable  functional 
Jo.'iu  for  f(r)  than  Bv(T(t))  .  In  Appendix  B  it  is  shown 
that  it  is  reasonable  to  assume  B v  is  piecewise  linear  in 
t  ,  so  that 


f(T)  a  Bv(T(t))  =  Bq  + 


(6.47) 


across  the  zone  of  figure  6.2,  where 

Bo  5  bv(T(t0)) 

B  •  -  •  By(TfT0)) 

T1  ■  T0  ‘ 

Expressing  this  in  terms  of  £.  (see  l-q.  (6.39)), 


f(CA)  -  Bq  +  B'(t  +  iAx  -  41 


-  B^  +  B '  iAi 


(6.48) 


where 


K  -  l!o  -  11 '  T- 


Ihc  sources  Sfj  >nj  for  the  first  2n  primary  layers  then  be 


come 


»n 


Z 


r4(n) 


Cl.n)  =  £  (B'  *  n't  At) 

£=1 
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while  for  the  adjacent  2N  primary  layers 

2n 

£(l»2n,n)  “  £  Vi0>)  i»i  *  B'(l*2")4tl 
ft=l 


"  S(l.n) 


+  2n  B'  At  Y4 
n 


(6.49) 


where 


A=  1 


The  Y~  refer  to  a  constant  source  f(x)  *  1  ,  and  so  may  be 
calculated  by  doubling  as  in  liqs ,  (6,37).  A  further  simpli¬ 
fication  is  that  symmetry  prevails, 


Y 


n 


This  is  intuitively  obvious,  since  a  constant  source  would  be 
expected  to  emit  the  same  at  -p  as  at  +  p  .  It  can  be 
proven  inductively  by  first  noting  that  the  I*  of  liq .  (6.38), 
which  for  the  Planck  source  are, 


ZP,plk 


£l  -  w(t)Jat  M"1 


=  [l  -  <*>(T)j 


At 


LlJ 

1/Mj 

}/vm. 


(6.50) 
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j 


arc  equals 


El»,plk  5ZP,plk  “  El>,pik 


Because  f(i)  =  1  , 


^P.plk 


(6.51) 


by  Eq.  (6.38),  so  that  symmetry  is  proved  (or  n  »  0  .  If 
symmetry  has  been  proved  for  a  general  n  ,  then  from  the 
doubling  formulas  (6.37) 


n+1 


n 


r  (r  Y 
n v  n  n 


♦  Y  ) 
nJ 


♦  Y 


n 


t  T(r  Y  ♦  Y  )  ♦  Y 
n  nv  n  n  nJ  n 


which  completes  the  induction  proof,  and  gives  us  the  simpli¬ 
fied  doubling  formula 

Yn*l  ■!>♦*■  r«C  ♦  rn»v»  ‘6'52> 


f°r  Y„  . 

The  two  layers,  each  containing  2n  primary  layers, 
arc  now  composed  just  as  for  the  exponential  source  (see 
Uqs.  (6.44)) : 


137 


3SR-1034 


'(l.n+l)  ‘  tn  rn[r„(£a.n)  *  2"  «’  V  *  za,„)] 


+  St.  n,  +  2n  B'  At  Y 
U.n)  n 


£(l,n+l)  =  tn  rn[rn  Z(l, 


n)  ^ (1 »n) 


+  2“  B'  Ax  Y  +2 


(6.53) 

(l>n) 


Together ,  Eqs .  (6.52)  and  (6.53)  form  a  doubling  scheme  for 
the  Planck  (and  any  linear- in -t)  source.  The  scheme  is 
initialized  by  Eq.  (6.51)  and  by 


(1.0) 


=  («'  +  B*  At)Y 


(6.54) 


and  iterated  from  n  =  0  to  n  =  N  -  1  . 

More  accurate  piecewise  polynomial  approximations  to 
»v(T(t))  arc  possible,  at  the  expense  of  more  complicated 
source  doubling  formulas.  In  particular,  economical  piece¬ 
wise  quadratic  and  piecewise  cubic  schemes  have  been  derived 
and  will  be  programmed  at  a  future  date  to  ascertain  if  the 
increased  accuracy  of  the  heating-rate  calculations  so  ob¬ 
tained  warrants  the  additional  computational  expense. 

6.3.3  Forward  and  backward  Passes 

By  the  procedures  of  sections  6.3.1  and  6.3.2,  we 
obtain  reflection  and  transmission  matrices  and  some  vectors 
for  every  2011c  of  the  atmosphere.  If  the  zone  structure  is 


0  .  t,  .  t2  ,  <  V] 


and  the  notation  is 
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+ 

u“ 

n 

=  u±(t  ) 

v  ^nJ 

n,n±l 

5  r(Tn-Tn±P 

n,n±l 

5  tfVTn±l> 

+ 

n,n+l 

*  E±(Wl 

(6.55) 


then  according  to  the  Grant  and  Hunt  algorithm  (whose  deriva¬ 
tion  is  lengthy  and  will  not  be  repeated  here)  the  solution 
for  the  intensity  vectors  u~  splits  into  two  parts. 

The  first  part  is  the  forward  pass.  In  the  forward 

p,ass,  the  mxm  matrices  E  ,  G  ,  and  II  and  the  m-vcctors 
(1)  (21  (31  n  n  n 

Vn  '  n  ^  ’  antl  n  '  arc  calculate^  recursively  as  follows 

Ej  =  0  v}1’  »  v'2)  =  v'3)  =  0  (6.56a) 


n+1  = 

t1  -  rn+l  ,n  V 

(6.56b) 

% 

n+1 

^n+1  ^’njn+l 

(6.56c) 

'n+1  = 

E  H 
n  n+1 

(6 . 56d) 

;n+l  = 

rn,n+l  +  *n+l,n  ^n+1 

(6 . 56e) 

-Cl)  - 

n+1 

?  (r  V(3^  +  y~  1 

‘n+1  ln+l,n  vn  Ln,n+lj 

(6 . 56f ) 
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y(2)  =  y(3)  +  y(l) 

n+1  n  Si  *n+l 

(6 . 56g) 

y(3)  a  y(2) 

n+1  n+l,n  vn+l  Si, n+1 

(6 . 56h) 

(n  3  1,2,  . ..,  ,\’) 

The  initialization  of  at  zero  reflects  the  fact  that 

there  is  no  incident  diffuse  radiation  at  the  top  of  the  at¬ 
mosphere,  t  -  0  (the  incident  solar  radiation  is  already  ac¬ 
counted  for  by  the  transformation  from  to  i  ).  Economy 
of  storage  is  achieved  by  performing  the  sub- layer  superposi¬ 
tions,  which  yield  the  r's,  t's,  and  Z's,  simultaneously 
with  the  forward  pass;  in  this  way  the  r's,  t's,  and  Z's 
need  not  be  saved  as  a  function  of  zone. 

In  order  to  identify  certain  quantities  in  the  back¬ 
ward  pass  formulas,  it  is  necessary  to  look  at  the  surface 
boundary  condition  as  derived  in  the  previous  report: I1! 

iv^TN+l,p^  “  *  2  f0 

(6.57) 

(TN+1  '  )  dp'  (-1  £  p  <  0) 


where 


V">  “  «v(»)  W  *  Ms„„  P„  ("sun*") 


Applying  our  half- range  Gaussian  quadrature  to  tlie  integral 
in  liq .  (6.57), 
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ill 


iv(TN+l  *  _|,i ^ 


fi»<"i)  *  2  £  cj  Mj  MVi/j) 


(i  -  1,  ...»  m) 

This  can  be  written  in  matrix- vector  form,  using  the  notation 
of  l:q.  (6.S5)  ,  as 


Vl  "  w  +  rG  Vl 


where 


w 


W 


(6.58) 


(rCJij  •  2ej  Vj  Pv(MJ,Mi) 


The  backward  pass  may  now  be  written 


u 


N^l 


(I  -  -V.)'1  »  *  v'JJ) 


'N+l  G'  K"H*1 

Vl  "  W  4  rG  UN+1 


(6.59) 


and 


141 


3SR- 1034 


%  3  (lnU  Vl  ’  V,(,M 
“n  *  Vl  Vl  *  'nil 


(6.60) 


n  *  N.N-l ,  . . . ,  1 

Fl-om  the  components  of  the  u*  ,  the  fluxes  arc 

obtained: 


Fv(tn)  ■  2 it 


/)  "W 


*Tn/0-„n 

»0<lR  ♦  Sc 


sun  v 


2*  f  lwiv(Tn,p)  -  Miv(Tn,-M)Jdp 


♦  u  S  e 
Msun  v 


Tn^wsun 


HI 

U  £  ci  ui l'v^Tn* V  ’  MWl 


i«l 


(6.61) 


♦  u  Sc 
psun  v 


Tn^sun 


This  completes  the  algorithm. 


6.4  Mil:  SC  ATT  1;  RING  TRliATMiiNT 

Additional  effort  has  been  expended  since  the  semi¬ 
annual  report ^  to  make  the  Mie  scattering  calculation  as 
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fully  automatic  and  Independent  of  ad  hoc  assumptions  as  pos¬ 
sible.  This  work  is  summarized  in  Section  6.4.1.  Secondly, 
when  the  phase  function  has  a  strong  forward  peak,  it  can  in¬ 
troduce  serious  errors  into  the  numerical  formulation  of  Sec¬ 
tion  6.3.3^**^  unless  a  considerably  larger  number  of  angles 
is  used  than  we  intend.  Therefore,  Section  6.4.2  de¬ 
scribes  a  method  for  eliminating  this  peak.  In  Section  6.4.3, 
the  azimuthal  integration  of  the  phase  function  is  treated. 

And  finally,  a  method  for  re-normalization  of  the  phase 
function  is  suggested  in  Section  6.4.4. 


6.4.1  Mic  Scattering  Computation 

Our  subroutine  (MIOSCT)  for  calculating  the  Mic 
scattering  coefficient  Gv  ^  ,  Mic  absorption  coefficient 
a  ..  ,  and  Mic  phase  function  P  ..  for  a  particle  of 
radius  a  follows  sti.nd.-.id  techniques,  l60*61*62*6-^65*66! 
with  minor  modifications  to  achieve  computational  efficiency. 
The  subroutine  was  tested  against  two  sets  of  tables ^ 

f<*  9  1 

as  well  as  against  certain  graphs  in  Kerker 
all  eases  with  the  published  results. 


and  agreed  in 


Techniques  for  the  integration  of  the  Mic  quantities 
over  a  size  distribution  n(a)  arc  considerably  less  stand¬ 
ard.  The  most  comprehensive  study  of  this  problem  is  due  to 
beirmendj ian, who  used  a  trial -and-error  method  to  dis¬ 
cover  just  how  small  the  integration  increments  Aa  ,  where 


7  w  *» 

a  -  —j~  (A  »  wavelength)  , 


needed  to  be  in  many  different  cases.  Deirnicndj ian  offered  no 
universally  applicable  method  for  picking  Act’s  small  enough 
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to  ensure  accuracy  but  large  enough  to  avoid  wasting  computer 
time.  Studies  by  Davc^^*^^  arc  of  limited  applicability 
since  they  were  restricted  to  the  visible  spectrum.  However, 
the  work  of  Dave  and  especially  that  of  Dei nnend jinn  estab¬ 
lished  certain  ground  rules  for  the  integration:  (1)  it  is 
best  to  use  trapezoidal  integration  since  the  Mic  functions 
oscillate  so  rapidly  with  a  that  higher-order  Newton-Cotes 
quadrature  schemes  lead  to  no  increase  (and  may  actually  lead 
to  a  decrease)  in  accuracy;  (2)  the  necessary  Aa  is  in¬ 
versely  proportional  to  A  ;  (3)  the  larger  a  (the  upper 
limit  of  integration) ,  the  finer  the  angular  mesh  needed  to 
properly  resolve  Tv  M  ;  and  (4)  the  larger  the  interval  of 
integration,  <*max  -  am.n  ,  the  larger  the  allowable  Aa  . 

To  these  observations,  we  added  several  of  our  own. 

First,  the  computational  expense  increases  enormously  with 
increasing  a  x  ,  due  to  the  fact  that  the  number  of  terms 
in  the  Mic  series  (summed  by  MlliSCT)  is  roughly  equal  to  a  . 
Therefore,  it  was  decided  to  split  the  size  integration  into 
two  integrations,  the  first  of  which  is  for  the  purpose  of  re¬ 
ducing  amax  as  far  as  possible.  To  reduce  the  cost  of  this 
preliminary  trapezoidal  integration,  a  coarse  integration  mesh 
(Aa  -  («max  *  «min)/50)  is  chosen.  The  initial  value  of 

o__v  comes  from  a  conservative  estimate  of  a _ „  which  ac- 

max  max 

companies  n(a)  .  The  integration  proceeds  step-by-step  from 
amin  the  relative  changes  in  all  quantities 

(0V  M  »  %  M  »  Pv  ,\j ( 0 i ) )  become  smaller  than  a  prc-spccificd 
number  6Q  (currently  0.001),  at  which  value  of  a  the  inte¬ 
gration  is  stopped  and  a  taken  as  the  new  ot  .  Iliis  pro- 

max  * 

ccdurc  proved  entirely  satisfactory  in  the  eases  in  which  it 
was  compared  with  Deirmendjian^l ,  although  the  scheme  tended 
to  choose  a  *s  somewhat  larger  than  his.  When  this  was 
the  case,  the  need  for  the  larger  a  *s  was  indicated  by  the 
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differences  between  our  results  and  his  tables;  when  we  used 
his  a|jjax,s  instead  the  differences  were  no  longer  found. 

We  also  observed  that  the  angular  interval  needed  to 
resolve  the  forward  peak  in  1*  ^(0)  is  inversely  propor¬ 

tional  to  amax  *  but  that  progressively  larger  angular  steps 
could  be  taken  away  from  the  peak.  The  sharpest  forward  peak 
will  result  from  particles  with  a  *  a  •  The  forward  dif¬ 
fraction  peak  is  proportional  to 


where  z  ■  asinO  ,  and  is  the  spherical  Bessel  function 

of  order  one.  The  quantity  is  shown  in  figure  6.3.  Since 
the  function  changes  rapidly  over  an  interval  Az  ■  c»A0  of 
unity,  we  choose  as  a  measure  of  the  angular  width  of  the 
sharpest  forward  peak. 


flu 


figure  6.5.  Graph  of  f(z). 
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(AO) 


o 


~ —  (degrees)  . 
max 


In  addition,  we  apply  the  restriction  that 


@.125°  <  (A0)q  <1°  . 


Ihc  tabulation  takes  place  for  enough 
cover  the  forward  peak.  The  remaining 
vals  of  a  pre-set  number  of  2(A0)  , 

8(A0)q  to  complete  the  set,  with  the 
AO  <  5°  . 


intervals  of 
angular  mesh 
then  4(A0)q  , 
restriction  th 


(A0)q  to 
has  intcr- 
then 
at 


Wlth  “max  and  thc  a,,Slllar  mesh  determined,  the  final 
Mic  calculation  is  performed  using  Romberg  integration . I69J 

Trapezoidal  sums  for  thc  integration  over  «  arc  calculated 

for  step  sizes  of  Aa  ,  Acx/2  ,  Aa/4  and  Aa/8  .  The  value  of 

Aa  is  determined  from  Na,  the  number  of  intervals,  which  is 
defined  as: 


I  25 

for  a_  < 

25 

N  «  J 

max 

a  1 

(25  ♦  0.55  («.„  -  25) 

for  a  > 

max 

25 

with  thc  constraint  that 


N  <  SO  . 
«  — 


Aa  is  then  taken  as 


Aa 


a 


max 


-  a 


mi  n 


K 


Mu 
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The  four  trapczoi da  1  sums  arc  then  used  in  the  Romberg  deferred- 
approach-  to-  the-  limi  t  scheme . 

The  above  scheme  has  proven  highly  satisfactory  in 
duplicating  Dei rmendj ian ' s  tables  without  the  need  for 
trial -and-error  determinations  of  the  Act's  .  The  adjustment 
of  a  few  simple  parameters  would  enable  the  scheme  to  work 
for  any  class  of  size  distributions. 

6,4.2  Truncation  of  the  Phase  function  in  the  Forward  Peak 

Hansen  has  shown  that  the  forward  peak  of  the 
phase  function  can  be  truncated  without  significantly 
changing  the  results  of  the  transport  calculation.  Thus, 
for 


0  <  0  <  0*  , 


where  6*  is  the  truncation  angle,  the  real  phase  function 

P  (z,Q)  is  replaced  by  a  truncated  phase  function  P*(z,G) 
v  t  . 

in  the  manner  of  figure  6.4.  The  logarithm  of  Pv  is  taken 

to  be  a  linear  function  of  0  on  [0,0*]  »  and  the  values  of 

P*  and  its  derivative  arc  matched  to  those  of  Pv  at  0  =  G*  . 

The  truncation  scattering  angle  G*  is  selected  ac¬ 
cording  to  the  following  rules.  If  Pv  for  0  =0  is  less 

than  a  prescribed  P  ,  P,  is  not  truncated  at  ail.  The 

*  max  v 

currently  used  value  Cor  P  is  5  .  We  also  require  that 

4*  4  Hill 

0*  £  20°  and  Pl(z,G)  <  P  (z,0)  .  Within  these  constraints, 
is  chosen  to  give  the  maximum  value  of  P  (z,0)  bounded 
above  by  x  . 

Truncating  the  phase  function  modifies  the  scattering 
terms  in  the  transport  equations.  Let  l:  be  the  fraction  of 


14  7 


Figure  6.4.  Example  of  truncated  phase  function. 


the  scattering  which  the  truncated  phase  function  fails  to 
take  into  account: 


*  fc/p  ' 


Pt)dfi' 


The  neglected  scattering  can  be  approximated  by  a  delta  func 
tion  in  the  forward  direction, 


- -  t J.yaiWJ  ft .  *-fUYXrKtii*  W!  -«*.  *».•-■« ! 


_  ^ 
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) 

| 


P  =  P  -  P*  +  Pt  , 

*  f  <s(3'  -  3)  +  p1 


The  approximate  phase  function  must  be  correctly  normalized, 
imposing  the  requirement 


77 •ft  «(S'  -  S)d8*  ♦  i_ Jf'itt’  -  1  , 

f?  *  1  •  F  ■  1  ■ 

f  *  4itF 

The  second  step  comes  from  the  definition  of  F  and  the  assump¬ 
tion  that  P  is  correctly  normalized  (to  unity).  In  this  ap¬ 
proximation  the  scattering  terms  of  the  radiative  transfer  equa¬ 
tion  become 

scat,  terms  =  /p(3,3' )  I  («'  )dn'  -  3  1(3) 

-f?/[47rF  +  Pt]l(^')dfi'  -  3  1(3) 

-  §¥ypt(3,3,)i(«,)d«'  -  Bd-F)i(S) 

=  irffa  I(3')dft'  -  B(l-F)I(fi) 

=  -  3'I(3)  , 


i 
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who  re 


e '  -  (i 


I'JB 


9 


1" 


Therefore 
P  and  3 
quantity 


,  the  scattering  terms  have  their  original  form,  hut 
lire  i  ep  laced  by  1*  *  and  p1  defined  above.  The 
P  can  be  calculated  as  follows: 


V  =  =  \  f  [I’COJ  -  P^OjJdo 

Jo 


i  r° 

=  J  I  [P(.0j  -  l,t:(0)Jd0 


since  ljt  =  P  lor  0  >  0*  .  Using  the  analytical  form  of 


log 


11) 


a  +  b  0 

s 


we  obtain 


P 


t 


e  '  0 
e 


s 


l<'a  ,  1)’  ---  b  f.n  11) 


ih  on 


P1  ( i.i  J  d  0 


jVv  V 


1  51) 


SMi  1  (»'.l 


Tlu*  integration  of  I’  over  tlu'  forward  peal,  mu  I  In- 
performed  numerically,  since  1*  is  only  tabulated  f « •  r  a  J  i  •• 
Crete  set  ol  (J  It  is  not  possible  t"  use  a  high-order 
integration  scheme  such  as  Romberg  or  (iaussian  due  to  the  uu 
equal  spacing*  of  the  0  points.  We  have  chosen  to  m.ihe  a 
piecewise  quadratic  fit  to  I’  ,  and  integrate  the  quadratic. 
If  I’  ,1’.,  and  1’t  are  three  phase  function  values  cor  res- 
ponding  to  Uo>0i  ,  and  t,  ,  then  the  quadratic  satisfying, 

Q(0j)  =  I’i  0=0, 1,2) 


is  given  by 


2 

Q(°j  =  ^  l\  1.(0)  , 

i  *0 


7 


j  =  u 


where  |  |  indicates  a  product  o  I’  terms  with  the  i**1  term 
deleted.  The  integral  is,  then 


I  hi 
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g( o) do  - 


d(e 


c 


°0>/» 


9 


d  “  (°o  4  °2)/2 


c 


■£ 


0. 


1 


^•4.5  Az i mutlia  1  Integration  of  Phase  Function 

The  phase  function,  !’v(z,ft,ft»)  ,  determines  the 

fraction  of  the  photons  within  a  volume  element  around  z 

traveling  in  the  direction  15'  that  will  be  scattered  into 

the  cone  dfl  of  directions  around  ft  .  p  depends  only 

on  the  angle  between  ft  and  ft*  ,  where  cosine  is  u 

*  * 

If  l’v  is  integrated  over  the  azimuth  angle  <j>  ,  an  aver¬ 
age  phase  function,  Fv(z,i.,u')  results  which  does  not  de¬ 
pend  on  $  or  $ '  ,  as  discussed  in  Ref.  1.  It  is  this 
azimuth-averaged  quantity  which  is  relevant  to  the  problem 
of  radiative  energy  transfer  in  a  plane  parallel  atmosphere: 


(  2  .U.U') 


1 

2  77 


/ 


2tt 


•*v(  2  .MjjJdG 


1_ 
2  u 


./Ti- 


(  I  I 


»  ) 


)  cosij  |ilv 
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p  and  p'  are  cosines  of  the  angles  0  and  0'  between 
the  positive  z-axis  and  ^  and  ,  respectively. 

If  P^(z,ps)  is  a  sum  of  Legendre  polynomials  in 
pg  ,  the  integral  can  be  done  analytically.  In  the  present 
case,  however,  P^  and  pg  are  given  only  for  a  finite  set 
of  points.  Hence,  we  must  numerically  integrate  the  tabular 

data. 

The  problem  is  to  generate  a  set  of  (P  ,<}k)  from 
the  given  set  of  (P  ,ps  )  for  0  <_  £  2tt  .  1The  trans¬ 
formation  is  J  J 

Ps  =  pp'  +  /(1-p2)  C 1  -  P  ' 2 )  COS  <j>  , 

*  Cj  +  C2  cos  <j>  , 


or 


<J>  =  cos 


(6.62) 


The  relationship  of  pg  and  $  is  shown  in  Figure  6.5. 

Note  that  y  s  ( <}> )  is  symmetrical  about  tt  .  This  means  we 
may  limit  the  numerical  quadrature  to  0  <_  4>  £  tt  .  The  set 

of  (P  ,$•)  is  generated  from  the  (P  ,p.  )  by  selecting 

1  •  j  j 

all  those  p.  which  give 


-1 


< 


C- 


< 


1 


and  computing  the  from  Lq,  (0.62).  There  may  be  no  pg 

which  yields  0  or  n  for  ,  yet  these  values  are  needed 
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-1  -J 


Figure  6.5.  The  scattering  angle  yg  vs  the 
azimuth  angle  <J> . 


to  complete  the  integration  at  the  end  points.  In  this  case, 
an  interpolation  is  done  for  as  a  function  of  yg  to 

yield  a  Pv(z,y*J  »  where  y*  is  the  value  required  to  give 

$  =  0  or  <J)  =  tt  . 

Given  the  complete  set  of  (P  ,  <Jk)  ,  a  quadrature 
method  may  be  applied.  Trapezoidal  integration  is  being  used 
currently.  Since  the  set  of  scattering  angles  has  been  chosen 
carefully  to  define  the  phase  function  well,  and  since  the 
rapid  change  in  the  phase  function  near  the  peak  is  eliminated 
by  truncation,  it  is  likely  that  the  trapezoidal  integration 
scheme  will  be  adequate. 


15  1 
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The  calculation  of  described  above  must  be  car¬ 

ried  out  for  each  non- redundant  pair  (y,y')  from  the  set  of 
Ui  .  Symmetry  relations  resulting  from  the  definition  of 
Pv  are  used  to  reduce  the  amount  of  computing, 

V^i’V  =  Fv(z’VPi}  (6*63) 

Fv(z,-yi,-p.)  =  Fv(z,p.,pj)  .  (6.64) 

By  combining  Eqs .  (6.63)  and  (6.64),  one  can  also  show 

=  Pv(z  , Vij  ,VJ±)  .  (6.65) 

From  Eq .  (6.63),  and  referring  to  Figure  6.6,  all  points  above 

the  positive  sloping  diagonal  may  be  derived  from  points  below 

the  diagonal.  From  Eq.  (6.65),  points  in  octant  6  map  to 

points  in  octant  1  and  points  in  octant  7  map  to  points  in 

octant  8.  Hence,  Pv(z,p,p')  need  be  calculated  only  for 

combinations  of  ( ij  ,  jj  ’ )  which  lie  in  octants  1  or  8.  In 

addition,  if  either  y  or  p'  is  ±1  ,  then  P  (z,ii  ) 

>  vv  ,Hs' 

has  no  4’  dependence.  The  integral  becomes 

1 

pv(z,F,u')  =  27  /  Pv(z,yy')d4>  =  P  (z,pp«) 

Jo  ' 

Pv(z,PU’)  is  found  by  setting 

Fs  =  FP ' 

and  interpolating  in  the  set  of  P(z,p  ) 

Sj 


155 


3SR-1034 


6.4.4  Renormalization  of  the  Phase  Function 

Grant  and  Hunt^^  stress  the  importance  of  using 
numerically  properly  normalized  phase  functions.  In  terms  of 
the  angular  quadrature  scheme,  this  implies 

m 

\  2  Cj  pjk  =  1  (k=l ,  . . . ,  m)  (6.66) 

'  j  =  l 


where 


156 


3SU- 1034 


P 


Jk 


Jk  =  1 v(z»pj  ’pk 
jk  3  Fv(z’VMk 


The  importance  of  Eq.  (6.06)  lies  in  the  fact  that  it  ensures 
flux  conservation. 

Eq .  (6.66)  would  be  satisfied  exactly  if  F^  were 
free  of  error.  Errors  are  introduced,  however,  from  the  fol¬ 
lowing  sources: 

(1)  the  truncation  of  the  Mie  series  summa¬ 
tion  in  MIESCT  (this  is  a  relatively 
inconsequential  source  of  error)  ; 

(2)  the  integration  over  sizes  (Section 
6.4.1)  ; 

(3)  the  approximate  integration  to  obtain 
the  area  under  the  forward  peak  of 

Pv  ,  when  truncating  (Section  6.4.2); 
and 

(4)  the  azimuthal  integration  (Section 
6.4.3) . 


Basically,  there  arc  two  ways  in 
can  be  applied  to  the  calculated  matrix 
additive  correction  can  be  made: 


which 

calc 


a  correction 
An 
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calc 

) . ,  =  p  •  i  ♦  e  . , 
jk  Hjk  jk 


or,  a  multiplicative  correction  can  he  applied 


■  “  *  V  ^klc 


In  either  case,  there  arc  unknowns  e.,  (the  number 

2  ^  J" 

is  reduced  from  m  by  the  symmetry  requirement  e.  .  =  e..) 

KJ  J  k 

and  there  arc  only  m  equations  (liqs.  (6.60))  to  determine 
them.  This  underspecification  can  be  handled  in  many  ways, 
no  one  of  which  seems  much  preferable  on  physical  grounds 
over  any  other.  Grant*  corrects  only  the  diagonal  elements, 
e..  *  ei  6  •  i-  >  and  so  *ias  a  determinate  set  of  equations  for 

J  K  J  J  K 

the  e.  .  He  points  out  that  the  phase  matrices  p...  arc 
J  J 

usually  strongly  diagonally  dominant  and  that  his  procedure 

thus  incorporates  the  entire  correction  where  it  makes  the 

smallest  relative  change. 

An  alternative  procedure  is  proposed.  First,  assume 


(1) 


and  solve 


m 


j  =  i 


=  i 


(k=  1 ,  .  .  .  ,  m) 


for  the  e  ^ ^  ;  then  assume 
J 


* 

Private  communication . 
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e 


jfc 


e 


(2) 

k 


and  solve 


1 

2 


m 


E 


c. 

J 


calc 

’jk 


(*  *  «*2)) 


1 


(k=l ,  .  . .  ,  m) 


for  the  e 


Finally,  to  ensure  symmetry  of 


form 


e 


•j1’  *  42) 

- 2 - 


which  is  to  be  applied  as  the  multiplicative  correction  to 
C  (1 1  c 

Pjk  •  No  decision  has  as  yet  been  made  between  the  latter 
procedure  and  that  of  Grant. 

6.5  CODIi  ORGANIZATION 

In  Figure  6.7  a  block  diagram  of  the  code  is  given. 
Those  sections  for  which  computer  coding  and  testing  arc 
complete  arc  enclosed  in  solid  boxes;  those  parts  which  arc 
not  yet  programmed  or  only  partially  programmed  arc  enclosed 
in  dashed  boxes.  A  brief  description  of  each  of  t lie  compo¬ 
nent  parts  is  given  below. 


INPUT  Provides  for  read-in  of  zoning  flags, 

atmospheric  specifications  (pressure, 
temperature,  mixing  ratios,  aerosol  con¬ 
tent,  etc,),  code  parameters,  option 
flags,  and  editing  flags. 

SliTUP  From  user  input,  generates  zone  struc¬ 

ture,  angular  groups,  and  frequency 
groups . 
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SOLCON 

Mil:  SCI 

Mil: 

TRNCAT 

phi i nt 

RAYLlil 

RLN'ORM 

I'LKAVG 

SPCTRM 

BUSS 

Rlil-MTX 

LOTRAN 


Irora  date,  time,  longitude,  and  latitude 
generates  solar  zenith  angle  and  kai th¬ 
rum  distance. 

Computes  Mie  scattering  and  extinction 
coefficients  and  phase  function  for 
a  single  spherical  particle. 

Integrates  the  Mie  quantities  from 
MlliSCI  over  a  size  distribution  of 
spherical  particles. 

Truncates  the  forward  peak  from  the 
Mie  phase  function,  ir  necessary,  and  re- 
rc* adjusts  the  Mie  quantities  appro¬ 
priately  . 

Integrates  the  Mie  phase  function  over 
r. :  1  mu  t  h  . 

Computes  the  Rayleigh  scattering  coef¬ 
ficient  and  phase  function. 

Re -normalizes  the  total  phase  function 
numerical ly . 

Computes  B/w  ,  the  average  of  the  Planck 
function  B^  over  the  frequency  group 
Av  . 

Computes  SAv  ,  the  average  of  the  solar 
constant  Sv  over  the  frequency  group 
Av  ,  correcting  for  the  actual  I-arth- 
sun  distance. 

Returns  a  surface  emissivity  e  (u)  for 
the  frequency  group  Av  accordYng  to 
user  specification  of  t he  surface  compo¬ 
sition. 

Returns  a  surface  reflection  matrix  r,. 
appropriate  to  the  frequency  group  Av  1 
and  t lie  user  speciiicd  surface  composi¬ 
tion  for  use  in  liq.  f(>.59)  of  t lie  Grant 
and  Hunt  algorithm. 

Computes  transmission  functions  T. 

Av 
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1  I TTRX 

TRNI’RT 

hi)  IT 


lit--;  the  transmission  function  T. 
by  a  sum  of  exponentials. 

Solves  a  monochromatic  radiative  trans¬ 
fer  problem  by  the  Grant  and  Hunt 
method  (Section  6.3). 

Prints  or  writes  to  tape,  at  the  user's 
option,  fluxes,  intensities,  and  heat¬ 
ing  rates. 
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7.  PARAMETERIZATION  CONSIDERATIONS  -  FUTURE  PLANS 

The  next  contract  period  should  see  the  completion  of 
the  family  of  mcso-scalc  codes  which  can  be  used  to  form  a 
detailed  parameterization  of  the  wave  drag  caused  by  flow  over 
mountain  ranges. 

The  test  problems  completed  to  date  have  given  some 
indication  of  the  effects  of  atmospheric  stability,  compressi¬ 
bility,  and  moisture  on  momentum  transfer.  In  particular,  the 
results  detailed  in  the  last  semiannual  report ^  indicated 
differences  in  wave  drag  edits  of  factors  of  greater  than  two 
due  to  atmospheric  stability  and  horizontal  wind  profiles 
alone.  When  compressible  effects  were  included  in  the  formu¬ 
lation  of  the  equations  and  the  results  of  a  test  problem 
compaied  with  those  from  the  Boussinesq  code,  differences  of 
the  order  of  25  percent  are  indicated.  Moisture  effects 
studies,  which  have  just  been  initiated,  indicate  small  dif¬ 
ferences  in  momentum  transfer  when  compared  with  the  Boussi¬ 
nesq  results. 

S3  intends  to  complete  two  additional  meso-scale  codes 
which  will  be  used  to  stud)’  the  effects  of  turbulence  and 
Corilis  forces  on  momentum  transfer.  With  the  family  of  HAIFA 
codes,  a  matrix  of  problems  will  then  be  completed.  The  re¬ 
sults  will  allow  a  designed  and  detailed  parameterization  to 
be  formed  of  momentum  flux  caused  by  orographi cal ly  induced 
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induced  flow.  Hopefully,  this  parameterization  will  then  be 
tested  using  one  of  the  global  circulation  models  presently 
available . 

The  next  several  months  of  work  on  the  atmospheric 
radiation  portion  of  the  contract  will  have  as  its  primary 
thrust  the  completion  of  the  radiation  code,  its  testing,  and 
comparisons  with  the  Mint z - Arakawa  radiation  package.  To 
this  end,  several  further  tasks  must  be  completed,  foremost 
among  these  is  the  exponent i a  1  -  sun  approximating  subroutine 
to  be  developed  in  collaboration  with  Prof.  Ivans.  Next  is 
the  completion  of  the  coding  of  the  Grant  and  Hunt  algorithm. 
In  this  connection,  a  basic  set  of  data  and  analytic  represen¬ 
tations  for  surface  cmissivities  and  reflectivities  must  be 
programmed.  The  assembled  code  must  then  be  debugged  and 
tested  against  any  known  solutions  (for  example,  the  ones  for 
pure  Rayleigh  scattering).  Comparisons  with  the  Mintz- 
Aral.awa  radiation  parkage  will  then  he  done  and  used  as  the 
basis  of  a  paper  to  be  delivered  at  the  AMS  Conference  on 
Atmospheric  Radiation  at  Port  Collins,  Colorado,  August  7-9, 
197,!.  further  studies  to  validate  the  primary  approximations 
in  the  code  are  also  planned,  but  may  Im  squeezed  out  by  lack 
of  time  between  now  and  the  conference. 
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APPENDIX  A 

THE  LEDESGUE  INTEGRAL 


The  Lcbesgue  integral  was  originally  devised  as  an 
extension  of  integration  to  functions  which  are  so  discon¬ 
tinuous  as  not  to  be  Kieiuann  integral) lc.  However,  t lie 
Lebcsgue  integral  concept  is  also  useful  in  the  integration 
of  continuous  functions  (such  as  line  spectra)  which  vary  so 
wildly  and  erratically  as  to  have  the  aspect  of  discontinuous 
1‘uuc  t  i  ons  . 

Viewing  integration  as  an  approximation  process,  the 
fundamental  partitioning  in  Kieiuann  integration  is  of  the 
a b s c i s s a ,  as  in  figure  A-l.  Figures  such  as  trapezoids  are 
then  constructed  using  the  partitioning  intervals  as  bases, 
and  the  sum  of  the  areas  of  all  tneso  figures  is  taken  as 
an  approximation  to  the  integral. 


y 


Figure  v -  1  .  Kieiuann  integration. 
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In  Lebesgue  integration,  it  is  the  o r d i n a t e  rather 
than  tlic  abscissa  which  is  partitioned.  Consider  the  shaded 
strip  between  y^  and  y  ^  +  ^  in  figure  A- 2.  The  x-intervals 
in  [ a ,1) J  for  which  f(x)  c  [>’i»>'i  +  1]  ;irc  oi'  lengths  A^ 

A A,,  A ^  respectively.  The  sum 


is  called  the  measure  of  the  set  {x  c  j  a  ,b  J  :  f  (x  J  f- ;ly :  ,y  ■  +  j  J  ) 
Multiplying  a  typical  value  fi  qyt  +  1 )  by  the  measure  a. 
and  summing  over  all  i  ,  gives  an  approximation 


(A.1J 


to  the  Lebesguc  integral. 
>' 


l'igure  A- 2.  Lebesgue  integration. 

A-  2 


3SR-J034 


The  formula  (A.J)  has  the  appearance  of  a  Rlcmann 
quadrature  rule  like  Simpson's  rule  or  Gaussian  quadrature. 

The  nulogy  is  valid,  of  course,  but  the  resemblance  is  only 
superficial.  All  Riemann  quadrature  rules  of  the  form  (A.l) 
are  based  on  polynomial  approximation  of  the  integrand  f(x)  ; 
no  such  approximation  is  made,  or  even  seems  possible,  in  the 
case  of  Lebesgue  integration.  Also,  the  a.  in  Riemann 
quadrature  rules  are  proportional  to  the  length  of  a  single 
x-partition,  whereas  the  measures  a^  represent  the  sum  of 
the  lengths  of  an  indefinite  number  of  disjoint  x-part i tions . 
Furthermore,  the  a.  in  Riemann  quadrature  rules  are  inde¬ 
pendent  of  the  function  being  integrated;  in  Lebesgue  quadra¬ 
ture  rules,  the)-  depend  on  the  function. 

The  advantages  of  Lebesgue  quadrature  are  obvious  in 
the  case  of  functions  which  vary  rapidly  over  (a,b)  .  Rie¬ 
mann  quadrature  would  demand  an  extremely  fine  x-partition, 
involving  perhaps  hundreds  or  thousands  of  points,  whereas 
Lebesgue  quadrature  could  furnish  equivalent  accuracy  with  a 
much  smaller  number  of  y-parti tions .  The  a.'s  arc,  of 

A 

course,  function-dependent  for  Lebesgue  quadrature,  but  within 
certain  large  classes  of  functions,  e.g.,  randomly- vary ing 
ones,  the  a •  ' s  would  he  essentially  func t i on - i ndopenden t . 
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APPENDIX  li 

THE  VARIATION  OF  BU(T)  WITH  x 

Consider  ;m  arbitrary  zone.  The  total  opacity  <  in 
this  zone  is  made  up  of  three  parts:  (1)  a  Rayleigh  scatter¬ 
ing  part,  (2)  a  Mie  scattering  and  absorption  part,  and  (3)  a 
line  absorption  part.  The  Rayleigh  part  varies  as  p  /T  ,  or 
substantially  as  p  ,  which  is  exponential  with  altitude. 

The  Mie  scattering  part  is  found  experimentally*  to  decrease 
exponentially  with  height  or  be  constant  with  height  (a  de¬ 
generate  exponential).  And  the  line  absorption  part,  accord¬ 
ing  to  the  approximations  employed  in  Sections  (> .  2 . 2  and  0.2.3, 
varies  as  pY  ,  which  is  again  exponential.  Hence,  a  plausible 
interpolation  of  k  across  the  zone  is 


where  z'  is  measured  from  the  zone  boundary  zq  (z=zo+z'l  . 
The  different  components  oi  <  have  different  scale  heights, 
so  that  the  11  of  lq.  (B.l)  represents  some  sort  of  mean 
va  1  tie . 

Computing  the  optical  depth  i '  (again  measured  1 rom 
the  zone  boundary  ;  -  i(zo)j,  and  solving  for  z'  , 

R.  I'enndorf,  (ieopliy  s  i  cs  Research  Paper  No.  2!>  (1'dbl),  Ceo- 
physics  Research  Directorate,  AI-CRt. . 
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(B.2J 


z 1  =  -II  An  1 


(i  -  o) 

o  ' 


it  is  assumed  in  the  current  version  of  the  code  that  k 
docs  not  vary  substantially  across  a  zone.  Therefore,  z'/H 
must  he  small  compared  to  unity.  hxpanding  liq .  (R.2), 


(■  -  s-»!(i'rf!) 


(B.3) 


To  an  excellent  degree  of  approximation,  the  tempera' 
ture  varies  linearly  across  a  zone: 


o 


T'  is  the  lapse  rate.  Since  T'  0°/kin  ,  T  200°  ,  and 
z'  <_  1  km  ,  it  is  valid  to  expand 


/  It v  V  hv  /.  T’  \ 

•x1‘Ut)  j  CX|’  i.iyi1  •  y;  •  ) 


li  v/  k  I 
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exp 


l'r  !’  11  •"('  ' 
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who  re 


6  = 


T' 

T 

o 


11 


The  question  now  is,  is  it  valid  to  expand  the  binomial?  That 
is,  can  we  expand 


exp 


a. 


hv/kT 

o 

e 


(» 


(B.4) 


who  re 


„  -  ,rT 

C  “  o  — 


O 


hv 

IT 


r 

o  To  Ko 


The  answer  is,  yes,  if  c  is  small  compared  to  unity.  Taking 

the  worst  possible  cases  for  the  various  quantities,  A  -  3p  , 

1 


T  =  200° K 


T'/T  =  8/200  -  0.0-1  km 

o 


leads  to 


e  •- 


t  ' 


K 


O 


Hut  from  hq .  ( B  .  3  J  ,  with  i '  <  1  kiii  ,  this  means 

c  <  1 

Hence,  for  mos t  cones  in  mos_t  I  requency  groups  the  expansion 
(B.4J  is  valid.  The  cases  for  which  it  is  not  valid  (small  \  , 
small  TJ  will  generally  be  associated  with  negligible  values 
of  11  (T.)  .  Using  the  expansion  (11.4)  in  the  l’lanck  function, 
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c  l’ 


A 


(1-ej  -  1 


a, 


Wk'T 


-  1 


hv/kT 


1  * 


liv/k' 


o 


which  is  1  in  on  r  in  t  because 

cut'  =  t  -  T 

o 


(B.5J 


Hence,  it  is  reasonable  to  approximate  By  as  piecewise 
linear  in  t  through  the  mesh. 

The  final  step  in  liq,  (B.5)  assumes  that 


rr:= 


If  hv/kT  <*  1  , 


B 


v 


*  A 


hi 

iiv 


which  is  linear  in  z'  ,  anJ  hence  approximately  linear  in  i 
by  liq  .  (  B  .  3 )  . 
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COMPRESSIBLE  CODE  LISTING 


c-i 


3SR-1034 


U  O  O  U 

a 

n 


ft 

3 

•  2 

• 

•  b» 

• 

ft 

Li 

•  > 

• 

Cl 

•  4J 

• 

o 

•  ll 

• 

ft 

o 

•  * 

• 

• 

ft 

L* 

ft 

•  VI 

ST 

cr. 

Ci 

•  u 

• 

o 

9 

9 

Li 

• 

2 

•  L! 

H 

ft 

•J 

o 

*: 

•  Cr* 

ft 

LI 

«f 

ft 

ft 

•t 

•  C» 

bi 

9 

0. 

L 

ft 

•£ 

•  L 

2  n 

H 

U 

o 

>• 

•  r 

•  V 

O  2 

L  L 

n 

ft. 

2 

Cl 

bi  o 

U.  O 

c 

CL 

9 

•  z 

•  L 

ft  c? 

9  Cl 

bi 

Cl 

• 

r*  •» 

t:  ft  cl 

ft 

•J 

Si 

»• 

• 

•  mi 

V*  ii 

-J  J 

9  X 

•  u 

• 

fc  ft 

2 

9  « 

ft  o  |,l 

a 

2 

•J 

•  a 

•  ft. 

b'  L, 

Cl 

C»  ft 

m f 

b* 

t 

L 

•  9 

• 

ft  b- 

H 

-i  ft 

Li  a  9 

2 

ft 

•  ft 

ft 

•X  ft  Im 

Li 

lift 

bl 

Li 

•  X 

• 

Z  Z  K  t 

X 

O  X  mi 

mi 

O  Li  *t 

n 

~i 

•  tc 

•  H  O  O  w  w 

•J 

ft  t> 

o 

>:  b-  fr 

b-i 

O 

v 

0) 

•  bl 

• 

bl  ft 

9 

f  ft 

ft 

bi 

ot 

H 

•  c 

•  a 

ft  b*  X  2 

Cl 

Li 

O  »•  u 

Cl 

tr  ft  C 

o 

r 

VI 

•  c 

• 

(DUO 

•J 

.J 

L.  o 

o  o  l 

• 

ft 

c 

ft 

•  L 

•  «* 

»  *  L.  ft!  ►; 

•i 

t 

ft  L. 

L 

b  L  1  • 

>* 

0 

X 

»• 

Li 

•  ft 

• 

ft:  cc  ft  ft 

O 

ft 

Vi  b  2 

•* 

b- 

r 

•c 

L 

ft.  ft 

r  • 

CC 

• 

•  > 

b '  b  Cl  C 

C» 

li<  c> 

Ci 

2 

2  ft  ft 

• 

V» 

•J 

1  • 

bi 

•  • 

• 

r  ti  ii  u 

a 

.1  Li 

O 

o 

*T  2 

a 

Li 

lb 

ft  K 

O 

a: 

•  u 

• 

fr  l. 

O 

ft 

C»  Ci  o 

u.  •  ■ 

fti  «y  «i  «| 

b 

Li 

J  L 

C» 

ft  Ci 

fr 

•  sc 

•  L. 

K  ft  H  h  bib 

#' 

a 

>•  i*  Z  U 

li  ft  ft 

b-  ft  St 

f » bi 

t/'  .i  1.1  Cl  ». 

2 

a 

•  I. 

• 

O  Ci 

l: 

!•■ 

'Ll  i » ;•  i,  L 

L 

c 

C> 

n  vi  bi 

b  ft 

C*  Cl 

»  ft  bl  VI 

o 

fr 

ss 

• 

•  O 

z  z 

bi  ft 

ft* 

o  cr 

L. 

Li 

V  X.  I  " 

V  9 

•t  ft 

bi  Cft  L»  J 

0* 

2 

ft 

•  L' 

• 

ft!  M  ft  ft 

ft  ft 

r: 

Ii  ►  L  ft 

ft 

f-  ft 

f '  O  1  • 

b-  ft 

J  C 

ft  JL'C  b 

bi  r 

• 

•  Cl. 

• 

> 

9 

;  >  n  ft*  u.  x 

•e 

b  • 

bi  U  O  It 

2  •( 

L 

2  2  o  >; 

ft 

U.  ft 

ft  Li 

9 

z  s* 

•  X 

•  if 

ft  Vi  2  Z 

>  ft 

c: 

2 

Ci  f.  Cl 

V 

r 

«) 

.  c 

rt  Cl 

2 

2  C*  o  ft  ft  fr  a 

L 

ft 

•4 

O  o 

• 

• 

J  J  ft  Mb.  J 

fr  2  bi  o; 

cr 

2  z 

ft  Li  ft  O 

Cl  U  ft 

9  ft 

ZHH 

b-  n 

ft  ft 

•  Li 

•  o 

..» 

0: 

L 

bi  9  O  X 

«c 

ft: 

►  rc<  ii. 

V  ft 

9 

n  ln 

Cl  fi  f ' 

Ii 

St  Ci  L* 

ft.  ft 

•  Cl 

• 

ML.i  J 

r  ci 

ci 

Jr  b 

bl  b!  o 

2  b 

ft  ft 

b  ►-  b  2  2  ft  ft 

l.  r*  b » 

w  b 

mi  fr 

s  «t 

•  c* 

•  ft! 

l.'  ti  J  J 

ci  y 

b 

L 

r  a  I,-  ft 

:« 

2 

/ 

i’  v. 

C>  Li 

►'MlCUUO 

2  Li 

ft  Li  mi  »* 

ft  IL  « J  ft 

Li  r. 

•  u 

• 

•.*  L' 

rr  m 

ft 

U  <  f  r 

*  • 

bi  9  9  bi  1 

C»  b 

r>  o 

r  j  j 

O  ft 

Ci  CL  1. 

CL  Lb 

ft 

»:  ft> 

•  ft 

L  L  C.  Ci  ft  b- 

r, 

l 

y.  L  9  2 

ft 

it 

i  b  :• 

1 

n  l< 

c.  t  t  ft  ft 

1  1 

b* 

2  9  >  C*  2  C  * 

L 

ft  * 

•  VI 

• 

O  C» 

bi 

CC  ft  bib* 

bi  ► 

>-  L.  b  b-*  r 

J»J 

a:  > 

I.- 

b* 

!)  b* 

O 

fi  ft 

•  >1 

•  ft 

L  L 

ft  r 

ft 

b 

rb-it 

• 

b- 

ft 

b  •  bi  ft 

•f 

;r  »t 

fr 

ft  fr  ft  9  9 

► 

b* 

Oft 

ft 

nr 

•  X 

• 

fi  ll  o  u 

2  9 

v>  x 

;>  ft  «•.  b.- 

UU)HHOO*» 

c:> 

•  St 

Li  Li  :>  :i  o  fi 

1* 

a:  cr 

tt 

»  a 

o  & 

•  ft 

•  z 

b,  L 

b  y 

lj 

1.  *-  1.'  V* 

ft 

c 

(i 

Ci 

t- » . 

b- 

b  Ci  c» 

Cl  X 

•f  Li 

’■*  ft. 

• 

• 

ft'  Cfi  b-  b' 

O  1  i  1,1  <1* 

b*  |.  2 

ft 

O 

o 

ft  ►  #• 

n 

r»  c. 

2  •*  b-  Li  L> 

b- 

I'd  2 

Oft 

b* 

b-  m 

•  X 

• 

r  >:  b  ft 

n  ft 

* 

Ci 

ft  b  »•  n 

ft 

*4 

*»  -j  c*  n  o- 

9  i*  b 

b 

b 

ft 

cr  a:  «€ 

o-  4  »  ft  1 

•  H 

• 

r»  9  b-  b« 

9  *i 

ft 

;i 

•f  V'  l  j  2 

r» 

j: 

i. 

I.i  f '  Cl  ti  J  fr  (r 

L 

fr  o  n  *c  •» 

O  C»  b-  I  i 

O  bl  L* 

9 

«t  Cl 

• 

•  L 

x ft  l* 

C>  1 

► 

i. 

ft  L  ft  Ci 

»* 

o 

v« 

»• 

b  bi  n  ii. 

C  Ci 

L 

ft.  c,  C  ft  ft 

•• 

»• 

»*.ft  1> 

L  ffe 

© 

fr 

ft 

•  V- 

• 

•  •  •  • 

•  • 

• 

till 

•  •  • 

b  • 

• 

• 

ft  1,1 

•  ft 

•  Li 

•  •  •  • 

•  • 

• 

>  •  •  • 

•  #  • 

b  # 

b- 

b  i  i  •  i 

l’ 

• 

b  •  • 

•  ft 

ft 

• 

M  X 

•  M 

• 

•  •  •  b 

•  • 

.i 

*t  •  •  •  • 

•  •  # 

•  • 

i- 

zi  b-  o  ft. 

ft 

b 

ton 

ft  U, 

fr 

• 

bi 

•  if 

•  fi 

•  •  t  • 

1.  ft 

c 

3  •  •  • 

b  b  • 

b  •  b 

1  ♦ 

bit*  C.'  >.  >: 

a 

c> 

Cl  L  o 

•  •ft 

•f  i 

X  ft 

•  9 

i  i  i  • 

1  II 

ft 

<  •  •  • 

•  ft  b « 

r>  > 

•  *. 

ft  j  j  y  :i 

i.* 

lr 

b*  a  -j 

I*  41  b 

c 

l  ft 

•  • 

• 

ft  ft  ft  ft  bi  *: 

ft 

c* 

ft  •  f  bl 

t»  ft  o 

b-  i. 

ft 

r.  L  i  ft  n 

b* 

C»  i>  Li  4 

ft  Li 

ft  n 

b*  3 

•  mi 

X  2  C?  Ci 

ft  b 

c« 

b* 

C>  L.  b-  ft 

9 

J  n  f»  ci  b- 

9  ft 

b  i 

ft  ft  ft  bib* 

ft 

;• 

<L  ft  2 

r>b 

Mft 

x  o 

•  mJ 

• 

X  IT. 

•  X 

• 

• 

»  • 

•  O  Oil  O  Cj  UoUU  C»  Cl  4  •  ()  O  II  u  (I  c  (•  (>U  o  I)  I*  li  (  <  IH*  ll  (MM>  U  C  I)  If  |i  <’  II  I*  U  U  U  U  |l  <' 


w 

ft 

bl  MNh 
fr 

fr 

r 

o 

o 


#  Wtf  ©•  4  m  m  ^  <t  v  r  ii  f  i  ►»  .#  w  w  ft  r  ft  f '  < 


ifin.»wirMi|  r*Nt  ftjftuftrftn i«tiM 

iHiMih  j  +  *  *  4  #  #  #  4  i>*  *1  r>  m 


3SR-1034 


n  i«u  •- 
t  ••  *1  I  1*0* 
H  .12  *  * 
I  1  <  .If 
4  *  to**  ft  •  J 
If  4*  *4  4  #*!• 
4  • 

»*  J.I 

•  «!** 
#.  «s  •-  4  *  to 
Ml  4  «  C 


» 

•  • 

4..  • 

»  • 
fe  • 

K  • 

•  1.1  • 

>  • 
to  • 

41  • 

• 

to  • 

II  • 

• 

•  • 

»•  • 

II  • 

4  • 

4  • 

M 

l»  It  • 

41  •  • 

M  • 

I*  4 

II  •• 

»•  •  • 

If  4*  • 

•*  •  • 

:»  • 

r*  •  * 

(t  •  44 

4*  ••  U 

r  iii  ^ 

It  •  _ 

I*  II  *i 

•  ►  •  I 

♦  ill  #H 

♦  t  4  III 

4  4  .4  • 

♦  *  I  to  #4  • 

4  to  4*4  #  • 

4  1  4-  til* 

4  »IC  It! 
•If*  4*1  4*  • 

•0  +•  *  Jfc  »  • 

•n.4  to  l 

•*  •  •«  4'  <  >  * 

•  • 

r  • 


«« 4*4114  U«««l 


y«M* 


«*yii  «* 


►  jh 

5  • 

i  k  •• 
u»  ti 

s.\.c 

£  04  44 
if  ?4 
4  *•  ^ 

04  II 

*  . 

•»  y  «4 

•  4/  4"  «| 

•4  4*  to  M  * 
4  4j  4  f  #■» 
>  »  »  14 

4  to  5*  to  t  • 
4*  to*  4*4.  •  a  4 
!•  i.r 
*44  4*  *4  4*  *4 
4*  4*  »  *> 

•  4.  •• 

1*4*  I  I  I 

•  ♦  •  !•«  f* 


Ml 

*  ii 
4  4*  4a 
I-  •> 
I  1  mt  «• 


S  ? 

5  41  " 

|w  at  t*  I* 

U  •*  4* 

4  i  W  tl 

%  •  to. 

<M  «#  O  S* 

to 

•I  .1  » 

40  «*  4ft  •! 

4  to  4*  *i 

4  4  4  4 

•  • 


i«4H»  «*y 


yUMI"  lllllj  ««UI*  Ml 


s  **  •<! *  44  4*  4*  ■#  m  *■  f  MliH*  ||  M»  #-  *>  *-  4  00  0 

*MWirl*il4l4B«0*«4.-«  4f*f'4»M'*4  M  f*4  4<'  f  4  to  t  to  W;  *■ 


•*■**<*  «  a***#*1 *»#-,£.«-  «  #«*  *•>**  >4  *?{-■«* 

00I000  0IM0*  00  00 


3SR-1034 


O 

V 

H 

* 

ft 
Cl  M 

*s 

tiki 

it  ft 
ft  ► 

*  r 

» t 


i  v 

■  •» 


% 


«c 

O  4 

ft. 

*  ft 

ft  II 

«? 

1.0 


n 

ft 

r 

•« 


o 

sf 

kl 

«> 

O 

ft 

V 

m 


w 


«S 

► 

II 


99 

▼ 


s? 


**H  3h' 

9P  c  r 

,10  3 


r 

O' 

o 

kl 

if 

9: 

O 


o 

s 

£8 

88 

p>  • 

n*- 

•  Vi 

O  • 


ft  * 

e-*r 

•  u 
► 

M  l> 

4 

9. 

9*  • 
t> 
V 
t  m 
n  • 
y 
« 

•  N 

4 

r 

X  II 
4  ft 

I  O 


8 

il 


ft 

u 

k 

a 


k 

• 

ft  •* 

It 

•• 

9  * 

k  Cl 

fc< 

®  X  mi  mi 

9 

9> 

ki  n 

ft  ft 

9  ft* 

9 

*  H  J3 

r 

9  ft 

ft; 

c» 

S  ft 

A 

ft 

k*  *1 

k«9- 

#•  MM 

ID  ll 

• 

r 

9- 

kit 

%•  •) 

i*  i  • 

U9.  kl* 

# 

•  *  9 

19 

k* 

» 

l» 

y  a 

•1  k 

•I9« 

U  U 

9 

9- 

a 

M 

trt 

• 

M 

M|» 

• 

W  9 
k  « 

•9  * 

b  ft. 

l»  *  » 
9“  ft  »*k 

ft- 

9* 

mi 

('ll 

»! 

19 

r 

a*. 

9 

ft  ft 

9  fc 

9  Ii 

*•9  f  • 

9 

t  ft 

ft 

ii 

a  c* 

•  k 

•*  9 

•» 

l»P  •!  S' 

9* 

•  •  99 

•9 

4» 

t* 

s 

•J 

• 

•■  y 

f  mi 

•Jfi  9-  #* 

9  9* 

9-» 

9= 

c 

Ii 

•»  9  • 

1  mi 

!»•*•»** 

ft 

90  90 

ft 

ft 

c 

•f 

1* 

J  1 

S'  •< 

•  k  k  k 

9* 

1  k 

ft 

C 

ft 

• 

• 

1ft 

• 

• 

r  k 

«  II 

•  pm  9*9* 

• 

• 

• 

9«9t 

• 

• 

1- 

ki 

ft 


•i 

4 

l3 


J  oo 

33H 

r*  ft  n 

•  a 1,1 

►  C)  J 
It  rl 
*•9  4 

•  Jl* 

•»  4  . 

1»  • 

*• 

•  ^1* 

i»i*  r 

►  irfci-  . 

M4*l«  J 
»  *h  IX 
k*  I 

•  •• 

•Ik.*- 
••  ;*n  n  k 
I  xpjr 
t'  I*  9  •! 
nwMwon 
*»  *  0.  k  •  *  y 

•  Mil  Ml'k 

*«. 

t* 

•  • 


It  UClIHi  UMUI*U 


nil 


lilt  til 


HUM  1*1* 


€-4 


iht 


3SR-1034 


« 

ft 

V 

• 

m 

V- 

ti 

mi 

u 

ft 

l> 

H 

M  Jt 

z 

• 

ft  u 

ft 

ft 

ft 

•  43 

• 

ft 

c> 

*J  & 

ft 

• 

• 

U  J 

Cl 

ft 

« 

ft  • 

H 

ft 

Cl 

U  ft 

• 

1* 

• 

M  ft 

• 

►* 

• 

•  r  n 

• 

• 

CI 

«UP 

• 

ft 

• 

N  »H 

• 

• 

r 

•  as  ci 

Jr 

ft 

•  M 

M  M  ft 

■0 

• 

►  ft 

II  t  «• 

• 

► 

M  Jb  « 

HUP* 

ft 

Ci 

•  ft  e  • 

•  •  • 

ft 

• 

ft  c?  • 

•A  ft  ft  * 

Cl 

ft 

i-  r  i» «» 

H  H  l.»  ft 

• 

n 

«>  J4  • 

M  M*  • 

ft 

•» 

ft  *11'** 

r  ft  s*  ft 

Cl 

ft' 

•  ft  ft  •  • 

ft  H  M  • 

• 

9 

•Aft  1(1 

• 

*■ 

N  r  cr  aft  f « 

•-  ft  ft.  ft 

• 

X  € 

Ml  L'  •  • 

H  • 

• 

c»  :• 

ca  tl*  K  J 

n  •  jr. 

• 

^  r  r  •  -  * ’  "  7  * 

a  ►  M  M  «.  O  • 

.1  m  r^rKui.ft  ir 

K  7  K  l'  b  iu'hf  rit 

*  O  •  *.*•••*•!  +■ 

e  M*«  U  •-  *  ft  10  .J  tt 

;»  nv.  if  k  i  k  •*  u  f  b 

»  ••»  . 

r 

b. cl.  ••  n  j  •• 

3  ;»  I  «•  «  «  4  f  «  <■  •» 

^mh«iiusii*v 

•  • 


r 

Ci 

9 


a- 

• 

• 

>  mm 

ft 

*•  c% 

ft 

ft  • 

ft 

i*  t: 

•: 

ft  .a 

ft 

1.  1 

41 

li.  o 

to 

i)  Z 

• 

to  hi 

I? 

Ci 


II  |I  IIUII 


U' 

z 

H 

I* 

W 

n 


ft 

n 

t 

». 

ft 

o 

CJ 


fi  m  ♦•»»#  in  •f  •*  I*  J; 

L  M  ••  II  ••  MM  MMHM 


ft 

ft 

«» 


c-s 


C'5«®°£SS»I3LE.5EriN 


3SR-1034 


r 

a  h 
m  n 
a;  • 

•ftl  CJ 
►.  M 
M 

lr.  © 


U-  I*  ►  X 

V*  7!  n  n 

©  ►»  o  cr  *• 

w  •  §  ti 

Ll  H  K  ►  ►  X  »» 


a  cj 

(tl  N 

► 

1/1  ( 


► 

^  C1 
H  Cf 


I 

.  Vi 

a 

L>  « 

r.  ■< 

H 

h  • 


V. 

o 

H 

l/> 

a 

u. 

u. 

H 

© 

0 

o 

I* 

ft. 

Ui 
©  H 
•  W 
H 

I  It) 
M  »: 

M 
>  ► 

**  r. 
on  a 
.  1  * 
r  o  * 
s  •  *f 


-  ~  h  »: 

L.  ~  > 


O 


I  ,«  f  j  f.'  K 


Ll  X  ()  •  It) 

on  *  u: 

CJ  *  ft  u  •- 
X  •  S 

ft  c  L  L:  n  » 

k  li  m  t»  »-i  r  «j 


► 


or 

o 

it. 


H  U 
> 


1  > 

N  S 

ft' 

ft 

3 

r> 

c. 

C  o  t) 

3 

3 

M 

i 

1 

i 

i  i 

i 

h 

•  4  ^4 

r ' 

• 

• 

v- 

H 

iii  a  n 

It 

o 

li 

I) 

• 

«f 

• 

.J 

X 

>* 

X 

>•  >■ 

X 

V- 

7  X 

* 

M  0) 

Cl 

a 

li 

n 

i>  r i 

Ci 

5' 

7. 

>2 

7.  7. 

it’ 

►  •  *i 

r. 

1  •  • 

• 

o 

3  M  r4 

►4 

M  t? 

ST 

it  >i 

II 

M 

I* 

r 

II 

1' 

II 

II  II 

n 

V) 

» • 

on 

♦  t 

cr 

-1  ll  l> 

%* 

M 

v>  w> 

>:  7. 

o 

» 

13 

*•< 

♦  4 

*t 

ft;  Cj  M 

y 

•?  tf 

o 

li  © 

II 

CJ 

o  x  ► 

>: 

r. 

X 

>■ 

? 

J 

s.  a: 

a: 

o 

w 

c 

»- 

3. 

a. 

a. 

a 

» » 

11  II 

•  t 

*- 

li 

r.  r. 

3 

3 

c» 

cr 

3 

:j 

I  t  V* 

n 

J* 

It 

i* 

X 

>- 

>  >• 

y 

» 

t 

► 

V 

1 

Vi 

►  4  CL  ft 

Ci 

© 

CL 

a 

ll 

o 

»-• 

• 

►  4 

• 

y. 

y. 

y.  y. 

X 

r  > 

H 

r 

n 

er  *" 

•  <  l 


y  i 
3  < 
y  c. 


i  t  n 
o  o 


IS)  ► 

I  t  • 

n  u 

it  >• 
t  j 

>  r. 
r  • 


5*  t* 
r.  m 
it  x 
>*  n 
o  r. 
>  *■ 


i: 

3 

n  ij 
•  it 
>  u 
n  l 
n  it 


In 

l/i 

U» 

CL 

CL 

X 


r.)  HN  K,  »  III  u>  K  O'  O'  ft  »4I!>I  h  tf  IO  10  K  M  ft  r  Hf)  W  *  Iflio  K  tr  ft  r  H  fl  M 
HHHHXHMHit  *4  N  f4  N  ft  CS  ft)  M  ftt  M  Cl  M  M"  h  MM  H  Kl  H  h  iM  tf  4 


ft 

0 

At 


C-6 


5'JKI  m  rTpOT*var*TP3K*V3X 
I«-Cn'j>*CICI  ,LT.  l.C£-20l  N0DI«T=1C 


3SR-1034 


i 

i 


s 


s 

l 


X 

t} 

N  ~ 


n 

to  • 

• 

.J  • 

to 

9  • 

tfl  •  • 

O 

U-  •  K 

X 

ft  •  H  ~ 

ft 

•  H 

• 

U  •  •  N 

m 

C  •  >  t/l 

•> 

•  T  ft 

• 

to  ft  H  M 

to 

fr  t  •  • 

•t  H  X  X 

K 

• 

1*  frj  H  H  M 

M 

9  H  •  •  • 

n 

b»  •  D  H  O 

• 

T  X  ft  in  r« 

r  HHJtbi 

*0 

H  •  •  Mft. 

X 

*  S  X  •  H 

0»  » 

x  •  to  X  • 

C-‘  M 

t'  i r.  h  r  v  x 

n 

f  1  •  •  H  H 

M  • 

•  fsj  ►  *  •  • 

M  M 

X  H  I/I  O  »• 

•  •« 

•  Lift  ft  U3 

to  > 

0*  ft  ft  ft  • 

*0  » 

•  m  to  ft  to 

W  M 

»  v 
9 


•  •  ... 

•  X  •  1.1 

•  u  to  x  a 

•  •  O'  to 

•  U  X  •  • 

•  r.  1 1  ft  x 

i  h  •  «.f< 

•  to  b.  •• 


X 

•  •> 

•  ft  ft  to 

to 

►< 

rj  • 

ft  to  to 

• 

to 

to 

ft  •  •  « 

X 

• 

4  M 

H  X  X  ft 

to 

j 

to  ft 

m  U  IT  U 

• 

o 

f  1  to  •  • 

to 

to 

*•  M  9 

H  •  1  X 

to 

(> 

•  t 

•  9  ft  «*> 

• 

H 

H  H  "  *)  H 

x  to  to  • 

X 

ft  0.  ti  to 

to  •  •  to 

to 

HI'  - - 

(i  i;  M  to  x  wi  ••  to 
IUHO  ft  5!  •  t*  V* 

i«  ci  •  •  u>  n  •  to 

b1  to  Cj  H  M  ^  * 

35  ft  M  II  It  to  1  L 

to  •  “>  >  4  to  •  •  ft  t* 

to  l.l  IQ 
9  to  **  t.  to 


7b  II  «.J  f]  Ob’  H  M  I 
J  I  I.'  H  r  1‘  H  trl  |i  I  i- 
•4.1  hi!  **  H  »* 

kkh  un«  to  *-> 

M  U  .J  C-  r ►  ♦  >  Jb- 


w  x  - 

N  It  to  M  to 
»l»H 

bH|«  • 
jTHH 
«,<  U  H  a  N 
-  10  H  ♦  o  7  -  -  O  - 

Ct  w  *•  i?  it  b  ti  >  h 


r  •  -  «t  x  «  x  •* 
►-  :>  j  a)  >  c  >: 


.  »• «'. 
O  L  •  «i 
i'b  I* 


n 

b) 

• 

Ul 

to 

A 


u  n 
©  r 


o 

n 

ti 


fr  O 
O  3! 

It  It! 

1 « 

1 1 


VI 

uHNMf^Ab'e'ACIHNKuMMChAA 
ft  rl  H  H  H  M  M  H  H  »<  M 

«L 

ft: 

o 

u 


to 

u 

tk 

M 


C-7 


3SR-1034 


> 

ft 

• 

a 

SS 

A* 

# 

ft 

M 

O 

h 

x 

► 

«c  a  to 

m 

ft 

a 

a 

a  ft  ft 

a 

a 

► 

X 

X 

9X9 

•k 

a 

* 

ft 

X  t  J 

H 

H 

a> 

• 

t 

•r  n  w 

et  t 

a  a 

a- « 

«t 

X 

a  a  a 

i»i  X 

t<  x 

tl  ft 

X 

tl 

f'  <f  ft 

►  V 

X  X 

•  a 

ft 

a 

x  a»  c 

9  %• 

«c  •• 

ft  a 

mi 

A. 

ana 

J  • 

mi  • 

»4  ft 

ft 

X 

a  •  *4 

• 

a 

a  x 

a 

11 

fiao  U. 

r  x 

st  a* 

a-  ti 

9 

o 

|.  X  a 

t>  • 

It  a 

t*  i:  m 

x 

X 

a  ft  4 

X  X 

a*  m 

•i  ft  i 

ft 

ti 

a 

a  a  »i 

iA  X 

v  x 

a*  a  •» 

a 

X 

•I  Ai  X 

ft  « 

«• 

ft  m 

L.  ft  ft 

X 

X 

X 

H 

!•’  9 

► 

li  a 

M 

r.t  m  x 

X 

•« 

X 

>  !• 

a 

>  a-  x  c  it 

•i- 

X 

A. 

a  •  a 

X  X 

••  2 

a 

A  a-  a 

a 

t> 

A. 

1-  HA 

H  • 

*  a*  a 

ft 

ft 

fi 

a 

X 

X  X  W 

• 

a 

a 

a- 

►  •*  « 

a- 

a* 

X 

a 

a  Cl  9 

:•  ►  _j  li  ►  u  o-  •  ►- 

ft 

m 

« 

X 

>  •  »J 

c>  • 

«f 

£»  a 

at  a«  a  ui 

a 

H 

•: 

O 

•IX  U 

j  n 

i 

O  M 

• 

uft  t 

H 

ft 

N 

o 

-1 

a  X  • 

i.i  x 

m  x 

ft  o  x  a 

x 

ft 

X 

c. 

A  a'  A# 

at  *x  e 

ft  1' 

c. 

J  a  a 

ft 

ft 

a 

u 

X 

X>l‘ 

a 

K 

a 

*  t-  a*  a* 

a 

a 

ti 

a 

a 

rz« 

|,t  M  i.i  ft.i  at  it* 

>ar 

x 

M 

S> 

X 

X 

•  i*  j 

a*  a 

a  a* 

15  ft 

•: 

M 

It 

ft 

o 

X  X  t 

3  •< 

D  «n 

3  a  | 

a 

ft 

x 

r 

mi 

ft  a'  a 

a>  it 

a-  it 

a-  a. 

ft 

X 

a 

x 

tm 

a  at  a* 

«a  •* 

ft 

ji:K 

a* 

> 

a 

a 

H 

»•  X  O 

ft  x 

tL  X 

«  •*  a 

ti 

9 

X 

a 

ft  U  9 

t;  mm 

ti  •* 

a«  ft  a 

tu 

a 

C 

ft 

»• 

a  A.  mi 

a  x 

b  a* 

a-  A 

a 

M 

ta 

ft 

9-  •  ll 

S: 

t. 

a*  x 

9 

ft 

U 

X 

«a 

n  9  a 

!*• 

t* 

T  •» 

a- 

M 

a 

V 

ft 

a  X  A. 

a* 

x 

m  a 

Cl 

ft 

M 

a* 

A. 

o 

t:  Cl  ft 

a  f . 

3 

a 

a 

•  « 

o 

• 

kl  X 

w  ft 

ti  ft* 

at 

ft 

:a 

1) 

a 

a- 

ft  '• 

J 

MX*? 

a*  ft 

fc-'  a* 

*  x 

a 

X 

a 

•  •  a-  ci 

t  A 

• 

a*  x  x 

a*  J 

a* 

a-  ti 

a. 

*• 

ft 

N  ft  a 

•  •  a 

•  •«.  r 

t  ft 

t  r« 

b  K 

a 

a 

99 

ft  ti  t*  ii 

a 

j!  x  •  a 

1.1  • 

t;  • 

tl  • 

N 

m 

ft 

a* 

ti  A  tl 

n  t  r  n  a 

r.*  # 

f.  ft 

i'  # 

X 

#• 

t» 

mJ  V 

cr 

t  i 

a 

•  -at  v 

t* 

x 

x 

n 

S? 

t* 

ti  a 

a 

l.f  a  »  t.  x 

a-  • 

O  ft 

c>  r 

o  n 

a 

ti 

tr 

V 

ft 

• 

ft  • 

3  • 

a*  * 

a*  a 

a*  • 

h  n  if 

a  a* 

a-  a*  v 

a* 

x  at  a 

k  it  a 

n 

r* 

t« 

► 

I*  ft 

t 

t*  ft 

l*  ti 

a 

ft  SX  J 

V*  ;•  • 

a* 

ft  x 

a<  a- 

a 

«• 

A 

• 

<i 

ft  t 

i  •  t  •  y  pi  ?,* 

ta  M  A 

t  ft 

a  n 

A  C 

r 

t 

•It 

ft 

•It  •  V* 

LI  • 

t  a 

ti  a 

• 

H  A 

ti 

j  •  c«  a 

i  t 

•  • 

to  »»  !•  • 

k<  »*•  • 

a-  fi 

a  a* 

a-  ti 

i 

a. 

a 

a.  a 

• 

ft  a  ft 

:» 

1  • 

x  s  a 

i.»  x 

t.  t 

i  x 

a 

j* 

i 

a* 

a 

•;  a  a* 

t  •  • 

•  ft 

t  cv 

9  C* 

ft  j* 

a* 

f 

ft 

»'  1  9  t* 

t*  i.i  9 

«  • 

ft  a 

•i  a 

•» 

a* 

•  4 

.1  •» 

a* 

a 

to 

t*  l  *i 

dpi 

a  »• 

Co.  ft 

A  M 

i* 

a  • 

•i 

a- 

» 

9 

*  a-  t 

•  a  t  ti 

on* 

x  a* 

•a  » 

••  »• 

a 

* 

j* 

a* 

•  a 

:» *i « 

*  »  ,!.»♦* 

9  J  I 

Cl  •« 

a  t» 

C  tl 

t 

a. 

a 

*  *  *  *  ft  f 

t 

f  ft 

to  1  *  9  A 

a  u  a 

• 

a 

a 

a 

a 

J 

;a  a-  a 

* 

a. 

ar  # 

»'  > 

L  /  a 

a 

x 

a 

a 

*||I|5 

i  *•». 

• 

I*  •* 

•  •  •« 

an  m  • 

ft 

K 

K 

t 

t»  f* tr 

*a< 

a 

a 

f 

a  *» 

a 

at 

•  y  , 

a  at  A« 

A* 

Co 

tl 

»• 

f» 

tl 

.1 

n 

tr 

•K 

ft  ft 

tin 

At 

•4  H 

A- 

AI 

Ci 

•t  r 

r.* 

X 

»•► 

9 

r 

•• 

o 

9  f* 

t> 

to 

>  ft 

ti 

tr 

mi  m 

r* 

d  r 

ti 

A. 

a  x  a 

n  tv 

tm 

• 

a  •  a 

»•  l 

9 

t» 

a  9  a 

i*  t 

AJ 

If 

a  t  a 

it  i 

a 

a  %•  a 

u  a 

•1 

aa 

a  i  •  a 

K 

a  ct  a 

o 

t* 

a 

x  a  *• 

li 

X 

a  m  a 

x  ••  c  a 

X  o 

mi 

a  A  a 

t*  j  tr  tf 

p  ft 

ft 

a 

a  r-  a 

U  lllb 

9  • 

A 

a  x  a 

1  ft-  t>  A 

a-  f 

X 

A  »•  t,  fc  t“  •- 

ft  • 

t» 

A 

9  •*  i*ft  •  Cl 

a  9 

ft 

X 

Mi*  •  •*  »• 

V  «* 

\ 

t 

A  lo  •!  ft  •* 

••  tj 

a 

t 

5  9  10 

t*  1 

Co 

«l  til*  »>  1*  X 

ai  a* 

ft 

«a 

*4  «  -  9  •• 

t  x 

»• 

X 

9  9  t  t  f  X 

A  » 

a 

•  a 

t  a  !•«.  5t  i 

a 

•  • 

a 

*■  V* 

i'  n  • 

i  o  Fi 
H  V  I  I 
t*  •  I  • 

•  m  •  m 

*• 

-  t  -i 
«»*«•«» 


F 

M 

i!»  ui'unuuuouou 


(•  ti  u  ti  t  l»M  «'«»«•  w  t» 


•  J €*•*  u 


I 


'  ^  k»wr»*  •  •  9  A*  tA  A*  4?*  •  f  *  •*  A*  X  jP  *  t*X  9s  t  #"»  aAA#(N 


t;*a 


3SR-1034 


► 

• 

m 

•• 

« 

ft 

ti 

► 

X 

.  • 

*1 

• 

o 

* 

H 

* 

•i 

ft* 

ft 

*  — 

•t 

tl 

r 

•  ft- 

1 

ft» 

« 

•J  • 

ft* 

1 

ft 

• 

«  « 

ft. 

■A 

•ft 

•  ft 

M 

ft 

► 

m 

•*  ti 

tl 

M 

ft 

ft 

•  Ci 

• 

«.f 

►* 

ti 

M  1 

m 

• 

t*  *• 

ft 

ft 

« 

• 

•  A 

«• 

ft 

t* 

0 

m  m 

• 

IB 

o 

ft 

ft*  t* 

• 

r 

«  • 

m 

c 

a 

•  m 

ft 

z 

« 

ft 

►*  ft 

•i 

ti 

to* 

X  Cl 

«• 

r 

r 

•  • 

« 

• 

ftl 

• 

•  ft* 

ft- 

ft 

ft 

• 

»•  r 

K 

« 

3 

• 

•  • 

• 

• 

► 

Cl 

• 

• 

ft 

• 

ft 

ft 

t« 

•> 

•  ► 

• 

• 

O 

ft 

ft 

►  ft*  • 

ft. 

•• 

• 

e 

r* 

fti.vi 

ft 

ft 

r. 

J? 

• 

•  «  ►« 

■f 

r 

it 

#»  ft  u- 
fi  ft*  a 

Cl  t  »u 
9  ».  •«  19 
tl  •  •* 

►  JK> 

N  ti  y 


••ft 

to* 

ft  • 

V 

*  ft 

t  *  r  * 

ft  ** 

ft 

•9 

• 

ft  • 

to* 

to* 

n 

IV  •  * 

ft  *1 

•  fJ  1 

nit.* 

•to 

•  •-ft  t> 

•  ft 

•  •♦  A 

ft 

ft  t 

•  *  •* 

n 

MK'to 

ft  ft  ft 

ft 

ft  ft  ft  N 

••  ft 

Aft  fto 

t  t  • 

V  4.  •* 

r 

•*  »|  tot 

u  ft  ft 

« 

11  1 

m  *r 

M|.  | 

* 

o  ti  «■> 

/  9  1 

■a 

•«  • 

«•  i.i 

> 

•  A  Tl  ll 

«  It 

e  *i  c. 

o 

•  i 

t*  t-  • 

« 

X  O  9 

•»  mi  ft 

r 

t  1  A  •• 

ft  a 

»• 

fl  ft  to¬ 

f.  to-  H 

r 

•  ft  • 

O  M 

to* 

N  to  1 

*  •» 

a*  1  •« 

ft 

ft  r*  t 

S  A 

M 

•  tot 

II  If 

• 

1  C  to  t 

to*  to* 

to*  a  • 

•• 

•  •  k 

„l  ft  1> 

»• 

ft  O  •* 

uri*a 

m 

•  A 

ti ».  ft  r. 

•  II 

•  M  ft 

ft 

ft  ;•  ii 

•I  •  • 

ft 

OliM 

f  f*ft 

94 

ft  4. 

tin?: 

•  O 

•  ft  .1 

1* 

ft  lit*  • 

M  *  A  * 

«• 

•  v> 

ft  tift  ft 

ft 

ft  ft 

•  ft  ft 

►  ft 

ft  ft  I 

tl 

r%  t  ft 

to‘111 

X 

A  A  A 

a  •  *«  n 

• 

•  • 

•  •  A  ft  •» 

• 

•  •» 

ft 

•  to!  to*  1 

•  *  A  A  • 

!» 

*»  •( 

•• 

lift 

ft  ••  !• 

ft  A 

ft  “  A 

r 

a  •  :> 

•*  t  ••  •« 

t. 

ft  •  A  ft 

If  ft  J 

• 

►  ► 

1  ft 

•  r  ?•  v 

V  ft 

rc.i 

»• 

ft 

r  c>  ri'i  v  n  a 

•' 

r  ►  c 

•  VII  ft 

•9 

ie  r 

<  t 

ft  ft  •• 

• 

•  A 

IJ 

•  I'MZ 

tl  w  o 

«C 

li  •  to*  iS  0 

ft  «J  to*  ft* 

ft  • 

•  • 

•  t  (i 

.i  ft  •  • «  »* 

t  • 

* 

mi 

• 

r  ti  9  t 

tu.  ft*  a 

1  * 

1*  l*«t  II 

•  1  in* 

t*r 

M  • 

t  •  • 

•  • 

•  1  * 

1  f*  ft 

• 

•» 

•  •  *•* 

i  •  •  •» 

•  c.  it  n 

r 

•  t*  1.  li  •  A 

•  1 1  ■»  fr  •  t  to 

ft  •  to  r  %  t*  • 

»  «*  I  ►  C»  v 

»»  I*  N  »*»r  U 

i*M«i  •  |  • 
M  F  h 

is  4  l« 


•»  fcS  »  ‘»*t  I  »r-  l  •  l.  L 

*1  i  n  ft  ••  ft  «j  ft  •  •  •*  •  *  ••  •  «f  •*  *  ^  i.i 

I  »•»  Mif  •  •  I  •  "  *  I 

I  •  »M^I«*U  t*M  »  t 

fc.  *»►  •  «!•«.►»,•»»  M  Vl  M]e  M  I  .1 

to  II  •!  tl  to  to  |l  Hto  »•  to  to  j  »j 

It  I'  1  ••  w  l  li  i  (  Ik  I  Ul  l»l'l  I  (. 

i  ft  ft  £»  MMHO  l»t  •  Ml  I'UMl.ll 


•  •  •*»i'  I' ••  *1  II  h  ||  t  I*  (I 

iNto  M  II  *»  .1"  *»"  II  VI  II  "  Ml 

Mu  ii/  a  M »’  • 

•*  a  a  9  I  •  •*  ••  •  ••  I  •••0*1 

rnKrt  •»’  jMtoM*.  n  •  «••*• 

o  »  o  d  *  •*  r» 

f.  «  i  ••  •>  i  **•'«'»  i 

r  »  •  f*  c.  •  v  n  i  (i  •  &  arc.  i 


I!  E 


O  «Mi 
v»  0  to 


CUIMiU 


V*  tt*  to*  «-  to*  %*  W 


iimi.v.  4v<«  p  r  i  v'l/'i'  nr  ►  t  r  Mfiv»  #v'if  ► 

*ft  m  v*i>  irv  vvi'itM  f.  ft  »  i  i  r  ref  r^<  r  •  rv  o  iiruri.r  n 


3SR-1034 


♦  r.  a 

■>  t  « 

*  z  — 


Li  ->  CJ  t*  l»  M 

N  H  n  *  H 

IT  **  >“  ft* 

r  **•  ► *  ;?  »4 

ft  H  •  *  • 

o  H  Z  H  H 

t  C'  •»  f.) 

•  n  :*  c»  x  » 


a 

O 

c. 

Cl 

n 

H 

w 

M 

O 

o 

a 

►- 

►- 

►* 

o 

o 

O 

•1 

e> 

o 

ft* 

n 

ft. 

** 

ft* 

ft 

a 

H 

H 

H 

*■4 

U. 

M  L. 

ft.  CL 

T  •* 

n  >*  n  k  n  ^  r 

*  V.  *  *  •  2  *  X 

r«  »4  *■<  tl  M  ft* 

•  w  •  w  •  *  —  1 1 

o  o  ?  o  ;»  o  o  »  "  i 

|l  Li  II  LI  II  LI  II  H  • 

—  <  ft  ^  ^  *; 

►>  ~i  T  *5  *■'  I 

*  »  H  •  *>  *  •  I 

w  to  H  **  ►'  H  il  ! 
*  **  L  ft*  L  w  L  ••  w  *1*  1 

IC  H  D  H  >  H  !•*  4  II  < 


ft*  >.  ft.  *.4 

Cl  S?  tl  •- 
O  •  CJ  c. 
I'J  rl  H  * 


->  i.j 

H  >.  *»'•  ** 

X  1-1 

UJ  to 

•  z> 

r4  ft. 

w  r-  ** 

r«  *: 

U 

L>  *t  L- 

• 

L. ( :  Li 

•-  L  > 

1/ 

(i 

s:  *t 

L*  ►  - 

r»  i- 

D 

ft: 

ft.  Ci  s: 

U3  J 

!•*  4  « 

J; 

H 

0 

O  t>  (1 

o-  O 

Ci  r  o 

O 

1. 

t< 

1  1.  O 

l’  i' 

4  r  .1 

t; 

L 

• 

L"-  •  L 

l>  l.i 

H 

• 

• 

M 

tl 

t) 

n  n 

a 

C> 

ij 

r  12 

«•« 

H 

ti  ct 

»4 

i 

? 

ft) 

■>  *)  tl  »«  M  to  .♦  ifnfl  ►  tr  O'  n  w  m  ►.  if  UI  tn  r-  <r  c>  r>  •■<  m  »  v.  ta  n--  C'  r  »<  .•  ►’  '•J  ^ 

F  tJ  »«  H  f«  •■•  »«  •<  «•<  M  r«  rt  t I  fg  f  •  f  g  M  fg  HU  tl  f  I  f>  M  f  f>  ►>  ►">  J  ^  f, 

H«HH  »(  H  Ht'H  Hrt  HHHr  ft  *1  M 


C-10 


StWO'JTTNE  L*P'_*C 

••THIS  SU’OOUTTNE  CCNTOOLS  THE  SOLUTION  CP  POISSON'S  EOTN 

L»“L*CIANtPST>rCONST*NT 


3SR-1034 


09 

09 

H 

90 

«« 

> 

H 

> 

1 

• 

n 

09 

» 

09  ri 

X 

r 

H  — 

z 

Ll 

Z  «!• 

^  ►J 

n 

0 

O  t  * 

• 

O 

•u  n 

M  •• 

a 

L  •  - 

r. 

ft 

Cl  • 

ri  • 

•  09 

1  x 

Lj 

09  09. 

1  7. 

4J 

M  r-4 

•  w 

o 

Z  w 

M  » 

•t 

-  n 

w  • 

X 

M  73 

31  ** 

fl 

Ll  • 

Ll 

♦  M 

c 

CL  *9 

S 3 

Cl 

•J  rl 

*>  -> 

Li  *- 

c» 

•  • 

55 

L  M 

L 

►  •  H 

« 

•  fi 

•*  90 

••  0 

Z 

3  3 

f> 

H  J 

CJ 

♦  ♦ 

M 

CL  L  • 

in 

*»  ^ 

>•  O 

fi 

*3  1 

fl 

y  — 

M 

•  » 

►i 

•  « 

ft 

H  ri 

X 

X  09 

0. 

1  ** 

X 

Z  N 

H  » 

n 

•*  fl 

ft 

W  W 

Li  L 

on  * 

a 

n  • 

f  J  ** 

a  ft. 

L 

♦  in 

•  r 

r»  n 

o  f » 

•*  V 

fl 

n  • 

►  • 

H  • 

r 

1  CJ 

09 

O 

ft  ** 

ft 

*>  i» 

p-4 

r 

>*  H 

i.i 

•  Cr. 

ft 

Li  *J 

7  f 

x 

ri  *f  ft 

X 

•X 

•  • 

1  n>  f» 

z 

«»  M 

i»  x  e 

►- 

M  »  X 

* 

CJ  X 

r  z  *• 

sr  x 

x  w  r» 

*•  ►< 

>4 

n  *• 

i>  CT  M 
ft  U 
J  Z  ? 

O  LI 
i<l  H  „i 
ri  L  >  •* 
73  Y  > 

~j  L* 

i*  >:  ^ 
Z  h  t’ 
►4  C«  Li 


r- 

LJ 

ft- 

r 

r  li 

cr. 

►  4  Li 

o 

C  D 

li. 

O  pJ 

in 

a 

z 

i.i  > 

a 

T 

ftl 

ft 

X  Z 

z 

X 

c* 

• 

l!  M 

fti 

c 

L*  X 

Y. 

z 

U 

• 

M 

fl  ^ 

H 

r  u 

• 

L.  C’ 

Z  L 

ft.  H 

9.J 

m  r 

LJ  Ut 

cr  ct 

O  9f 

fti  O 

o 

C.  Li 

X  Y. 

Cl  X  09 

ft1 

fti  • 

ft  <S  Ll 

It  X 

t>  ft 

ft  Li 

n  li¬ 

r« 

M  Cl 

n  ft  3 

ft' 

N 

ft*  CJ 

H  «T  *i 

X  L- 

ft* 

a  z 

CL  •< 

li  li  9.1 

• 

o  • 

LJ  Z  > 

V  t-  o 

ri 

f*  o 

M 

?rn 

X 

•J 

CJ  z 

li  X 

•* 

O  Li 

*t  z  z  a  **  •  h  •» 

X  •  •  ~  f  4  O'.  *6. 

ft.  M  ft  I  •  «X  Ml' 

1 1  1 1  U  *  •  pi  *-  r* 

OiOL  ri  ;«  I.'  ft  ri 

o  •  i.i  l«  ri  •'  r. 

fi  ci  <j  •  z  »*  ► 

s  r;  t.  n  ->  »•*  h  :: 

Cr  ri  ri  Cr  M  *  ft-  •  •  ►  i>  ri  r*  ( 

t) 

L  O  C  r  ||  «  O  L  L  V*  O  c«  I 

•  (I  G  :«  M  oo  (•  C  t>  H  •  l'  O  1 

•  » 

fi  in 

n  n 


O  M  •  •  I 
I.i  M  M 
•  L<  l»  I* 
5*  M  ^ 
J  ^  H 
Z  ft  b  L* 

o  i.j  ft-  x 


Li  X  O 
C.  -J  ft* 

cl  5  x 
r.i  ft>  i» 

H|JZ 
L  *  O'.  13 
U  L. 

r-  r. 
z  ~j  «* 

31-1. 

:«  o: 
Lj  ► 
3  f*  L 
X  ►  » 

.T  O 
Ll  ft  ►  - 
:»  M 


ft- 

H 


C  ft 

•S  ►  IT: 
L’  O 


ii : 


**  fi 

•N  C>  CJ  u  c 

o  «r  OH  . 

*•  .J  iO  Z  i 

ft>  n  rr  • 

»  .  t"  L  L.  M  I 

C.  L-  •  lit 

ft.  fr  Li  m  J 

X  o  I'.  ft*  i 

X  I’  I  3  (J  I 

OP  ti  * 

.j  i.‘  i*  y  m  i 

*J  c».  »-  n  i 

•i  •*  i.  l.'Li 

li  •  M  •  Cl  I 

•  • 

CJ 

c.* 

fti 


•*  RH  3 


li. 

—  M 

I  fs*  L* 
►  ft. 
M  H 

*  !-« 

I  i.J  M 


i.l  i.J 
Z3  h 


ft.  *t 

•J  ?• 

Li  in 

n  r  u 

•  o  o  •* 
o  h 

*•  L  *  •  • 

ft-  *f.  H  M 
_J  4-  Mm 
M  X 


I.i 


1  ft  <  ft 
II  -J 
I  f .  L- 
i  |i  O 


ft-  3. 

Y.  C 


3  1.  Ci  f 
•■•Ci? 
-I  «'  <r  ft* 

ri  O  il 

•x  3  Ci  fl 
o  •  c.  ft* 


a 

L 

X. 


ft" 

L< 

> 

-I 


M  X 
ft*  c*  . 

•  L- 
H  w 

M  Z  Cli 

•»  *•  r  n 

ft>  cr  ft  *• 

•  o  m  in 
M  il  ►*  M 
w  ~  O  CJ 
XML-ft.fi- 

Li  m  :>  x 

C»  •  Z  X  X 

•  :  ri  ft.  ,| 

M  M  ••  ft  .J 

II  w  ;r  ft.  .j 

t  .  ft*  C  J  91  «1 
M  O  O  4  O 
• 

C) 

1. 


J  ll  o 


O  11  O  U  c.1  o  ll 


I 

Ct  H  N  M  f  Olfi  M  OOCMCl 
Li  H 

fi 

•< 

X 


H  N  cur  >0  Mft  0  Cl  M  N  M  <c  m  X  w  ct  CJ  #4  fti  Ki  ^  in  in  X  ft>  cr«  fj  H  ftl  M  in  U»  ©  CM  C)  *4  fti  M 

n  M  H  H  "  f  I  fl  H  H  ftl  Cl  N  M  M  N  N  M  M  N  M  M  •!  M  M  H  M  I-'  M  M  tf  it  a  it  it  it  it  it  it  it  V>  u*  W*  in 


C-ll 


SEV  STPrAKPUNCTTON  V SLUES  TK  THE  •»•••»  SPP»T 


c  f T  *  J) =0°  f I  *  J  J 


U 

| 


3SR-1034 


!■ 

j 


in  U>  f-  »  r  o  H  N  h  ff  IfllO  h  »  (T  L)  »<  N  fn  <f 
i/t  4P  ui  u>  m  u?  io  is  is  uj  10  iu  io  io  is  f*  r-  f-  ^  r» 


Cr 

Cl 


C-12 


zucum 

CTfCNriCN  3(KX.KYPlt.RU«tN;>.VXI2> 

"9v?'»*l?»icc  fO.sn.iRtmii.uuiii  i.ivxtii .vvtin 
•  infers  clcbal 


V» 

U 

3 


O 

Li 

v« 

u 

H 

C» 

Ui 

f. 

V* 

n 

z 


v> 

u 

o 

H 

o 

z 

Vf 

ui 

ti 

o 


~  x 

►»  H 
•*  L 
X  VI 
H  VI 
U  CL 
VI  • 

VI  «• 

0.  IJ* 

«  M 

Pi  w 

•  X 

M  UI 

**  r 
x  z 
UI  vi 
O'* 

5?  w 
M  vl 

•  •  t) 
-*  V<  UI 
M  w 

•  X  Cl 
M  VI 

•»  r» 
x  ?o 

Ul  V*  o 


x 

H 


X 

Vi 

cc 


«c 

a: 


UJ 

K 


2 

►* 

fc‘ 

O 


X 

6? 

V- 


2 

Ui 

►J 

u 

2 

V* 


X 

> 

* 

a: 

D 

ft: 


pj 

¥* 


1+ 

*n 


CL 

X 

2 


•* 

X 


X 

0 


Vi 

X 


X 

* 


X 

CJ 


x 

x 

< 

CD 

U 

o 

O 

z 


-J 

V 

X 

V- 

v> 

(Tj 

CJ 

X 


X  ft* 
v*  -> 
•  11 
ID 

H  T. 
•  C> 

x  r 
.  V.  U 


•  x  n  h  o 


VI 

ftl 

o 


2 

i* 

* 

z 

X 

Vi 

M 

O  CI 

-J 

ft 

VI 

ft 

n 

o 

z 

M 

ft* 

(.1 

VI 

z 

a 

z 

ftl  VI 

-I 

X 

1- 

ft 

Vt 

X 

vt 

UJ 

Vi 

vi 

n 

m 

ft 

»* 

u. 

ft 

u-  V 

> 

Ui 

0 

ft 

o 

01 

Pi 

m  mi 

ftl 

•i 

o 

o 

** 

O 

> 

V* 

0.  N 

ft 

:: 

ft 

o 

**  ft.  1 

r  o 

n 

O 

u 

z 

*. 

ri  Pi 

O 

h- 

? 

z 

H  S 

V 

Ui 

H 

ft 

t* 

z 

rl  pi  ft 

c>  *t 

Cl 

M 

z 

z 

o 

ft  v 

» 

UI 

VI 

ft 

ft  v 

V*. 

M 

>• 

ft 

ft 

v  ft  H 

t'  V 

l 

t) 

M 

fl 

u 

X  N. 

r< 

v . 

Cl 

t 

V-  ft* 

c 

ft 

0 

St 

ft 

UJ 

H*)  •)*) 

p:  Vi 

ti 

f  1 

• 

II 

Vi 

ii 

o 

•  I 

ui 

CJ 

V-  .1 

r 

X 

!■ 

Vi 

ft 

is. 

ti 

•-  ft-  ft- 

•  ui 

V4 

• 

o 

c.1 

ft 

X 

V 

X  X 

ft. 

o 

I.i 

a: 

♦  or 

l' 

Pi  M 

r.’ 

UI 

CJ 

ft 

H 

*) 

X  x  Cl  III 

«r»  o 

• 

10 

UJ 

L) 

ft. 

I.I  V 

Pj 

>■ 

.1 

v- 

*1  -1 

V* 

ft-  w 

r* 

•1 

ft 

V 

uj  ui  ti  r> 

Ui 

• 

►J 

CJ 

•1 

f  o 

►  • 

LC 

Ii 

•t 

►  V 

» 

V 

ft* 

r.. 

ft 

D 

ti 

P  r.  j  z  z 

l  i  (  i 

►  1  *- 

(•1 

a 

t  * 

ci 

s*  z 

Vi 

> 

ft* 

> 

•f  «( 

y. 

•l 

1' 

ft 

CJ 

ri 

it  V  *>  V  CJT 

V-  •« 

♦  c:. 

V- 

Vi 

V- 

;• 

l-l 

-J 

M  »•» 

M 

-1 

•J 

V 

Cr 

*<  mi 

13 

ft" 

r  >. 

C' 

a. 

ft- 

0 

ft 

Cr 

rl 

vi  h  •  v-  ra 

t!  VLiHK 

M 

w 

V 

..J 

l«  »• 

•  ' 

•  I 

»J 

i» 

O 

ftJ 

V 

o 

r  <t 

cr 

ft 

o 

ft 

Cl. 

II  II  y  «f  o 

ti  ft  C.* 

1'  i.i 

ft 

u 

ft 

ci 

n 

fl  N 

•i 

•i 

C  i 

1. 

X  «f 

L 

L 

f»  Ci 

Ci 

U 

V 

» 

o 

VI  p:  ft-  tl  u  2 

h  .t  • 

H  U 

:« 

VI 

It- 

• 

o 

U  VI  Vi  CJ 

u 

1 1 

u 

ft 

>  l* 

Cr 

ft 

Ii  u. 

u. 

ft 

V. 

ft 

V 

n 

*>  *»  L>  t»  ft.  ill 

• 

ft 

ft 

• 

v< 

ft 

ft 

O 

CJ 

C) 

p  n 

o 

o 

ft 

Cl 

*4 

in 

1.) 

r  Cl 

VI 

Ui 

ft 

CJ 

H 

t>  tl 

Cl 

N 

ft 

H 

H  Pi 

M 

ft 

PI 

V* 

• 

v» 

w 

IU  I 

(V 

ft 

r 

p 

u 


H,#U  hit  IflUlVpPllHNV  if  Ui  ifi  P*  if.  CP  f"?  H  tj  M  it  UI  |,*»  P  (p  (T.  «1  rt  M  Ki  .1  W  If  f*  ft'  01  n  ♦*«  M  M  if  IP  IP  ft*  «T  01  C  i  H  N 
H  H  »<  H  H  H  »4  »<  H  »l  fti  N  f  i  f  J  1 1  M  f  J  M  f  J  h  M  Ki  h  M  M  h  11  h  M  tf  it  if  if  St  if  if  ip  :f  »  UI  Ui  u> 


N 

a 

Ui 

o 

«r 

o 

« 

Ll 


01 

M 


VI 

v> 

Ui 

a 

ft. 

X 

o 

v 


¥• 

ft. 

Cl 

til 


C- 13 


DTfftN'ICM  XXtNX*NTI 


m 

m 

#• 

m 

n 

m 

• 

n 

H 

• 

w 

X 

X 

s? 

X 

1 

l/t 

X 

** 

z 

X 

n 

o 

1 

• 

M 

fi 

H 

n 

4 

H 

• 

M 

C 

H 

M 

•• 

z 

X 

•> 

o 

X 

X 

• 

u 

X 

X 

«« 

4  *• 

z 

► 

4 

>  3 

w 

ft* 

> 

C  * 

o 

«r 

O 

t)  H 

X 

a 

u 

►  4  ~ 

1C 

• 

CJ  o  O 

V 

b 

13 

a  •  m 

o 

• 

•  Cr 

CJ 

m 

X  M 

Z  1  s 

> 

n  i 

H 

t’ 

O'  ^  • 

•>  n  ^ 

4^ 

H 

#  n 

•  *  Cl 

u. 

•J 

**  • 

Hl<  Z 

1 

C> 

•*  X 

♦  -  1^ 

3 

X 

*)  z 

►  i  X  1 

U 

o 

•  M 

cxo 

z 

n  x 

o.  ♦  -J 

u. 

III 

•*  X 

1  «-  o 

•• 

t'lfl  ♦ 

1  li 

♦ 

3 

fl  •  « 

*  n 

♦  »<  ♦ 

H  • 

3  *  f 

• 

X  ♦  *<  H  •  H  1 

►  ><*)•-  r  e«  **  • 

I  Z  Z  •  X  X  ♦  X  I" 

|  X  »  II  H  X  Z  M  *  ~ 

n  h  >  •  *-  m  o  *  —  •-  x 

|,i  H  ll  6  •  J  ♦<  ft.  •  X 

/  fi  *>  r j  **  i*  o  (i  w  r*  ii 

o  •  ;►  *  o  L.  H  *  a  ^ 

U)  f  .j  rf  fj  t  '  n  UH  •> 

o  *  <  i  iJ  ii  lit  u*  « •  H  • 

lif  &  tt  h  n  n  >  t-  H  >  I  i 

X  ii  II  :»  »  J  X  >  I  I  w 

v-  o  n  h  c>  c  iio  i-  •'  > 
•  i.i  t  ,  •  Ci  ii  ii  r.  u  I*  x 


x 
X  x 

»  :*  *» 

1,1  CJ  (l  ~  ll1 

i-  :>  1 1  *i  r> 
ii  »:  •  z 

II  ►«  i.*  x  *-« 
n  I-  i*  X 
►J  Z  *-  z 
PU|‘  >  c 
ll  O  •  x  t! 
• 

IA  O 

t*  C» 


o  tt 

t>  111 

x  •  * 

ODI*. 

Cir  •  I 

•  N  fl 

~  I  .1 

*•  **o 

M  “>  U 

**-•) 

H  X  • 

♦  X  ►  » 

Im  4  •- 

x  aK 

0  cl  X 

I  •  • 

»*  T  »* 

H 

*  t  H  • 

X  ►l-M 

?  Z  •  X  *• 

>.  *  •  1 1  K  O 

fl  H  ».«*-»»  X 

i-  ii  n  o.  •  o*  r« 

OI«t)  ~i  •*  >  "  C:  ! 

•  •  *  O  "  i*  ; 

I,,  c .1  t,  «r r  U  «  •  -»  ll  < 

•  I  ,  ll  ►  IJ  It  •  HI 

r>  f»j  r*  h  ii  ?»  n  <v  I 

i<  j  -i  ; 

o  c>  o  €j  x  c.*  i 

C»  Ct  ll  C»  t'  I*  >•  I  i 


i  •  t  *  o 
»>!•»•  3 
*-  jr  i  -  n 
X  O  t  V 

X  «  »  ill 


»  «r  V>  U>  X  •»  «'  tl  e4  fe  *  *  IA  U»  X  «n  O  •!  M  *t  in  10  X  •>  9*  CJ  «*’  *J  M  •» '10  tf.  X  ■*  £ 

hhhhhhhhm  »i^n  f*J  m  t-i  r  i  m  ti  im  m  ►  h^HXMXhx 


C- 14 


5®»T 


3SR-1034 


I 


i 


> 


N 


>  b 

•  t>  ►  £ 

*IM»*  ft  ► 

*•  1>U  Hd 
««.*«>* 

»'U 

•  i.»  ft.  i.*  ft  ft 

Mft  * 

r#  .  i  *h. 

*•  D  t*  ft  ft  ft  »• 

b  s  X  y  a  f  * 

K 

ft 

NllAMMId 

ft  »•  fi  n  »« t§  t« 


ft 

•4 

!• 

ft* 

U 


fi* 


Wft* 

x  m 


tl* 


> 

it 

• 

ft* 

f’ 

»« 


► 

r 

p 

•* 

£ 


b 

* 


•  * 


V 

M 


M  • 

•  ft 

ft'  r 

ft  • 

•  ft 

14  V 

•  %• 

i*5 

ft  If  «" 


ft* 

V 


..r: 

•4  ••  ft  O 

•  •  r  x 

M  #•  •  ft: 

ft  ft  M  • 

•  •  ft  • 

It  II  ft  t« 

•  *  if  • 

•  ft  ft  ft 

mi  5»f» 


f'J*  *t 
•*  ft  .< 


,  •>  ft 
till  itll 
•  •  »#  • 

SH  44 
•  •  *JX 

•  •-  «  «f  ft  a?  • 
ft  ft  «»•«  ft  ft  Oft 
aift  •  r 

n  •  ••nr  *  m>  ft 
•  ***««*•  ft  It  ft. 
W',  !</  •  ift  *!  f 

•  •  *1  l«  II  »  w  •  » 

m  •  ft  ft  •  ft  #*<••£.  ft  t»tr 
ft .  *  *  fe*  *  -■  S# 


ft  • 


a  •••••••  lu* 

ft  A  MH  •(•<MltS  t*  *1  » 

•  I  ft  II  ft  ft  ft  ft  ft  ft  ft  ft 

ur  «*i*fti>ftftft  »  6  •  t 

Oftft  •  •  •  •  •  • 

%  >**••••  ftolf^l 

ftft  •  ftl  «tft  Ml  4*1 

X  I  *  ft  C 

ft  (tftftftf  ftti  •  ft.lft  b 
►  lIMSl  J*  MM  I  •  Ml 

9m  ftiM  V  M- 

o  ;t  r  •  «••*««»  . 

ar  ft  ft*  *  «i  ft  •«  ft  ft  t,  •«  t«  *#•  ft  f  +  £*  ft  ,*  «#♦>  *4 

t  Ml  I  ft  I  I  «•*•*>  «*  *S  ft  «>  ft 

ft-  S»  !•  In  ft  ••  ft  «t  ft  ft  ft  *  ^  ftl  dt  l,  f  lift*  »  |  »J  I  ft» 
JlftMMM  l.MII  fiMIlUftftf  fit  ftM-ftfeftb 

*  ••  ».  *•  r  •»  •* 


pi  t*  w 

’!>»*®i*|ftr  ail  • 


fc  * 

If 


4* 


•  « 


c*ii 


j^o  Tmrt-iM 


3SK-J054 


(C*lf 


3SU-1034 


**  *  C 

V  £  P 

c  •  •» 

•  * 

#  » 

IS  ie* 

»*•*  »  «;  I 

sr  v‘=vi 

irfti j 

i  *  i 
;  t  •«*.  • 


•*  #ft  • 

•  r#*  • 
I  *  #/ 

»«*•  Mf* 

•  »*»  r 
«m  • ■  m  ft - 

•  •« 

tJ  •  •» 


fl.*£  i£S 

•  i#»  I 

•  *;»»•>  i 

«  •  «» »•  •  »)« 
ft  # »  *  l  i 

•*  •  «•$  •*•  i 


°  S  e 

•  m  it 

**  «i 

1 1*  ?  « 
V  •  «*  ** 

j.r  ^  . 

<<  •  ?  m 

*  t  »•  * 

»  t.fc  « 
»  ft  If  ft 

*>  ••  «  i$ 

b  •  •  *w- 

F  t*  a  •»  •  • 

ffft»f*  i»*f 

•H  •:  *i! 

.jint,,! 

\»  *3  4  •  *"  £  y 

9  *  **$-  '$  *a  \ 

I  •  ft  *j  »  •»*>  < 

I  ««t  •*  ft  «►  1 

I  I  ft-  ii'1 6  **%.  * 
>  §»•»#->»*•« 4«1 

’v  s? 


F  . 
k!!  k 

+  —  i 


.••ft-  f 

.r  #&K 
»usrr 
!ir.fc*s 

«*#«*=  t» 
’  -9  «>  •» 

4 •  til  •• 


«*»*#  fc> 


•*#»«!(•  I#  *  *.•?*. 


£»iststssiisrns*issrr.snrsi 


«*»  f 


C-*« 


3SR-J034 


m  — 

i  !i| 

m  V  »*%# 


sn 


I 

li  i 

*-■*•*  i 

*y  •-  a*  i 

•  Si?  < 

**  #•*  m  J 

r|*  «•*  < 
**«»«  <| 
•t  4*  2«  J 

#  •  •'  #*  I 
#  •  %•  » 


*£ 

tui 

«P| 

rtf 

••  #  * 

«  <#  •'#  ♦• 
*  •  *-* 

S*$J4  J 


,  =  t 

}§. 

iff 

tit  i 

\  [i  ’ 

*>!»** 


*  »  •  I  t  i  * 

■  -  f  -  1  \  *  *» 

#“  T  *■  *-  *  *  *  * 

*•***'  "  ■  l 

m  m>  h  •“  ■  ..  ifa  .«..  •» 

hrrrHii 

UUiii!i 


S  I. 

f  Ji- 

f  •  f 
i  ii  f 

E  i>  » 

i  M  3 

I  ‘t  t 


f  f 


•  •  •  v*4  i * 

ijfllfl 

%  §  fe  ♦  ♦  ♦ 
#5  **’  w*~-  »V 

«u)  t  •  «  at?  *  *.- 

*  •*  *  -%  « *.  « 
&&*£}*'*  ^  %  % 


£  tUii 


•*t«  »>  4 

■■*■  ■» 


c*« 


3SR-I054 


3SR-1034 


: 

m 

> 

n  • 

• 

9* 

• 

at 

m  r 

«•  •  9t 

• 

• 

• 

*S! 

•  jl 

n 

«•  • 

•  m 

Vi 

•  C* 
»l  M  r+ 

sap.:; 

Z 

••  r 

•  mm 

mm 

►  1* 

• 

m 

•  «» 

1  « 

mm 

r  m 

••  •  •» 

• 

# 

gr  •• 

•  M 

mm 

•  #t 

»  M  • 

• 

•  • 

r  t* 

H  W 

m  m 

*•1* 

1*.  • 

• 

•  •* 

*■  • 

•  H 

•  «• 

•*  • 

•  #t  ► 

•  r 

»  • 

Ml 

•«  •• 

V9  • 

§  v  vt* 

•  V 

•  •• 

*-  • 

•  M 

•  •« 

m  t>n 

•  «• 

•m 

«s- 

»  • 

»'l» 

Si  Mfl 

•  •• 

•*  • 

•  ** 

*•  • 

•«  i 

•  *9  • 

tl*0 

•  an 

Jt  • 

•  • 

»  • 

^  • 

t  / 

»  •  • 

•  •» 

••  ft 

•IM 

•  It 

►  • 

•••«  ► 

at  • 

p  r 

•  •» 

mm  • 

r  •• 

•  ►  r 

S  Sf» 

►  •» 

It  M 

fill 

•  •• 

fi  illl'A 

•  tin 

•  •» 

»•  •• 

*•  1* 

•  •» 

m  t» 

•»  |*  *•  I* 

•  #v 

f|M 

•  %• 

••  • 

••  V 

»  •If* 

•  ill#  9 

■»  %* 

fiSi  i 

•  m 

*>  %* 

l«  3>  • 

•  •-» 

•  l»n  u 

«5l  • 

m  m  m 

;•  • 

w  •  ~ 

"Hit 

::  * 

»•*  •  » 

mm  • 

vi*  •  •  i* 
•  *1 
I  ff 
i  * 
mm  •  r- % 
l>«  •  J  •• 


It  «i  M 

Ml 

r  »»  *t* 

•  ••»!* 
•  •**««*»> 
•*  •  >  *  t  • 

•  4,  • 

M*  ••  • 

V*  i*  •«.  •  M 

t  I  I  M 

t  r  »*  • 


*  i  « 

^  ^ 

•  •  ••  * 

«•*»  i«  *  i 

•  •  > 

«-  tt  •  r* 

t  m  n  «t  • 
*  * 


•  (l  «  ►» 

•  i  •  n  •*  • 

•  III'  •  M 

•  •  •»  #»•* 

•»•«!•  «  • 

*•  *'■  I  !•  %» 

*>«t  M  •* 

l.f  « I  « 

*• 

•  if-  •••«»■ 
It  l«  L"i"  It  i* 
t-  *»  •  •  • 

*  •  I  •»  I  I 

It  «  *t  M  •»  t 

•  *  •*  b-  t* 

I*  •  i  ••  •  t* 

•  »»«  ••  <r*  it  • 


►>  o  V  if  *•»■»• 

•  •  ?  • 


•  §•  •  •«  •  •  • 

t,  •«,  •>*•••  • 

W  f  •  M  •  I 

^  Vi  ;•«•••«  •  i 

•*•*;»►  • 

iff.  f.i  •' 

#  f*  *•  #•  ?*  M  M  • 

«•'«•*«  •  !•!*• 

i  tv  tf  vii  »  *  •*  ;•  ,i  • 

•  i  k  k  i  *»  *i.  i*  • 

•  •  «  m  i  • 

*1  t«  it  i  it  it  m  •  •  • 

4*  I  >  ^  • 

%,  >  •  «»Jtl  M  M»*  •*  • 


*•  • 
r  %i 
.*»  • 
n  o 
•  •  • 

» «>  •  n 

•  I  M  • 

>  •«  *»  M 

mm  •» 

I'ltl*  II 
>«••*>  t- 
"  •*  • 

i  •  •  •  n 

|l  M  • 

!*•»•»  i 
■  ••'i 
•»  i 

•  •  «.  •  k 

•  v 

i  •  1 1«  it 

.  |i  -  *.  r.  *i 

I 

•  *l  n  •  I  it  i 

w  I  »•  w 

i  k  »•  “  •• 

Ml  t»  • 


««•»  nun 


%••• 


3SR-1034 


L' 

r 


& 


o 

r 


o 

T 

or 


D 

f* 

r 

o 


a. 

r 

LI 

K 

p» 

■)  •• 

• 


X 

2 

M 


rl 
It  • 


rl 

•  I. 
>• 
tii! 
*>  CJ  * 

•  Hi 
► .  ft* 
**  O  > 
>  fr-  I 
I 


a 

ft 

ft 

H 

ft 

♦ 

*. 

O  H 

•  • 

n 

n 

C>  fl 

♦ 

• 

r* 

CJ 

0 

♦ 

*1 

N 

M 

H 

l: 

M 

1 

N 

(1  ft 

•  • 

M 

M 

CSi 

o 

M 

■> 

• 

►  i 

►  1 

►  I 

■» 

*» 

ft* 

f 

• 

►i 

• 

• 

>- 

•  • 

«« 

• 

♦ 

Pi 

ft 

• 

ft 

ft 

• 

> 

M 

rt 

1 

mi  z 

H 

•J 

r- 

•  • 

f 

*• 

■> 

•• 

ft 

»■! 

•4 

• 

-l 

** 

fti 

*-4 

r. 

Sr 

e> 

• 

L 

M 

n 

X 

ft. 

• 

>* 

• 

•  • 

V* 

CJ 

ft 

»- 

•r 

ft 

4 

0m 

• 

L 1 

o 

X 

ft 

♦ 

• 

• 

T 

M 

>  L. 

Lt 

M 

7 

►  « 

f« 

H 

H 

•  • 

Cj 

L’' 

ft 

ff 

e. 

> 

ft 

1 1 

• 

> 

IA 

•t 

►I 

frt 

H 

f  i  fr 

ft* 

M  fr* 

• 

w 

ft- 

ff 

•  t 

f 

• 

f'. 

V. 

M 

H 

• 

X 

• 

ft* 

to 

It 

It 

H 

H 

Z 

C. 

• 

• 

;■» 

CJ 

> 

• 

r* 

„f 

ft* 

•  • 

r- 

c. 

f. 

f  X 

• 

r* 

I . 

ft 

ft 

w 

• 

1 3 

• 

3  « 

rt 

•» 

M 

V- 

L* 

Z 

f» 

OS 

O 

It 

Cl 

ft 

I.* 

ii 

i. 

9. 

•  • 

ii 

tt 

v 

II 

L 

n 

H 

r» 

s 

It 

Vi 

>: 

• 

z 

H 

IM 

a 

II 

a 

it 

■> 

II 

D 

f) 

L* 

i 

Li 

H 

L< 

*W 

H 

z 

*» 

> 

;> 

•  • 

ft 

»_ 

•• 

it 

♦ 

H 

»» 

ft 

3 

• 

o 

*• 

if 

*• 

f4 

( 

Cj 

• 

7  ST 

ft 

U 

X 

• 

♦ 

♦ 

c . 

f 

C 

• 

ft 

c  V 

•  • 

►  • 

f 

c 

t  f* 

1 

CJ 

f. 

r 

if. 

• 

M 

> 

♦ 

• 

1. 

K  tj 

(1 

i< 

»t 

-) 

•  •  ft* 

►• 

II 

ft 

r 

»■  ■ 

•  • 

•  4 

ft 

f* 

f  ♦ 

* < 

n  rt 

•  • 

»~ 

c 

Li 

fl 

r4 

w 

o 

C' 

ft 

it 

ft 

• 

fr¬ 

O 

f’ 

M 

ft 

It 

<>  H 

M 

N 

** 

• 

K  » 

I  t 

-1 

H 

V 

H 

fi 

It 

fr-' 

w 

f  • 

V 

ft* 

•J 

k< 

K 

•J  U 

•  4 

Li 

►• 

a 

w 

ft 

►* 

ft 

•- 

f- 

f- 

U 

• 

it 

►  • 

V 

II 

fr  4 

ft* 

S  ft* 

a 

ft 

7  ►* 

57 

o 

> 

•  1 

w 

fl 

ft* 

s 

> 

>. 

*•  >: 

l.t 

L  fr* 

•  • 

if 

w 

>. 

w 

» 

*-• 

>: 

LI  I'*  f- 

• 

it 

:t  n 

k- 

> 

>: 

C 

C> 

** 

O  L 

i; 

s* 

O 

r. 

> 

L 

ft 

L 

<: 

f 

II 

i 

4 

•  •  l 

•  • 

t> 

r> 

**i 

L 

U 

c* 

1 1 

l 

ft 

L 

• 

D 

L. 

c 

ft 

L 

13 

O 

ft  r> 

rt 

O  ft 

fr- 

II  ft* 

1 

f  » 

ft 

• 

Mlt  WVt 

n 

<•« 

• 

t. 

ft 

c- 

• 

f4  ft 

•  • 

ft 

«•• 

• 

r 

1  « 

f' 

« 

• 

r-. 

• 

ft 

ft 

• 

It 

Ill 

• 

M 

H 

C; 

1“ 

•  r 

• 

• 

• 

•  * 

• 

H 

• 

• 

• 

• 

• 

a 

n 

a 

CJ 

•  • 

o 

n 

n 

• 

CJ 

rj 

n 

*4 

L' 

C) 

•  • 

fl 

in 

r1 

• 

fl 

it 

N 

fl 

•  • 

• 

L> 

L» 

u 

u 

U 

Li  U 

iJ 

o 

L3 

if 

O 

u 

o 

c 

s 


»*  V 

n  * 


M 

o! 


♦ 
H 
** 
a 
.  i 
** 

m 

ft* 

ft. 

♦ 

n 


t 

v 

* 


O  M 

1  C.  r 


m 


M  w 
H  D 
M  » 


S 

I 


»  ft  o  r«  N  «  *  ifl  ui  n  »  r  n  h  n  >’  4  in  «'  i*  o)  p.  t,  m  n  ^  ^  i>*  ^  i'  h  fi  n  *  ir>  i*1  r*  «  ^  ti  h  Kitf  if  «  ^  ^  tj  h 

tl  f)  MHMMH  H  ft  H  ft  H  It  t:  t*  M  N  f|  M  ft  (>j  f;  l'i  M  h  h  h  ►  fl  f*  h‘  M  if  i  it  it  it  a  t/  it  if  ;f  It  i  U>  4A  If  If  l/''  Li  111  Ii3  if  L>  IP 

M  M  H  ft  M  HH  H  ft  HltftHfl  H  H  H  H  H  *4  M  rf  H  ft  M  ft  M  H  M  H  ft  H  f  I  fl  ft  H  H  *|  M  ft  M  f  I  fl  fl  H  M  fl  1 1  t  I  rl  M  H  M  ft 


zC-25 


16*  COf'Ciri=-XT»*Pl«U<I»l»  31  •UCT*1»2J»U<T  •1.111 

1SS  30  TO  130 

i6G  ?o  cuMcrrri=xTPAPitucT*i.Mrr2*.urT*i»KYmi«oci*i»MY*i 

icl  p  yori 


3SR-1034 


C-  26 


3SR-1034 


•» 

p 

m 

a 

i 

* 

7 


5  S 

•  O 

SB  - 

•  C*  H 

»•*  * 

•»  ft  H 

O  •  »P 

•  y  c?*> 

n#  •  •  • . 

•  • 

M  Cl  M  *»klW 

*•  •  f  J  • 

n^C'i  c» 

T  I*  wk 

«*■  *  9 

♦  •>»**  e* 

•  Mr  »cr 

I  H*«  |  3  m  HO  •* 

»  IJ  R  6  *1  14  M  I  •  *• 

.  V  K  •  r.  H  ft. 

.  MJ*  y  ^  ► 

b  •  •«>••:*  cj  m  x 

►.  »«  HU  X 

•  X  »•  l»  T  •  •  «»  t  t 

i  n  ►•  «  <*  r  MU' 

.  V.i  M  ft*  •••  4*» 

•  ic;  •  • 
j  ariir  m 

■  U  |  w  *.  >;  •«  w 

>  if  r>  ft  r  t.  i  i*  ft  ». 
♦  *  r  f.»  r»  ♦•  ••  c-  •  « 


ii  t» 

•  •  • 

#  i  «*  ftj 
•i  i*i  ii  s» 
r»  •  »*  if  Z 
*1  *)  M  O’ 

•  •  <  M*  3 

i!  H  (> 
»-*  L  *•  iV  •  i! 
3»*  i»  lift  fti 


w  •«  f  <  •*#  <#  ft»  «>  r*  S'  +  c  •  •«  *  *  n  *  ft'  » 

•  #«*«••«  ««  M  H 

V 


C-27 


3SK-J054 


m 

l 

I 


I 

u 


4  hi 

r:r. 


s:1 


r.|:r. 


Ikw 

& 

m  i* 

;r.  * 

1.1.  ? 

*  »  •  -  .. 
•  *  • 

•*Ih 

t  Jll  h|  »> 

i  r« »*mm 

s*r  •  •  *  r 

-  KM *#«■!* 


ikp 


** 

#1 


•*  flk 

t« 

•  • 

M  S 

W  2 


M 

VI 

•  ► 

If 

nk 

;t 

4*m|» 

•»#  f 
•  *  W 

h  I, 

VZ  V 

s>% 

!: 


•I 

•• 

M 

IT 

» 

M 

r* 

# 

; 

u 


u  in# 


iM|«n#«I|f4  #f  MU  h 
MMMH 


f. 


C-2* 


